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ABSTRACT 


The data obtained from the flight of a parachute-born 
blunt probe through the D-region of the ionosphere had long 
been reduced by hand. This process was tedious and time- 
consuming. The objective of this work was to study the 
feasibility of performing computerized data reduction 
techniques using the Hybrid Computer Facility at The 
Pennsylvania State University. 

A presentation of the theory of blunt probe operation 
is included with emphasis on the equations necessary to 
perform the analysis. This is followed by a discussion of 
computer program development. Included in this discussion 
is a comparison of computer and hand reduction results for 
the blunt probe launched on 31 January 1972. The comparison 
showed that it was both feasible and desirable to use the 
computer for data reduction. 

The results of computer data reduction performed on 
flight data acquired from five blunt probes are presented 
in the fourth chapter. Four of these probes were launched 
at Wallops Island, Virginia, and the fifth from White Sands 
Missile Range, New Mexico. 

The computer programs were developed in a manner that 
would insure ease of operation and flexibility. Operation 
of the programs is not complicated. Appendix A is the 
operations manual which provides the user with a means of 
becoming familiar with computerized data reduction tech- 
niques . 
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CHAPTER I 


INTRODUCTION 

1 . 1 Historical Background 

In 1964, a program was initiated in which the primary 
objective was to develop a subsonic probe which would 
measure the characteristics of the D-region of the 
Ionosphere (L.C. Hale, 1969). Early probes were launched 
using the standard ARCAS meteorological rockets. At about 
80Km the probe separated from the rocket and fell, in the 
inverted position, suspended from a parachute. Accuracy of 
measurements required the probe to remain subsonic through- 
out the entire range of data acquisition, 40 to 80Km. 

The theory of parachute born blunt probes indicated 
that the probes collected charge according to a simple 
mobility mechanism (D„P. Hoult, 1965 and A.S. Sonin, 1967). 
Because of the simplicity of the theory involved the blunt 
probe collected data which were independent of the angle 
of attack, with current proportional to probe potential. 
Since those early years the blunt probe has gone through 
several design modifications. Figure 1 shows the present 
blunt probe used for data acquisition. The collector sur- 
face is the small inner disc at the top of the probe, and 
it in turn is surrounded by the guard ring. 

While the blunt probe was undergoing improvements in 
design the theory of operation was continually being 
studied. The theory of the blunt probe was originally 
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developed by D.P. Hoult in 1965. In 1968, L„C. Hale et al., 
gave a summary of the theory associated with the probe and 
in more recent years the application of probe theory was 
performed by J.D. Mitchell, 1973. 

1 . 2 Previously Related Work 

Development of the theory of operation was only the 
beginning of the efforts to gain an understanding of the 
D-region. The theory was only a tool which had to be 
applied to the data acquired during the flight of a probe. 
Each of the investigators mentioned in the previous section 
were faced with the task of data reduction. In order to 
obtain results from the probe flight, the data had to be 
hand scaled several times, averaged, and tabulated. This 
process took approximately two weeks to complete, and then 
the tabulated data were used in the appropriate equations 
for determining positive ion and negative conductivities. 

The final phase of data reduction consisted of searching 
numerous tables and applying the same formulae innumerable 
times. After approximately four weeks of intensive 
calculations the end product of positive ion, negative ion 
and electron densities was obtained. 

1 . 3 Specific Statement of the Problem 

Although there have been continual improvements in the 
blunt probe and its associated theory, the method of data 
reduction has remained somewhat stagnant. The methods used 
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have proven to be both time consuming and laborious. They 
are time consuming from the standpoint of taking three to 
four weeks to reduce data from a rocket launch, and 
laborious because of the repetitive nature of the analysis. 
There is a certain amount of error introduced because of 
these two factors. Errors introduced in the scaling phase 
of the analysis arise because no two people scale the data 
in the exact same way. These errors, coupled with those 
introduced m the repetitive applications of the formulae 
reduce the accuracy of conductivities, ion densities and 
electron densities. These effects could be decreased 
considerably if the amount of human interaction with the 
data could be reduced. One obvious way to achieve this 
reduction of human involvement is through the use of 
computer aided data reduction. 

The purpose of this work is to examine the feasibility 
of computer aided data analysis and to develop such a 
system. To achieve these objectives a series of computer 
programs are developed. These programs i) digitize blunt 
probe data, ii) analyze flight slopes, iii) calculate 
positive and negative conductivities, and iv) calculate 
positive ion, negative ion and electron densities. Before 
these programs could be developed a decision was made on 
the type of computer facility to be used. 

There were two possible choices of computer systems 
available for this study: the digital computer or hybrid 

computer facilities of The Pennsylvania State University. 
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Hybrid computer analysis has a definite advantage over the 
purely digital computer approach. The analog portion of the 
hybrid system allows direct input of the data stored on 
quarter inch magnetic tape, and also facilitates implementa- 
tion of signal sampling procedures and filtering techniques. 
In effect, the analog computer takes the place of the strip 
chart recording and provides the necessary interface between 
raw data and the digital computer phase of data reduction. 
The digital computer then applies the known mathematical 
relationships to the information received from the analog 
computer . 

Using a hybrid system for reduction of data obtained 
from a rocket born probe has definite advantages over the 
present hand reduction system; i) appreciably reduces the 
time between data acquisition and results; ii) frees 
personnel from mathematical computations which can be better 
performed by a computer; and iii) provides a greater 
accuracy in results. 



CHAPTER II 


BLUNT PROBE THEORY 

2 . 1 Introduction 

The theory of operation of the blunt probe can be 

described by the voltage applied to the collector surface. 

This applied voltage varies between -10 and +10 volts, 

Figure 2(a) . By using either a calibration or flight plug 

two modes of operation can be observed as the voltage varies 

through its range. One of these modes of operation is the 

calibration mode, Figure 2(b) . The second is the flight 

mode, Figure 2(c), which is an example of data obtained at 

one point in the flight of a probe. Data obtained from 

either of these two modes of operation represents the 

current of the collector surface (I ) as a function of 

c 

time. The ratio of current obtained xn the flight mode to 
that obtained in the calibration mode is directly pro- 
portional to conductivity. 

Presently there is some consideration being given to 
the possibility of reducing the probe sweep voltage to a 
lower range of -5 to +5 volts . The nature of the data 
reduction process insures that such a change would have 
little effect on the theory of blunt probe operation if it 
were accounted for in the calculation of the probe E-field. 

The theoretical analysis of data acquired by the blunt 
probe was designed to obtain knowledge about the chemical 
composition of the D-region of the Ionosphere. In order to 
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Figure 2: Blunt Probe Sweep Characteristics 
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gain this knowledge it is necessary to extract positive ion 
and negative conductivities and densities from the data. 
Blunt probe theory dictates that the conductivities be 
obtained from the recorded flight data. Once the con- 
ductivities have been calculated they are used to determine 
ion and electron densities. 


2 . 2 Conductivity Calculations 

By making use of the slopes m Figure 2(b) and 2(c), 
both positive ion and negative conductivities can be cal- 
culated. The positive conductivities are claculated by use 
of the formula: 


cr ' 
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The formula for calculating negative conductivities is: 
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In the above formulae R and r are the guard ring and 
collector radii, respectively, 
calibration resistance. The subscripts FLT1 and FLT2 refer 
to the time in the flight that the data slopes were taken 
(see Figure 2(c)). 

During design modifications of the blunt probe the two 
radii mentioned above have undergone changes. Throughout 


most of these changes in design the R/r ratio has remained 


nearly constant, (approximately 2.0 per cm). At the present 
time there is some thought given to miniaturizing the blunt 
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probe for use with a Super LOKI Dart rocket as opposed to 
the currently used ARCAS and Super ARCAS rockets. This 
change would not affect the method of data reduction. 

However, if the ratio is changed the new values of collector 
and guard ring radii can be typed into, the computer,,, before,. „ , 
conductivity calculations are performed. 

Initially the purpose of this work was to examine the 
feasibility of using the Hybrid Computer Facility of The 
Pennsylvania State University to calculate positive ion and 
negative conductivities. The basis for this decision was 
that this phase of the data reduction had been the most time 
consuming and tedious portion of the entire data reduction 
process. However, after examining the remainder of the data 
analysis computations it was decided that, with certain 
assumptions, the computer could provide the accurate results 
needed for total data reduction. The assumptions made were 
that in the altitude range of 35 to 80Km temperature and 
pressure variations could be described by the data of 
1 September 1965 (CIRA, 1965). These assumptions were shown 
to be valid when J.P. Cipriano and J.D. Mitchell, 1972, 
compared the results of two separate data reductions of a 
probe flight. One of these data reductions used the data 
from CIRA, 1965, and the other used actual temperature and 
pressure data obtained from a meteorological datasonde, 
which was launched on the same day as the blunt probe. 

These assumptions were made primarily for convenience 
in generating a table of constant values of temperature and 
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pressure in the altitude range of interest. This table 
permits elimination of one step in the data reduction pro- 
cess and any errors interjected by these assumptions are 
considered negligible. If there is ever a need for greater 
accuracy there would be no difficulty in generating a new 
table. This could be done by using the computer and data 
acquired from a meteorological datasonde. 

It is now possible to discuss the methods used to 
calculate positive ion, electron, and negative ion 
densities. The method used for determining positive ion 
densities is fairly straightforward and will be presented 
next . 


4. 

2 o 3 Calculations of Positive Ion Densities, N 

Calculations of positive ion densities are greatly 
simplified if it is assumed that mobilities of different 
ion species are comparable in value. This assumption allows 
the ion mobilities to be represented by an effective small 
ion reduced mobility model (A. Dalgarno, 1962). The 
equation used for determining positive ion densities is: 


N 


cr 


e u 


e u 


Oh- 



CS) 


There are several constants in this equation: i) the 

positive ion mobility u + and the reduced positive ion 

2 

mobility Uq + , which was assumed to be 1.8 cm /v . s . 

(R.K. Cole and E.T. Pierce, 1965), li) the standard atmo- 
spheric temperature and pressure at sea level, Tq and Pq 
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respectively, and iii) the charge of an electron, e. Values 
for T and P are found in the table which are stored on disc 
and DEC tape. These tables were generated using model atmo- 
sphere data from CIRA, 1965. The positive ion conductivity 
is obtained from Equation (1) of the previous section. 

Positive ion densities and electron mobility are two 
major variables required to determine the electron density. 
Because of the complexity of calculating electron densities 
their analysis will be divided into two parts: 1) electron 

mobility calculations, and 2) electron density calculations. 


2.4 Electron Mobility, u g 

Electron mobility varies over the altitude range of 
interest and is found to be dependent on the ratio of 
electric field to atmospheric pressure (E/P), E.S. 
McDaniel, 1964. The equation used for calculation of 
electron mobility is: 


u 


v d 

T~ 


1 

1 


a E 


+ p 


(4) 


in which v^ is the electron drift velocity and is pro- 
portional to alpha and beta. Both alpha and beta are 
dependent upon the particular value of the E/P ratio at 
the point of interest. Table I is a tabulation of the 
values calculated for alpha and beta for given values of 
the electric field to pressure ratio, J.D. Mitchell, 1972. 
Once the values for these two variables are calculated 
the electron mobility can be calculated by using both the 
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Table I: Parameters for Electron 

Mobility Calculations 


E/P 

(v . /cm . -mm. Hg) 

a 

2 

(cm. -mm.Hg/v.s.) 

P 

(cm . /s . ) 

0.2 

1. 23xl0 6 

3. 2xl0 5 

0.4 

9.50xl0 5 

4.0xl0 S 

0.6 

8 . OOxlO 5 

4 .8xl0 5 

0.3 

6.00xl0 5 

6.0xl0 5 

1.0 

5.63xl0 5 

6.4xl0 5 

1.2 

4 . 91xl0 5 

7.2xl0 5 

1.4 

4.36xlO S 

7 . 6xlO S 

1.6 

4.36xlO S 

7.6xl0 5 

1.8 

4.36xlO S 

7.6xlO S 

2.0 

4.36xlO S 

7.6xl0 5 

4.0 

4. 14xl0 5 

9.0xl0 5 

6.0 

4 .OOxlO 5 

l.OxlO 6 

8.0 

3 . 95xlO S 

l.OxlO 6 

10.0 

3 . 82xlO S 

1 . lxlO 6 

12.0 

3.77xl0 5 

1.2xl0 6 

. 14 .0 

3 . 63xl0 5 

1.4xl0 6 

16.0 

3 . 59xl0 5 

1.5xl0 6 

18.0 

3 . 50xl0 5 

1.6xl0 6 

20.0 

3 . 50xl0 5 

1 ,6xl0 6 
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pressure at a given altitude and the electric field, which 
is determined from the blunt probe data. 

The method used to determine the electron drift 
velocity introduces uncertainties in later calculations of 
electron densities. Because of these uncertainties greater 
confidence is placed in electron densities considered 
relative to each other than in their absolute values 
(J.P. Cipriano, 1972) . 


2.4.1 Calculation of Electric Field, E 
The electric field referred to in Equation (4) is the 
field produced by the blunt probe as the probe potential 
varies from 0 to +10 volts. More precisely, it is the 
field produced at the point of measurement of the slope 
corresponding to the negative conductivity (see Figure 2(c), 
[A I /A t] „) . The electric field is directly pro- 

C r J-iiz 

portional to the voltage applied (V ) to the collector disc 

P 

at the time of the slope measurement and is determined from 
the formula 


E 


t r R 


(5) 


where R is again the radius of the guard ring. 

At this point in the data reduction the values of 
positive and negative conductivites , positive ion 
densities, and electron mobility are known. It is now 
possible to discuss the method of determining the electron 


densities . 
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2.5 Calculation of Electron Density, N g 

In the previous section it was shown that the electron 
mobility is a function of the E/P ratio. Since there are 
some uncertainties in the exact values of the electric field 
and pressure there will be a slight error introduced into 
the calculated values of electron mobility. This error will 
also be reflected in the calculations of electron densities. 
However, the data analysis techniques are such that the 
primary concern is not the absolute values of electron 
densities. The essential interest is the relative values 
obtained. Considering the consistent manner of data 
reduction and the interest in relative electron densities 
leads to the conclusion that the errors introduced can be 
neglected (J.D. Mitchell and L.C. Hale, 1972). The electron 
densities can be determined from an equation relating the 
density to negative conductivity. 

Negative conductivity is proportional to the sum of 
the negative ion and electron densities, N - and N g , 
respectively. The equation for negative conductivities is: 

= N~e u + N e u (6) 

e e 

where e is the charge of an electron, u is the negative 
ion mobility, and u g is the electron mobility previously 
discussed. Application of the principle of conservation 
of charge yields 

N~ - N + - N (7) 

e 

and substitution of Equation (7) into Equation (6) gives 
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- (N + - N ) e u + N e u (8) 

e - e e 

Because u g »u , Equation (8) can be written as 

cr~ = N + eu + N eu (9) 

- e e 

and solving for the electron density results in the 
equation 


N 


cr _ 


N + e u. 


e u 


( 10 ) 


e 

The only term which remains to be defined in this equation 
is the negative ion mobility. 

This mobility can be determined from the following 
relationship: 


u 


u 


T P 
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0- T 
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(ID 


where u Q _ is the reduced ion mobility, which is assumed to 
2 

be 2.3 cm /v.s. (R.K. Cole and E„T. Pierce, 1965). All 

the meteorological data in this equation are obtained in 

the manner discussed in Section 2.3 (Calculations of 

+ 

Positive Ion Density, N ). 

At this point the only remaining unknown in the data 

analysis procedure is the value of the negative ion 

/ 

density. The negative ion density can be determined by 
using the principle of conservation of charge, which was 
introduced in the previous section (Equation (7)). Sub- 
stitution of the calculated electron and positive ion 
densities into this equation yields the negative ion 


densities . 



CHAPTER III 


COMPUTER PROGRAM DEVELOPMENT 

3 . 1 Introduction 

Before a definite computer programming approach could 
be decided upon it was necessary to determine the type and 
amount of machine and program operator interaction. There 
were two possible ways in which the data could be analyzed: 
(1) total machine analysis, through complex mathematical 
and signal analysis procedures, or (2) standard hand data 
analysis procedures which could be modified for imple- 
mentation of computerized data reduction methods. The 
second method is the most advantageous approach because it 
allows observation and human interpretation of data 
(L.C. Hale and J.D„ Mitchell, 1972). This method of 
analysis can be implemented by making complete use of the 
hybrid computer facility. 

The equipment available at the computer facility 
consists of i) Digital Equipment Corporation’s (DEC) 

PDP-10 with 80K of core and two 5 megaword disc units, 
ii) a DEC type TU-20 magnetic tape unit, iii) a DEC type 
340 digital display scope, and iv) an Electronics 
Associates Incorporated (EAI) model 680 Analog computer. 

In order to make efficient use of the computer facility the 
data analysis was divided into three phases: (1) digiti- 

zation of data, (2) conductivity analysis, and (3) positive 
ion, negative ion, and electron density calculations. 
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3 . 2 Digitization of Data 

The data obtained from the flight of a blunt probe is 
recorded on standard quarter-inch magnetic tape. In order 
to perform the data analysis, the information stored on 
this tape must be digitized and transferred to magnetic 
tape which is compatible with DEC'S magnetic tape unit. 
This requires that there exist some type of interface 
between the two tape recording units. The analog computer 
was chosen as the interface. 

Because of the nature of the data recorded from a 
blunt probe flight, the analog computer interface had to 
perform two functions: (1) convert the raw data to an 

analog signal, and (2) convert the analog signal to a 
digital signal. The first function could be obtained by 
developing an analog data tachometer. 

3.2.1 The Analog Computer Data Tachometer 

The information stored on quarter-inch magnetic tape 
consists of variable pulse rate data. As the blunt probe 
sweeps through its range the information recorded takes 
the form of pulses. In order to convert this pulse data 
to the desired analog signal the analog computer data 
tachometer must consist of a pulse shaper, an integrator, 
and a signal filtering system. The pulse shaper chosen 
was a monostable with a two millisecond pulse width. The 
filter network consisted of a sharp cutoff low pass filter 
(A„S„ Jackson, 1960). Figure 3 is the schematic diagram 
of the analog computer data tachometer. 
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Figure 3: Analog Computer Data Tachometer Schematic 









19 


In order to understand the operation of the analog 
computer data tachometer, it is helpful to consider what 
takes place when a series of noisy data pulses are applied 
to the input. These data pulses trigger the two milli- 
second monostable, which in turn produces a series of clean, 
well-defined pulses. Following the monostable is an inte- 
grator with a potentiometer connected from the output of the 
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integrator to its OJ (OPERATE summing junction) and 
(INTEGRATE summing junction) terminals. The purpose of the 
potentiometer is to provide a discharge path for the inte- 
grator's capacitor, which enables the integrator output to 
return to zero before the next series of pulses arrive at 
the input. After the integration process the data takes 
the form of an analog signal. This new signal must now be 
applied to a filter to remove unwanted noise. Prior to the 
digitization the data must undergo an amplification and 
inversion. This is the purpose of the potentiometer and 
amplifier at the output of the low pass filter network. 

The performance of the analog computer data tachometer 
had to be evaluated before it could be used in the digiti- 
zation phase of data reduction. This evaluation consisted 
of optimizing the analog computer data tachometer and then 
comparing its output to the output of an existing data 
tachometer, which had been used in the past. 
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3.2.2 Optimization and Evaluation of the 
Analog Computer Data Tachometer 

Optimization consisted of adjusting the integrator 
decay time constant (potentiometer #4, Figure 3) and varying 
the cutoff frequency of the filter network. Observation of 
the output signal from the data tachometer revealed that the 
tachometer performed best when the integrator decay time 
constant was 50nsec. and the filter network had a cutoff 
frequency of 600 Hz. These criteria were satisfied when 
potentiometer #4 was set to 0.1, all potentiometers asso- 
ciated with the filter network were adjusted to 0.3768, and 
the integrators in the filter network were operating in the 
millisecond mode (MS) (refer to Figure 3). Upon completion of 
this optimization phase it was necessary to evaluate the 
performance of the analog computer data tachometer. 

The evaluation consisted of simultaneous observations 
of calibration and flight data processed by the analog 
computer data tachometer and the standard data tachometer 
previously used. Tests involving the calibration data 
revealed that the slopes reproduced by the analog computer 
data tachometer were approximately 15% lower than the slopes 
reproduced by the standard data tachometer. This reduction 
in slopes caused no great concern because the mathematics 
involved in conductivity calculations utilize the ratio of 
flight slope to calibration slope, Equations (1) and (2), 
Section 2.2. A 15 percent reduction in both calibration 
and flight slopes therefore would not alter the resultant 
conductivities. However, analysis of the flight data 
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revealed more than a simple 15 percent reduction in slope. 
The analysis disclosed that the analog computer data 
tachometer could not respond accurately enough to the fast 
rise times of data obtained early in the flight. It can be 
seen, Figure 4, that the data reproduced by the analog 
computer data tachometer (lower trace) was not an accurate 
representation of the actual flight data reproduced by the 
standard data tachometer (upper waveform) . Efforts were 
made to correct this problem by increasing the filter net- 
work cutoff frequency. This approach was not successful 
because noise became a limiting factor as the cutoff 
frequency was increased. The only conclusion which could 
explain the observed effects was that the capacitors 
involved in the integration and filtering process limited 
the rate at which the analog computer data tachometer 
could process data. Because of this inaccuracy the idea 
of an analog computer data tachometer was abandoned and a 
system was developed to make use of the standard data 
tachometer . 

3.2.3 Data Digitization System 

This new system of removing data from tape and convert- 
ing it to a usable form was similar to that used in the past 
(J.D. Mitchell, 1973, J.P. Cipriano, 1973); however, there 
was one difference. In the method previously used the out- 
put of the data tachometer was connected to a strip chart 
recorder. In the new system the output is connected to the 
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Figure 4: Analog Computer and Standard Data 

Tachometer Outputs 
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analog computer, which acts as the interface in the digiti- 
zation process. Digitization takes place by using the 
digital computer to control i) the analog computer's mode of 
operation, ii) the rate at which the analog to digital 
conversion of data takes place, and iii) the initialization 
and termination of the digitizing process. 

In the hybrid computer system the digital computer's 
control over the variables was obtained by writing a pro- 
gram, DIGPGM.F4 (see Appendix B) . This computer program 
provided the operator with the ability to control the data 
digitization process and allowed the data to be catalogued 
before storing it on magnetic tape. Operator control 
consists of answering questions about the setup of the 
analog computer and providing information about the type 
of data to be digitized, i.e., calibration or flight data. 
After this preliminary information has been processed the 
data sampling rate must be set. Then the data digitization 
process is started and its status is continually monitored. 
These two steps were performed by incorporating commands 
from HFOS (Hybrid Fortran Operating System, 1972) in the 
program. The digitized data is placed into an array which 
is periodically transferred to magnetic tape by MTIO 
(Magnetic Tape Input/Output, 1972). Figure 5(a) is a block 
diagram of the digitization system and Figure 5(b) shows 
the patching of the analog computer. 
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After the data has been catalogued, digitized, and 
stored on magnetic tape the task of conductivity analysis 
can commence. 

3.3 Conductivity Analysis 

To be able to appreciate the use of the hybrid computer 
facility in conductivity analysis one must consider the work 
involved in the standard hand reduction method; In this 
method the data recorded on quarter-inch magnetic tape were 
played through the standard data tachometer and a Sanborn 
Model 320 strip chart recorder. One disadvantage of this 
system is the need for a minimum of three rolls of strip 
chart recordings per flight. The first roll is a high 
speed record required to reduce the steep negative con- 
ductivity slopes in the first third of the flight. A 
second roll, at a slower speed, enables calculation of 
positive and negative conductivity slopes which occur 
during the middle portion of the flight. And the third 
roll, obtained by increasing the recording sensitivity of 
the Sanborn recorder, enables analysis of slopes late in 
the flight. Other disadvantages of this system are: 
i) the errors made when drawing slope lines on chart 
recordings are reflected in the values obtained for 
conductivities, .ii). the repetitive application of the 
formulae to tabulated slopes is not only time consuming 
but also tedious. However, this standard" method - of data 
reduction does have one major advantage: one has the 
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ability to easily extract meaningful results from data which 
has some noise present. The object of the computer program 
was to make use of this positive trait while minimizing 
errors and eliminating the disadvantages of the classical 
data reduction methods. 

To incorporate the ability of a human to discriminate 
between meaningful data and noise into the program requires 
the capability to observe the data. This requirement can be 
fulfilled by proper programming techniques and the use of 
DEC'S 340 digital display scope. Implementation of the 
display scope not only allows observation of data but it 
also eliminates the need for strip chart recordings, allows 
manipulation and analysis of data, and reduces errors in 
the scaling of data slopes. Manipulation of data includes 
the ability to i) amplify the observed signal to facilitate 
scaling of slopes which have small gradients, ii) perform 
time base expansions of slopes which have steep gradients, 
and iii) calculate the slopes of the signal. 

3.3.1 Conductivity Computer Program 

After becoming familiar with data acquisition and hand 
reduction methods the task of developing computer analysis 
methods began in May 1972. At that time the hybrid computer 
facility had 48K of core and did not have a disc. Conse- 
quently, the analysis program had to be stored on a DEC 
tape and operated from the DEC tape system. This type of 
system was slow and limited the size of the analysis 
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program. While the program was in the developmental stages 
it became apparent that its size was approaching the 
maximum usable core limitations imposed during timesharing 
hours. To avoid this problem the program was divided into 
subprograms called Chain Links (Digital Equipment 
Corporation, 1970). One subprogram was designated as the 
main program (LINK MAIN, see Appendix B) . 

The purpose of the main program is to control and keep 
track of the data reduction process. It consists primarily 
of questions which ask the operator what phase of the 
analysis he wishes to perform. All of the questions specify 
two possible responses. An incorrect response simply 
causes the program to type the question again. Upon 
recognition of an answer the main program either types out 
the next question or calls into core the chain link which 
performs the requested phase of the data reduction process. 
After the new chain link completes its operations, control 
is transferred back to the main link. The processing 
continues in this manner until the resultant conductivities 
with their corresponding altitudes are obtained. In 
July 1972, the analysis program was completed and the 
results obtained from it were ready for evaluation. How- 
ever, the evaluation of the program was delayed because the 
hybrid computer facility was shut down for the installation 
of the first 5 megaword disc pack. 

When the disc system became operable it was obvious 
that it was faster and more efficient than the DEC tape 
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system. For this reason the evaluation of the analysis 
program was delayed while the chain links were revised for 
disc operation. The modifications to the programs included 
installing the option of running the program from either 
system. This option insured that the data reduction pro- 
cess could be performed even if the disc system became 
inoperable. Before discussing the results of program 
evaluation the name and purpose of each individual chain 
link should be mentioned. 

DATAN (Appendix B) was the name given to the data 
analysis program. It consists of six chain links: 

(1) MAIN, previously discussed, (2) CALC, calibration slope 
analysis, (3) CALD, evaluation of flight slopes, (4) CONC , 
conductivity calculations, (5) MTAH, matches times-in- 
f light and altitudes, and (6) SIGP, automatically plots 
altitude vs. conductivity curves on the digital display 
scope. Links CALC and CALD require manipulation of the 
magnetic tape unit to bring the digitized data into core. 
This was achieved by using the Magnetic Tape Input/Output 
routine (MTIO) . Once the data were in core the program 
automatically scaled it before storing the data in an 
array for display purposes. The array had to be constructed 
in a manner that was consistent with specifications out- 
lined in the display routine package (DISUBS) . Several 
other arrays were constructed and displayed on the scope. 
These arrays included a grid, decision section, and display 
of time-in-flight information. 
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Figure 6 shows a typical cycle of calibration data with 
the operator-scaled slope. The slope was obtained by first 
placing the Lite-pen on the words TAKE SLOPE. This caused 
all other decision sections, except the ACCEPT region, to 
be ignored until a second ACCEPT was detected by the Lite- 
pen. The Lite-pen was then moved around until the asterisk 
was on the beginning portion of the calibration slope. After 
that, the Lite-pen was moved to the ACCEPT region. When the 
Lite-pen detected the light from this region, the X and Y 
coordinates of the asterisk were stored in core, and the 
asterisk was returned to the READY position. The next step 
was to select the second point on the calibration slope. 

Upon detection of light the asterisk was moved to the curve 
and a line was drawn between the first point chosen and the 
new location of the asterisk. The asterisk was moved until 
the line matched the slope of the calibration data. When 
this match was obtained the Lite-pen was placed over the 
ACCEPT region a second time. The new X and Y coordinates 
of the asterisk were used in conjunction with the first set 
to determine the slope of the data. This slope, along with 
others, was stored in an array until all the calibration 
data were analyzed. When the analysis of slopes was com- 
pleted, the Lite-pen was placed over the END region of the 
display scope. The program then proceeded to calculate the 
average calibration slope and store all the slopes and this 
average value on the disc or a DEC tape. Analysis of flight 
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data proceeds in much the same manner; however, there are 
more options available when scaling the data. 

Link CALD provides the ability to expand the magnitude 
or time base of the data. When the slope of the data is 
small it can be expanded by a factor of four. The flight 
data is displayed with a ten second time base. If the rise 
time of the data is less than two seconds the operator has 
the option of calling for a two second expansion. Both of 
these options increase the accuracy of data reduction and 
enable complete analysis without the need to process the 
data more than once. Figure 7 shows a typical cycle of 
flight data. The data were displayed and scaled when it 
was i) multiplied by a factor of four, Figure 8, and ii) 
expanded under a two second display, Figure 9. Slope 
analysis was performed the same way the calibration data 
were reduced; however, when the slopes were expanded they 
were multiplied by an appropriate scale factor. After the 
analysis was completed the slopes used in calculating 
positive ion conductivities and their times-in-f light were 
stored on disc or a DEC tape. The slopes used to calculate 
negative ion conductivities were stored with their 
corresponding times-in-f light and data required to calculate 
the E-field of the blunt probe. Once all the information 
was stored, the results were printed on the line printer. 

The next step in the analysis was to calculate the con- 


ductivities . 
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Figure 8: Flight Data Multiplied by Four 
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Conductivity calculations were performed by link CONC. 
This link uses the average calibration slope and the slopes 
obtained from the analysis of flight data to determine 
conductivities. These are proportional to the flight slope 
-d-i-v-i-d e d— by— t he— a-ve r.ag e— c a.l i b ration slo pe , Equations. (1) > 

and (2) in Chapter II. After the conductivities are 
calculated, the data is stored and then printed on the 
line printer. Once this phase is completed the times-in- 
f light must be matched to their corresponding altitudes. 

There were two possible methods of translating time- 
in-flight information to altitudes. One method was to 
extract this information from a radar plot. The second 
procedure was to read the data from a computer printout. 
Both these sources of information can be obtained from the 
launch facility. Use of this data in conjunction with 
link MTAH permits the transformation of time-in-flight data 
to corresponding altitudes. Upon completion of the match- 
ing process the computer program provides the option of 
automatically plotting altitude vs. conductivity data. 

The plotting process was performed by link SIGP. This 
link generated semi-log scales for these displays: (1) a 

plot of positive conductivities, Figure 10, (2) a plot of 
negative conductivities, Figure 11, and (3) a combination 
plot of both positive and negative conductivities’ 

Figure 12. Each plot remained on the display scope until 
a Lite-pen detection was received. The third Lite-pen 
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Figure 10: Computer Calculated Positive Conductivities 

for 31 Jan 1972 
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Figure 11: Computer Calculated Negative Conductivities 

for 31 Jan 1972 
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Figure 12 :l| Combined Positive and Negative Conductivity 
I Plot for 31 Jan 1972 



39 


detection terminated the program and concluded the con- 
ductivity analysis. 

3.3.2 Evaluation and Correction of 
Conductivity Program 

The changes from one large program to chain links and 
from operation off of the DEC tape system alone to either 
the disc or DEC tape system were successfully completed. 

In September 1972, the program was ready to be evaluated. 
Preliminary tests were encouraging. Except for errors in 
the times four expansion, everything appeared to be working. 
The expansion process did not display the data accurately 
when there was a slight bend in the data. An example of 
this problem can be seen by referring to Figure 13. 

Consider Figure 13(a) as an example of the data to be 
expanded. Figure 13(b) shows the expected expansion of 
the data, while Figure 13(c) depicts the observed results. 

An analysis of the problem revealed that the error was due 
to a peculiarity in the subroutine LINE of the display 
program package. The display scope is divided into 1024 
raster, units in the X and Y directions. (1 raster unit = 
.0095 inch) Whenever the AX or AY between two digitized 
points was less than one raster unit the subroutine would 
round off the delta to zero. Since AX was always one 
raster unit this error occurred whenever AY was less than 
one raster unit. -This error was cumulative but somewhat 
compensated for by AY's which were one raster unit or 
greater. The overall effect was to de-emphasize slight 
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curvatures in the data. The problem was corrected by using 
another display package subroutine, POINT. A statement was 
added before each call to LINE. This statement checked to 
see if a AY was less than one raster unit. If this 
situation occurred the program would perform a call to 
POINT before the call to LINE. The effects of this were 
to eliminate round-off error and display a point prior to 
generating a line on the display scope. All the chain 
links which were involved with analysis Of data on the 
display scope were modified to insure that no errors due 
to AY's less than one raster unit would effect the display 
of data. Upon completion of the modifications, tests 
revealed that there were no longer any difficulties in 
expanding data by a factor of four. 

The first conductivity analysis of flight data was 
performed on data obtained from the rocket launch of 
31 January 1972. A comparison of the conductivities 
obtained by computerized data reduction methods with those 
obtained by hand reduction methods revealed a factor of 
two error. This error implied that the ratio of flight 
slope to average calibration slope obtained by the computer 
reduction was less than the hand reduced ratio. Tests 
showed that this was not the case; in fact, the ratios 
obtained from the computer were larger than the hand 
reduced ratios. Further analysis revealed an error in a 
formula in the link which calculates conductivities (CONC) . 
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The formula was corrected and the data for 31 January 1972 
was again reduced. 

Resultant conductivities obtained from both the hand 
and computerized data reduction methods were plotted for 
comparison, Figure 14. The curves show that both reduction 
methods yield approximately the same results. Computerized 
results are more accurate than the hand reduction results 
because i) the display scope scaling of slopes is more 
accurate and consistent than hand scaling of strip chart 
recordings, ii) computer reduction techniques require that 
the data be scaled only once whereas the hand reduction 
technique requires analysis of at least three separate 
strip chart recordings, the latter increases the probability 
of scaling errors, and iii) the computer automatically per- 
forms all calculations, so another possible source of 
errors in the hand reduction method is eliminated.. 

The consistent nature of computerized slope scaling 
was observed when the previously mentioned ratio test was 
performed. Table II lists the values obtained for the. ratio 
of flight slope to average calibration slope and their 
corresponding conductivities. The tabulation includes 
computer and hand reduced results for two independently 
performed tests on the same flight data. One conclusion 
that can be obtained from this table is that the 
differences between conductivities of the two data 
reduction methods are primarily attributable to inaccuracies 
in the hand reduction of slopes. 



HAND REDUCED 
COMPUTER REDUCED 
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Figure 14: Smoothed Conductivity Curves for 31 Jan 1972 


Table II: Ratio Test Results 
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After determining that the computerized data reduction 
methods were feasible and accurate it was decided that 
computer calculation of densities should be attempted. 

3 . 4 Density Analysis 

The theory and equations associated with the calcula- 
tion of ion and electron densities were presented in 
Chapter II, Sections 3 through 6. Before calculations 
could begin, two interpolation programs had to be written. 
The first program ATP (Appendix B) was developed to 
interpolate altitude, temperature, and pressure data of 
1 September 1965. This information was obtained from 
CIRA '65 tables. This program included the capability of 
generating two data files: (1) ZTP.DAT, and (2) WIZTP.DAT. 

These files contained some of the information required to 
calculate densities from data obtained at White Sands 
Missile Range, N.M, , and NASA's Wallops Island, Va . , 
station, respectively. The reason for using the data from 
CIRA '65 tables was that constant data files proved to be 
accurate enough for the data reduction calculations 
(Chapter II, Section 2.1). If greater accuracy should ever 
be required ATP could be run again using meteorological 
data obtained from a datasonde launched on the same day 
as the blunt probe. The second program ALPBAT (Appendix B) 
interpolated the data in Table I, Chapter II. Its purpose 
was to determine the values of alpha and beta, in 
Equation 4, for 0.01 incremental changes in the E/P ratio. 
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Once this information was obtained the program generated 
the data file ALBA.BAT. This file, along with the files 
generated by ATP, were stored on the disc for use in the 
program which calculated densities (BENSIT, Appendix B) . 

3.4.1 Bensity Analysis Program 

Because of the dependency of electron densities on the 
positive ion densities (Equation (10) , Chapter II) BENSIT 
was divided into two parts: (1) positive ion density 

calculations, and (2) electron and negative ion density 
calculations. The calculation of positive ion densities 
is simple and fairly straightforward, Chapter II, Section 
2.3. Equation (3) was used to calculate these densities 
where the values of T (temperature) and P (pressure) were 
obtained from one of the files generated by ATP. The 
values of T and P correspond to the temperature and pressure 
at the altitude at which a positive ion conductivity 
measurement was made. Calculation of electron densities 
is considerably more complicated than positive ion density 
calculations . 

It was necessary to calculate the values of the reduced 
ion mobility (g-) and electron mobility (m- e ) before the 
electron densities could be determined. The reduced ion 
mobility was calculated by using Equation (11), Chapter II, 
Section 2.5, and the appropriate file generated by ATP. 
Calculations of the electron mobility were performed by 
implementation of Equation (4), Chapter II, Section 2.4. 
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The program first determined the E/P ratio and then used 
this ratio to obtain the appropriate values of alpha and 
beta from the file ALBA.DAT. These variables were then 
used in Equation (4) to determine - The only unknown 

remaining in Equation (10) was the value of the positive 
ion density at the point at which the negative conductivity 
was measured. This variable was obtained by interpolation 
of the positive ion densities which were previously 
calculated. Once all the variables were known the program 
proceeded to calculate the electron densities. The 
negative ion densities were then calculated using 
Equation (7). The results of positive ion, negative ion, 
and electron density calculations were stored on a storage 
device (disc or DEC tape) and then printed on the line 
printer. This printout concluded the data reduction pro- 
cess and an analysis of the validity of the densities was 
all that remained to be performed. 

3.4.2 Evaluation of Density Analysis 

The positive ion densities obtained from the computer 
reduction of the data from 31 January 1972 were in good 
agreement with those obtained by hand reduction techniques. 
Figure 15 is the plot of both the hand and the computer 
reduced positive ion densities. The differences between 
these plots can again be attributed primarily to the 
increased accuracy erf computerized data reduction methods 
(Section 3.3.2) . 
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A plot of hand and computer reduced electron densities 
(Figure 16) also revealed a discrepancy in the observed 
results. The electron densities calculated by the computer 
were found to be larger than those obtained from hand 
reduction. A portion of the difference can again be 
attributed to the greater slope scaling accuracy of the 
computer. However, there are two other factors which 
influence the differing values obtained for electron 
densities. One factor is that hand reduction processes 
ignore the term involving positive ion densities in 
Equation (10) . The second factor is that the hand reduction 
process uses a graph to determine the electron mobilities. 
This graph contains a curve which is an average plot of 
the electron mobilities obtained from data acquired from 
five rocket launches in 1971 (J.D. Mitchell, 1973); whereas, 
the computer calculates the electron mobilities for each 
set of data. 

Once it was determined that all the computer programs 
were yielding accurate results the task of data reduction 
began . 


HAND REDUCED 
COMPUTER REDUCED 
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Figure 16: Smoothed Electron Density Curves 

for 31 Jan 1972 



CHAPTER IV 


COMPUTER REDUCTION OF DATA 


4 . 1 Introduction 

The reduction of data is divided into three phases: 

(1) digitization of data, (2) conductivity analysis, and 

(3) density analysis. Appendix A is the operation manual 
which includes a detailed description of how to perform 
each phase of data reduction. It is not necessary to 
perform all three of the phases at one time. 

Because of the length of phase two it was subdivided 
into five operations: (1) calibration slope analysis, 

(2) flight slope analysis, (3) conductivity calculations, 

(4) matching of times-in-f light to altitudes, and (5) 
altitude vs. conductivity plotting. DATAN, which is the 
conductivity analysis program, was designed to provide the 
user with the option of performing one, several, or all of 
the five operations at once. If one of the first two 
options is chosen the results obtained at the time of 
interruption of the data reduction are stored for future 
use. When the user decides to return and complete phase 
two, DATAN asks a series of questions to determine what 
operations must be completed. After ascertaining which 
operation to perform next the data processing continues. 

Computer analysis was performed on data acquired from 
the flight of five blunt probes. These probes were 
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launched during the period of 6 January 1972 to 2 February 
1973. Table III lists the date, launch site, and 
statistics of each flight. 

4 . 2 Positive and Negative Conductivity Analysis 

The results of the conductivity analysis performed by 
DATAN are presented in a series of figures which consist of 
altitude vs. conductivity plots. There are three figures 
for each flight. The first figure of the series is a plot 
of positive conductivities. This is followed by a figure 
which depicts the negative conductivities. The final 
figure is a combined plot of both positive and negative 
conductivities. Each series of figures deals with the data 
for a particular probe flights i) Figures 17 to 19, 

06 January 1972, ii) Figures 10 to 12 (Chapter III) , 

31 January 1972, iii) Figures 20 to 22, 05 December 1972, 
iv) Figures 23 to 25, 16 January 1973, and v) Figures 26 
to 28, 02 February 1973. These series of figures are 
followed by Figure 29 which consists of smoothed curves 
of positive and negative conductivities for all five probe 
flights . 

The rocket launched on 16 January 1973 reached an 
altitude of 94.2 Km. This was the highest altitude ever 
attained by a Super ARCAS. It was believed that for the 
first time meaningful data could be obtained above 80 Km. 
DATAN was run a second time for this particular flight, 
in an attempt to extract the data at these higher altitudes. 
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Table 

III: Parameters 

for Blunt Probe 

Launchings 

Date of 
Launch 

Launch 

Site 

Local Time 
(Hours) 

Apogee 

(Km) 

06 Jan 72 

WI 

1337 

66 . 6 

31 Jan 72 

WI 

1307 

82.1 

05 Dec 72 

WI 

1314 

83.6 

16 Jan 73 

WI 

1304 

94.2 

02 Feb 73 

WSMR 

1315 

66.0 
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Figure 17: Computer Calculated Positive Conductivities 

for 06 Jan 1972 
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Figure 19<: Combined Positive and Negative Conductivity 

for 06 Jan 1972 
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Figure 21: Computer Calculated Negative Conductivities 

for 05 Dec 1972 
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Figure 23: Computer Calculated Positive Conductivities 

' for 16 Jan 1973 
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Figure 25: Combined Positive and Negative Conductivity 

Plot for 16 Jan 1973 
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Figure 27:1 Computer Calculated Negative Conductivities 
for 02 Feb 1973 
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Figure 29 : Smoothed Conductivity Curves 
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The resultant smoothed altitude vs. conductivity curves 
are shown in Figure 30. These curves indicate that the 
conductivities tend to decrease above 80 Km. This trend 
caused some concern so an investigation was initiated to 
determine if the results were accurate. 

The investigation revealed that data acquired above 
80 Km were not accurate. This conclusion was based on the 
following factors: i) photo-emission from the collector 

surface affects measurements when the probe face is exposed 
to the sun as the payload swings from side to size, ii) the 
rate at which the probe is initially falling can introduce 
errors, and iii) the probe was originally designed to 
measure conductivities in the altitude range of approxi- 
mately 40 to 80 Km; above 80 Km system saturation introduces 
errors. Although all three factors affect the data the 
third is probably the most significant. The reason for 
stressing the third factor is that a two second expansion 
of flight data revealed a probe sweep time less than 50 ms. 
The expansion made it possible to accurately scale these 
slopes. This accuracy made it possible to conclude that the 
observed conductivities were attributable to frequency 
response limitations of i) the electrometer used in the 
blunt probe, and ii) the data tachometer used in the 
digitization phase of data reduction. 
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re 30: Smoothed Conductivity Curves for 16 Jan 1973 
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4 . 3 Density Analysis 

The conductivities obtained in the previous section 
were used to determine the positive ion, negative ion, 
and electron densities. These densities were calculated 
by using the density analysis program (DENSIT) . Figure 31 
is a plot of the smoothed positive ion densities and 
Figure 32 shows the smoothed electron densities. A 
comparison of these two figures reveals that N » N 

i G 

between 30 and 80 Km. When this observation is combined 
with the requirement of charge neutrality (Equation (7) , 
Chapter II) the negative ion densities can be found by the 
relationship 

N_ - N + (12) 

Thus Figure 31 also represents the negative ion densities. 
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Figure 32: Smoothed Electron Density Curves 


CHAPTER V 


CONCLUSION 


5.1 Summary of Results 

A bank of computer programs was developed which 

i) digitized data, ii) performed conductivity analysis, 
and iii) calculated ion and electron densities. These pro- 
grams were used to analyze the data obtained from the blunt 
probe launched on 31 January 1972, The results of this 
analysis were compared with those obtained from hand 
reduction procedures (Chapter III) . This comparison showed 
that computerized data analysis was not only feasible but 
also desirable. It was feasible from the standpoint of 
yielding conductivities approximately equal to those 
obtained from hand reduction procedures. The desirability 
is exemplified by the increase in accuracy, which is 
attributable to i) greater consistency in scaling slopes, 

ii) the use of equations instead of approximate curves to 
obtain several variables, iii) increased accuracy in 
performing calculations through the elimination of possible 
human error, and iv) the elimination of the tedious nature 
of hand reduction methods. Another desirable feature of 
computer analysis is the significant reduction in the time 
required to perform the data reduction. The typical amount 
of time required to perform the hand reduction of data was 
three to four weeks. Computerized analysis has made it 
possible to reduce this time to four to six hours. 
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An important byproduct of the new data reduction system 
was the revelation of possible frequency response limita- 
tions of the blunt probe and the data tachometer. The blunt 
probe was found to be too sensitive above 80 Km and not 
sensitive enough below 35 Km. These sensitivity problems 
are believed to be attributable to the electrometer used in 
the blunt probe. Frequency response limitations of the data 
tachometer also influence the reproduction of data acquired 
above 80 Km. 

5 . 2 Suggestions for Further Research 

Even though the new computerized analysis system has 
considerably improved the data reduction process, there is 
the possibility for further improvement. The time required 
for reduction would be decreased by eliminating the need 
for human interaction in the processing of data; also, 
increased accuracy would probably be obtained through the 
implementation of statistical analysis of data. This type 
of analysis would allow the reduction of extremely noisy 
data tapes. 

Plotting of data could be extended to include the 
display of ion and electron densities. The plotting 
routines could then be expanded to include a subroutine 
which would determine and plot a curve which would best fit 
the calculated conductivities and densities. 

Current data acquisition equipment should be improved 
to insure accurate measurements of conductivities above 80 Km 
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and below 35 Km. To achieve this, the probe would have to 
be desensitized for measurements above 80 Km and made more 
sensitive for measurements below 35 Km. 

Presently there are three standard data tachometers 
which can be used in data reduction. The characteristics 
of the data tachometers differ; therefore, the resultant 
analog signals vary slightly. To insure consistency and 
greater accuracy the possibility of an analog computer data 
tachometer should be studied further. 
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A . 1 Introduction 

This manual was designed to familiarize the user with 
the digitization and analysis programs, DIGPGM and DATAN, 
respectively. Since its purpose does not include a detailed 
description of the hybrid computer facility it is recommend- 
ed that the user take the Hybrid Computer Short Course. 

This course is offered in the fall, winter, and spring 
terms. Its purpose is to acquaint the user with the hybrid 
computer facility and its capabilities. The course includes 
descriptions of the EAI 680 Analog computer and the DEC 
PDP-10 digital computer. Before discussing the use of the 
programs some introductory information as to how the manual 
is written is necessary. 

Words, characters, or statements which are underlined 
are typed onto the teletype (TTY) when the computer is in 
the monitor mode; e.g., PASSWORD , and NO EXECUTION 
ERRORS DETECTED . When the user types two control C's 
(CTRL C) the program is interrupted and control is trans- 
ferred back to the monitor mode. The control C appears 
on the user TTY as C. Questions asked by the program are 
preceded by several blanks so as not to confuse them with 
monitor responses. 

The manual is divided into three sections: (1) digi- 

tization, (2) conductivity analysis, and (3) density 
ana lysis . 
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A , 2 Digitization 

Before any of the programs can be operated the user 
must sign up for time on the sheets outside the hybrid 
computer office. The user must sign up the equipment 
needed under the project programmer number (PPN) . This 
number is the same for all programs (1001,422). The equip- 
ment required for digitization is i) the analog computer, 
ii) the magnetic tape unit (Mag Tape) , and iii) one TTY, 
When signing up for the analog computer, the user is 
required to supply the amount of core needed (7+5) and the 
words LOCK IN CORE, This informs the hybrid computer 
laboratory operator that on this particular day there will 
be a hybrid job which must lock into core. Upon arriving 
at the computer facility tell the operator that there is 
going to be a job locked in core. This insures the user 
that the core will be available for the digitization of 
data. The digitizing process is automatically terminated 
if the computer types a Hybrid Fortran Operating System 
(HFOS) error. If this happens, depress the control button 
(CTRL) and then two C's. This returns the user to the 
monitor mode. The digitizing process must be restarted 
from the beginning. If this error occurs twice inform the 
operator that something is interrupting the digitizing 
process. Other user programs which require line-printer 
listings or heavy disc usage can sometimes be the cause 


of HFOS errors. 
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A. 2.1 Logging Onto the System 

The logging-on procedure enables the user to address 
the computer from the monitor mode. The first step is to 
turn the TTY switch to LINE. A typical log-on procedure 
is initiated in the following manner: 

_^LOG 1001,422 

JOB 3 PSU HCL 504H37 II TTY7 

PASSWORD: LCH (the LCH will not be printed) 

2019 20-MAY-73 SUN 

• 

The user is now logged onto the system and is in the monitor 

mode. Before running any program it is advisable to assign 

the devices to be used to the job. The purpose of this is 

to insure that no one else can address a piece of equipment 

and interrupt the processing of data. Equipment can be 

assigned to a job in the following way: 

_jAS MTAO: (causes the magnetic tape unit 

to be assigned to a job) 

MTAO ASSIGNED 

_^AS HYB (assigns the analog computer 

to a job) 

HYB ASSIGNED 

Once the user has logged onto the system and assigned 
the necessary devices, the analog patch panel must be 
wired . 

A. 2. 2 Patching and Setup of the Analog Computer 

Figure A.l shows a typical digitizing system wired on 
the analog computer. The schematic diagram of this 
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Figure A.l: EAI 680 Analog Computer 
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digitizing system appears in Figure A. 2. The pushbuttons 
shown in Figure A. 2 are located to the left of the patch 
panel. The output terminals are located at the top of the 
patch panel in the logic section. Both the normal and 
inverted pushbutton outputs are available; however, the 
normal outputs are the only ones used. 

Wiring of the analog portion of the patch panel is 
slightly more complicated. The best way to become familiar 
with this wiring is to sit down at the analog computer or 
refer to the EAI 680 instruction manual which can be signed 
out from the hybrid computer laboratory office. Wiring of 
the analog computer proceeds as follows: (1) connect a 25 

volt power supply to the data tachometer, (2) wire the out- 
put of the tape recorder to the data tachometer, (3) connect 
the output of the data tachometer to the input of a summing 
amplifier with a gain of one (the output of this amplifier 
is normally connected to the Hi side of the potentiometer 
associated with the amplifer) , (4) connect the P terminal 
of the potentiometer to the input of another summing 
amplifier with a gain of ten (the low side of the potenti- 
ometer is normally connected to ground and is not available 
on the panel), (5) wire the output of the second amplifier" 
to an ADC (analog to digital converter) and to one of the 
plotter outputs, (6) ground the ADC's on either side of the 
one used, and (7) connect the plotter output on the oscillo- 
scope cart to the oscilloscope input. After these steps are 
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PUSH BUTTON 2 




Figure A. 2: Analog Computer Wiring Diagram 
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performed the user is ready to set oscilloscope controls and 
the potentiometer value. 

The oscilloscope adjustments involve the setting of 
the volts/division and time base controls, 1 v/div 
(calibrated) and 200 msec (uncalibrated, turn the calibrate 
knob completely clockwise) . Once this is completed the 
user can now set the mode controls on the keyboard of the 
analog computer. The keyboard is divided into three 
sections: (1) digital mode control (left side of the 

keyboard) , (2) device address (center) , and (3) analog mode 
control (right side of the keyboard) . On the digital mode 

g 

control section, depress the 10 and S (stop) buttons. On 
the analog mode control section, depress the N (normal) , 

SEC (seconds), PS (potentiometer set), and H (hold) buttons. 
The final step before digitization can begin is to adjust 
the potentiometer in Figure A. 2. In order to set this 
potentiometer the following steps must be performed: 

(1) make sure the PS button is depressed in the analog 
mode control section of the keyboard, (2) depress P 
(potentiometer) button in the device address section of the 
keyboard, (3) depress the numbered buttons, located in the 
device address section, which correspond to the number of 
the potentiometer used, (4) push the GO button in this 
section, (5) select the potentiometer setting by depressing 
four numbered buttons, and (6) depress the GO button once 
again. The value to which the potentiometer is adjusted to 
appears on the digital readout meter above the keyboard. 
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If a red light appears on the meter after the GO button is 
depressed the second time, the potentiometer was not set 
correctly. Should this occur depress the CL (clear) button 
and adjust the potentiometer to a value which differs from 
the fourth digit previously used. For an initial setting of 
the potentiometer follow the above sequence and use .5000 
as the value. When the potentiometer is set depress the 
PC (potentiometer coefficient) and OP (operate) buttons in 
the analog mode control section of the-keyboard . The next 
step is to advance the tape recorder to the point at which 
the probe separates from the rocket. At this point the 
launch noise diminishes greatly and the probe sweep has its 
steepest gradient. Play the tape and observe the data on 
the oscilloscope. The vertical deflection should be 
approximately four divisions. A deflection greater than 
four divisions will cause the data to be cut off when it is 
displayed on the display scope. If there is a drop out of 
data as the probe sweeps through its range increase the 
volume of the tape recorder slightly. Rewind the tape to 
the point of separation of the probe and observe the sweep 
again. If it is appreciably less than four divisions the 
value of the potentiometer must be increased. Reset the 
potentiometer to a larger value by performing the steps 
previously mentioned. Once the potentiometer is adjusted 
it should not be reset until another flight is analyzed. 

Both calibration and flight data for a given probe flight 
must be digitized with the same potentiometer setting. 
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A. 2. 3 Digitization of Data 

To digitize both calibration and flight data the 
digitizing computer program must be run twice. The digi- 
tized data is stored on the mag tape, one storage file 
being used for each type of data. 

DIGPGM is the name of the program which digitizes blunt 
probe data. Before running the program make sure the start 
and stop pushbuttons are not lit. If either of the buttons 
are lit depress the button to its right. To operate this 
program type the following words: 

_^RUN DIGPGM 

The program proceeds to ask the user questions about the 
data to be digitized and it requests information about the 
setup of the analog computer. Keep a record of the magnetic 
tape on which the data has been stored and what file number 
it has been given (see Table A. I). After the calibration 
data has been processed the program asks if data digitiza- 
tion is completed. A response of "Y" causes the program 
to be terminated. If another file is to be digitized reset 
the pushbuttons before typing the letter "N." Figure A. 3 
shows a typical series of events which occur when DIGPGM is 
being run. If the user does not desire to analyze the data 
at this time the magnetic tape should be removed from the 
drive and placed in its storage location. After this is 
done the user must log off the system by typing the 
following: 

.K/F 
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Table A.l: Mag Tape Storage Data 

for Rocket Launchings 


Mag File Statistics 

Tape Num. Date Type Day Type Data Launch Site 


DD-04 

3 

31 

Jan 

72 

L 

Ca libra tion 

Wallops 


4 

31 

Jan 

72 

L 

Flight 

Wallops 


5 

05 

Dec 

72 

L 

Ca libration 

Wallops 


6 

05 

Dec 

72 

L 

Flight 

Wallops 

DD-05 

1 

16 

Jan 

73 

Hi 

Ca libration 

Wallops 


2 

16 

Jan 

73 

Hi 

Flight 

Wallops 


3 

06 

Jan 

72 

Unknown 

Calibration 

Wallops 


4 

06 

Jan 

72 

Unknown 

Flight 

Wallops 


5 

22 

Jan 

71 

Unknown 

Calibration 

WSMR 


6 

22 

Jan 

71 

Unknown 

Flight 

WSMR 

DD-07 

1 

02 

Feb 

73 

Unknown 

Calibration 

WSMR 


2 

02 

Feb 

73 

Unknown 

Flight 

WSMR 





88 


^LOG 1001,422 

JOB 3 PSU HCL 504H37 II TTY7 
PASSWORD: 

2019 20-MAY-73 SUN 

.AS MTA0 : 

1TA0 ASSIGNED 

^AS HYB 
HYB ASSIGNED 

.RUN DIGPGM 


TYPE DATE OF LAUNCH AS MM-DD-YY 01-16-73 

TYPE ”CAL M FOR CALIBRATION DATA 
OR M FLT" FOR FLIGHT DATA CAL 

TYPE THE CHARACTERISTICS OF THE LAUNCH DAY. 

EITHER L, HI, H2 HI 

TYPE ANY COMMENTS OR NOTES SUCH AS NOISY, BALLOON, MAX 
ALTITUDE, UP TO 10 CHARACTERS IN LENGTH 94 . 5KM 

TYPE # FOR START SENSE LINE & # FOR STOP SENSE LINE. 

TYPE ##,## 00,01 

# OF ADC BEING USED- 1 

FILE STORAGE # ON MAG TAPE- 1 

TO START, SET THE START SENSE LINE HIGH & START RECORDER. 
TO STOP SET THE STOP SENSE LINE HIGH & STOP RECORDER. 

ARE YOU THROUGH DIGITIZING DATA? Y OR N Y 


CPU TIME: 14.92 ELAPSED TIME: 1:43.55 

NO EXECUTION ERRORS DETECTED ” 


EXIT 


- K/F r i 

JOB 3, USER L 100 1,4 22 J LOGGED OFF TTY7 
SAVED~ALL 40 FILES (695. DISK BLOCKS) 
RUNTIME 0 MIN, 15.87 SEC 


2021 20-MAY-73 


Figure A. 3: Example Run of DIGPGM 
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If the data is to be analyzed immediately after digi- 
tization the user must first deassign the analog computer 
and assign the display scope: 

^DEA HYB 
_^AS DIS 
DIS ASSIGNED 


A . 3 Conductivity Analysis 

Positive and negative conductivity analysis is per- 
formed by the data analysis program, DATAN. A description 
of this program was presented in Chapter III. The program 
was designed to allow the user to halt the data reduction 
process after any phase of the analysis and then re-enter 
at a later date. This option was incorporated to provide 
flexibility and continuity in the analysis of data. To 
operate the analysis program type the following: 

^RUN DATAN 

The program starts by asking if the operator desires a 
program description. If the operator is using the program 
for the first time the response to this question. should 
be "Y," This response causes the program to list all the 
questions aske.d and describes the action of the program 
upon recognition of an answer. If the operator types an 
invalid answer the program asks the question again. The 
listing was designed to enable the operator to analyze the 
data with a minimum of difficulty; therefore, this manual 
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only presents information on scaling of calibration and 
flight data. 

When the program is requested to perform the analysis 
of calibration data link CALC (see Chapter III) is placed 
into core. Before the analysis of data begins the program 
asks two questions: (1) what mag tape file is to be 

analyzed, and (2) what is the sweep voltage of the probe. 

At present the value used for sweep voltage is 22.0 volts. 
The operator is also instructed to estimate the sweep time 
of the voltage ramp. To enable the operator to estimate 
this time the display scope has one second interval marks 
displayed (see Figure 6, Chapter III). Upon the completion 
of the scaling of slopes the operator is requested to type 
in the estimated sweep time. Scaling of calibration slopes 
is performed in the following manner: (1) use the Lite-pen 

to initiate any operation, (2) observe the data and deter- 
mine if there are any scaleable slopes (if not, select the 
REJECT region) , (3) decide if the data should be shifted 
left (SL) or right (SR) , (4) place the Lite-pen over the 
TAKE SLOPE region, (5) move the Lite-pen to the beginning 
portion of the data slope, which will move the asterisk to 
the curve, (6) once the asterisk is on the curve select the 
ACCEPT region with the Lite-pen (once the first point has 
been accepted the asterisk automatically returns to the 
READY position) 3 (7) select the second point on the data 
curve (the computer automatically joins the two points with 
a line), (8) match the line to the slope of the curve, 
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(9) move the Lite-pen to the ACCEPT region, (10) when the 
asterisk is again in the READY position decide if another 
slope can be obtained by shifting the data (if there are no 
more scaleable slopes on a given display move the Lite-pen 
to.^NEXT SET; however, if there is more than one scaleable 
slope for a given display, repeat steps 3 through 10), 

(11) after all the calibration slopes have been analyzed 
make certain the line printer is available for use, and 

(12) place the Lite-pen over the END region. The END 
causes the program to store the data on a storage device. 
When the data has been stored the results are printed on 
the line printer. Scaling of flight slopes is performed 
by using link CALD (see Chapter III) . It proceeds in much 
the same manner as the analysis of calibration data; how- 
ever, there are different options. These options include 
the times four expansion and the two second display. 

When the flight slope analysis commences, the scope 
displays a ten second time base. Figure 7 in Chapter III 
shows the display of a typical cycle of flight data. 

Advance through the data by holding the Lite-pen over the 
REJECT region. Continue this process until scaleable data 
appears on the display scope, approximately 150 to 170 
seconds into the flight. From this point on the data 
analysis proceeds as follows: (1) select the number of 

slopes to be scaled, (2) place the Lite-pen over the "+" 

(for positive conductivity slope) or (for negative 

conductivity slope), (3) if necessary, choose a times four 
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expansion (option available for positive and negative con- 
ductivity slope analysis) or a two second display about a 
point (option available for negative conductivity slope 
analysis only) , (4) perform steps 5 through 9 of the 
calibration slope analysis, (5) if in step 1 a number 
greater than one was selected, perform steps 2 through 4 
again. After the chosen number of slopes have been ana- 
lyzed the program automatically displays the next ten 
seconds of data. If a given display of data contains 
unscaleable slopes simply place the Lite-pen over the REJECT 
region and the next ten seconds of data will be displayed. 
The expansion by a factor of four is initiated after 
selecting a + or - by performing the following steps : 

(1) place the Lite-pen on the vertical axis at the right of 
the display scope (this causes the asterisk to be moved to 
the line) , (2) move the Lite-pen over the X4 (if the data 
is not displayed correctly move the Lite-pen to NORMAL) , 

(3) place the Lite-pen over the axis again (if the data was 
off the bottom of the screen move the asterisk up the line; 
if it was off the top of the display move the asterisk down 
the line), (4) repeat step 2 and scale the data. To obtain 
a two second display perform the subsequent steps after 
selecting a - on the scope: (1) move the Lite-pen over the 

time below the point on the curve about which the expansion 
is desired, (2) place the Lite-pen on the 2 of 2 SEC D1S, 

(3) if the display must be altered move the Lite-pen to 
NORMAL, repeat steps 1 and 2, and (4) scale the data. 
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After the data has been completely uanalyzed.. make certain the 
line printer is ready for use and then place the Lite-pen 
over the END region. This concludes the analysis of slopes 
and the conductivities can be calculated. 

The conductivities are calculated by link CONC (see 
Chapter III) . Before the calculations are performed CONC 
asks if the user desires to change the probe radii. In most 
cases this will not be necessary. However, if the ratio 
discussed in Chapter II, Section 2 is other than 2.0 the new 
values of probe radii must be entered at this point. After 
the conductivities are printed on the line printer the 
times-in-f light can be matched to their corresponding alti- 
tudes . 

This matching process is performed by link MTAH (see 
Chapter III) . When the matching process has been completed 
the user has the option of obtaining an automatic plot of 
altitude vs. conductivity data or terminating the program. 
Calculation of ion and electron densities is the final phase 
of data analysis. 

A .4 Density Analysis 

The analysis of densities is performed by running the 
program named DENSIT (refer to Chapter III, Section 4) . To 
run this program type 

_o_RUN DENSIT 

When the program starts it requests that the operator type 
in the launch facility: WSMR (White Sands Missile Range) 
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or WI (Wallops Island) . From this point on the program is 
entirely automatic. The positive ion densities are calcu- 
lated first and the results are i) stored on a storage 
device, and ii) printed on the line printer. After this 
phase is completed the negative ion and electron densities 
are calculated, stored, and printed. At this point a 
complete reduction of data has been performed. Prior to 
logging off the system the data should be catalogued and 
stored on a DEC tape. 

A . 5 Storage of Data 

When the analysis is completed the results are stored 
either on the disc or a DEC tape under the following names: 
(1) CAL .DAT , calibration data, (2) AVCAS.DAT, average 
calibration slope, (3) VOST.DAT, probe sweep voltage 
divided by sweep time, (4) POS.DAT, positive conductivity 
slopes, (5) NEG.DAT, negative conductivity slopes, 

(6) PSI.DAT, positive conductivities, (7) NSI.DAT, negative 
conductivities, (8) PDEN.DAT, positive densities, and 
(9) NENI.DAT, electron and negative ion densities. If these 
data files are not coded and stored for future reference the 
data from a probe flight will be destroyed when the next set 
of data is analyzed. For this reason a coding scheme has 
been devised. 

Twelve DEC tapes have been set aside for the storage 
of coded data. They are located in the Hybrid Computer 
Laboratory in row B of the DEC tape storage cabinet, 
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tapes B-03 through B-15. Mount one of these tapes on a 
DEC tape drive and Write Enable the tape. Code and transfer 
the data from the storage device to the DEC tape in the 
following manner: 

^AS DTAN: 

DTAN ASSIGNED 
_^R PIP 

♦DTAN: MMDDYY . CAL/B «- DEV : CAL. DAT 
♦DTAN: MMDDYY. VST/B «- DEV: VOST.DAT 
♦DTAN: MMDDYY. AVC/B «- DEV : AVCAS .DAT 
*T>TAN: MMDDYY. POS/B «- DEV:POS.DAT 
j*DTAN: MMDDYY. NEG/B «- DEV:NEG.DAT 
♦DTAN: MMDDYY. PSI/B <- DEV.PSI.DAT 
♦DTAN: MMDDYY. NS I /B «- DEV :NSI.DAT 
♦DTAN: MMDDYY. PDN/B <- DEVlPDEN.DAT 
♦DTAN: MMDDYY. NDN/B «- DEV: NENI.DAT 
*fC 
• 

In the above example i) the N in DTAN indicates the number 
of the DEC tape drive to which the data is to be transferred, 
ii) the MMDDYY is the month, day, year code of the day on 
which the rocket was launched, and iii) DEV: is the device 
from which the data is to be obtained, either DSK: (for 
disc) or DTA# (the DEC tape drive) . After the coding has 
been completed type the following statement: 

♦DIR DTAN: 

This statement generates a listing of the information stored 



96 


on the DEC tape, 
taining the DEC 
the information 
is completed the 


Place this listing i„ the cannister con- 
tape. The listing serves as a catalog of 
stored on a particular Me tape. When this 
user can log off the system by. typing K/P. 



APPENDIX B 


COMPUTER PROGRAMS 
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•c ♦****+ ** * **** ****♦♦**♦***•♦* *♦ * »*** ******* **** ***♦*■♦•* **♦♦♦*♦ 

c 

c MAINPGM: PROGRAM TO DIGITIZE BLURT PROBE DAIA. 

C 

C *********************************************************** 
EXTERNAL ICRAC 

INTEGER DATA < 1200) , OLD, SIZE, DONE 
LOGICAL TRSL , STOP 
DATA DIV/10O0.U/ 

DIMENSION FI LI D ( 1 0) ,TOMP ( 1 00) 

C 

c INITIALIZATION & SETUP OF ANALOG COMPUTER 

c 

1 CALL HINIT 

CALL RSL 
IERCT=0 
INDsI 
C 

C REQUEST FOR INFORMATION AS TO DATE OF LAUNCH, TYPE 
C OF DATA, CHARACTERISTICS OF LAUNCH DAY, AND ANY 
C COMMENTS UP TO 10 CHARACTERS IN LENGTH. 

C 

type 300 

300 FORMAT ( * type DATE OF LAUNCH AS IW-DD-YY ',$> 

ACCEPT 310, FILIDC2), F1LIDI3) 

310 FORMAT (2A 5) 

TYPE 320 

320 FORMAT ( * TYPE "CAL" FOR CALIBRATION DATA*,/,* OR 
1 ”FLT" FOR FLIGHT DATA ',$) 

ACCEPT 310, FILID(A), FILID(5> 

IyPE 330 

330 FORMA T C * TYPE THE CHARACTERISTICS OF THE LAUNCH DAY.’,/ 

1, ’ EITHER L, HI, H2 ',S> 

ACCEPT 310, F1LID (6) , FILIQI7) 

type 340 

340 FORMAT ( ’ TYPE ANY COMMENTS OR NOTES SUCH AS NOISY, 

1 BALLOON, MAX ALTITUDE,',/,* UP TO 10 CHARACTERS IN 

2 LENGTH *,S) 

ACCEPT 310, FILIDC8), FILID19) 

C 

C ZEROING OF THE DATA ARRAY, AND QUESTIONS CONSERING 
C ANALOG CONTROL LINES AND FILE SYORAGE INFORMATION. 

C 

00 5 IB =1 , 100 
5 TDMP(IB)=0.0 

10 FORMAT (A 3) 

15 TYPE 20 

20 FORMAT ( * TYPE * FOR START SENSE LINE 4 t FOR STOP SENSE LINE 

1. ',/, ’ TYPE «« t *t \S) 
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,c *********************************************************** 

c 

C MA1NPGM: PROGRAM TO DIGITIZE BLUNT PROBE DATA. 

C 

C **************************************,******************** 

external icrac 

INTEGER DATA ( 11302) , OLD, SIZE, DONE 
LOGICAL TRSL, STOP 
DATA DI V/l 000.0/ 

DIMENSION FI LI D ( 10 ) ,TDMP < 1 00) 

C 

C INITIALIZATION A SETUP OF ANALOG COMPUTER 
C 

I CALL HINIT 

CALL RSL 
I£RCT=0 
INDrl 
C 

C REQUEST FOR INFORMATION AS TO DATE OF LAUNCH, TYPE 
C OF DATA, CHARACTERISTICS OF LAUNCH DAY, AND ANY 
C COMMENTS UP TO 10 CHARACTERS IN LENGTH. 

C 

TYPE 300 

300 FORMAT ( ' TYPE DATE OF LAUNCH AS MM-DD-YY ',$) 

ACCEPT 310, FIL1D (2 ) , F1LID(3) 

310 FORMAT (2A 5) 

TYPE 320 

320 FORMA T ( ' TYPE "CAL" FOR CALIBRATI ON DATA 1 ,/, ’ OR 
1 "FLT" FOR FLIGHT DATA *,$> 

ACCEPT 310, F1LID (4) , FILID<5> 

TYPE 330 

330 FORMA T ( * TYPE THE CHARACTERISTICS OF THE LAUNCH DAY.*,/ 

I, * EITHER L, HI, H2 ’,$) 

ACCEPT 310, FILIO (6) , FILIDI7) 

TYPE 340 

340 FORMAT ( * TYPE ANY COMMENTS OR NOTES SUCH AS NOISY, 

1 BALLOON, MAX ALTI TUDE , ’, /, ’ UP TO 10 CHARACTERS IN 

2 LENGTH \S) 

ACCEPT-310, FIL1D (8) , FILID(9> 

C 

C ZEROING OF THE DATA ARRAY, AND QUESTIONS CONSERING 
C ANALOG CONTROL LINES AND FILE SYORAGE INFORMATION. 

C 

DO 3 IB =1 , 100 
5 TDMP (1 3 ) =0 . 0 

10 FORMAT (A3 ) 

15 TYPE 20 

20 FORMAT! ’ TYPE # FOR START SENSE LINE & t FOR STOP SENSE LINE 

I . ',/, * TYPE ##,## ’,$) 
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ACCEPT 30, LSTAR, LSTOP 

IF(LS TAR .GT.7.0R .LSTOP .GT . 7 .OR .LSTAR ,LT .0 .OR .LSTOP . 

1LT.0) GO TO 15 
30 FORMAT (21 ) 

27 TYPE 26 

28 FORMAT ( * t OF AOC BEING USED: *,S) 

ACCEPT 30, IADC 

IFCIADC .LT.0.OR.IAOC .GT.7) GO TO 27 
TYPE 35 

35 FORMA T ( * FILE STORAGE # ON MAG TAPE = *,S) 

ACCEPT 30, I FI N 
C 

C INITIALIZATION OF MAG TAPE TO ACCEPT DATA 
C 

CALL MTlNIT(IFIN.l) 

FILID< 1 > = 1 FI N 
CALL SETDEC 

CALL MTVRI T (FIL1 D , 1 0 ,1 ERR ) 

IFCIERR .EQ.0) GO TO 25 
TYPE 26, 1 ERR 

26 FORMA T ( ’ ***MTIO ERROR: ’,12, ’ PROCESSING 

1 CONTINUED***',///) 

CALL MTSKIP (-1,0) 

CALL MTSKIP (1,0) 

IERCT=I ERCT+ 1 
I F (IERCT .GT .5) GO TO 74 
25 IERC T -Q 

C 

C SETTING OF ANALOG COMPUTER MODES AND PRINT OUT 
C OF STARTING AND STOPPING INSTRUCTIONS. 

C 

CALL STCO ( 'NSEC') 

CALL S AMO ( ’OP') 

CALL SLMO( 'RUN ') 

TYPE 40 

40 FORMAT ( ’ TO START, SET THE START 

1 SENSE LINE HIGH 4 START RECORDER.’,/,' 

2 TO STOP SET THE STOP SENSE LINE HIGH 4 STOP 

3 RECORDER. ’,/) 

C 

C SET UP DATA ACQUISITION TIMER, 4 CHECK FOR ERROR 
C IN SETUP. THE CALL TO "I NTC NT" TELL THE COMPUTER THE 

C NUMBER OF INTERRUPTS TO PERFORM, A ZERO IMPLIES 

C CONTI NI OUS PROCESSING. 

C 

CALL DTM (2000, I£RR,ICRAC, IADC , 1000, DATA) 

IFUERR.GT.0) TYPE 45, I ERR 

45 FORMATC' - DIM ERR = * ,1 27’“*>>* PGM TERMI NATED ****’,////) 

IFOERR .GT.0) GO TO 100 
CALL INTCNT(0> 
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C 

C TEST TO SEE IF SAMPLING SHOULD BEGIN. 

C 

50 I F (TRSL (LS TAR > ) GO TO 60 

GO TO 50 
C 

C START INTERRUPT 
C 

60 CALL STINT 

sizErioaa 
FREErSIZE 
OLD:0 
C 

C TEST TO SEE IF SAMPLING IS OVER 
C 

65 STOP =TRSL (LS TOP ) 

CALL DTMSTA (STOP, IS TAT, FREE .DONE) 

C 

C TEST TO SEE HO'a! MUCH DATA HAS BEEN PROCESSED. 
C 

I F (DONE -OLD ) 66,67,68 

66 IF (OLD -SIZE) 69,70,70 
C 

C PROCESS DATA AND STORE IT IN ARRAY "TDMP" , 

C THEN DUMP ARRAY "TDMP" ON TO MAG TAPE. 

C 

69 DO 71 1 =0LD+1,S1ZE 
TDMP(IND):FLOAT(DATA(I ) )/DIV 
IFdND.EQ .100) CALL MTWRI T (TDMP , 1 00 ,1 ERR ) 
IF(IND.EQ .100.AND.IERR .NE.0) ISTAT=-1 
IFdNO.LT. 100) GO TO 213 

C 

C ZERO ARRAY "TDMP". 

C 

DO 200 18:1,100 
200 TDMP (IB):3.0 

I ND : 1 
GO TO 71 

210 I ND :I ND+ 1 

71 lFdSTAT) GO TO 67 

70 OLD=0 
C 

C PROCESS DATA AND STORE IT IN ARRAY "TDMP", 

C THEN DUMP ARRAY "TDMP" ON TO MAG TAPE. 

C 

68 DO 72 I =OLD+ 1 .DONE 

TDMP(IND)=FLOAT(DATA (I))/DI V 
IFUND.EQ . 100) CALL MT VRI T (TDMP , 100, IERR ) 
IF(IND.E3. 100. AND .IERR .NE .0) ISTATr-1 
IFUND.LT. 100) GO TO 220 
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C 

C ZERO ARRAY "TOMP". 

C 

DO 230 13:1,100 
230 TDMP (IB )=0.0 

IND:1 
GO TO 72 
220 UD:INQ+l 

72 IFU5TAT) GO TO 67 
FREE :DONE 
OLDiDONE 

67 IF(.NOT.ISTAT) GO TO 65 

C 

C CHECK THE STATIUS OF THE DATA ACQUISITION TIMER 
C AND REPORT ANY ERRORS. 

C 

CALL DTMSTA (.TRUE., I STAT, FREE .DONE) 

74 1F<IERR.GT.0) TYPE 73, IERR 

73 FORMAT ( ' **** ERROR IN PROCESSING: MTIO OUTPUT 

1 ERR: ',12, * ****•,/,• *,4X, ’CONTROL TRANSFERRED 

2 TO START OF PGM. *,///) 

I F (IERR .GT .0 ) GO TO l 

C 

C ENDS OUTPUT ON MAG TAPE. 

C 

CALL MTWRIT(T3MP, 100, IERR) 

I F (IERR .NE . 0) TYPE 79, IERR 

79 FORMAT! * ***MT10 ERROR ',11,' ON FINAL OUTPUT***' ,/, * 
1THE VERY LAST DISPLAY MAY BE IN ERROR, NORMAL 

2 TERMINATION IN EFFECT*,//////) 

C 

C CLOSE THE MAG TAPE FILE AND TYPE THE FINAL 
C QUESTION. 

C 

CALL MTCLOS 

80 TYPE 90 

90 FORMAT ( ' ARE YOU THROUGH DIGITIZING DATA? Y OR N 

1 \S) 

ACCEPT 10, ANS 
I F(ANS .EG. *N *) GO TO 1 
I F(AN3 .NE. *Y ') GO TO 80 
STOP 
END 


100 
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C * 4 ***************************************************** 

c 

C BLKDAs BLOCK DATA 

C 

C *************i*t******* *********************** ************ 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


BLOCK DATA 


COMMON STATEMENT FOR DISPLAY SCOPE. 

COMMON LP,I5H0W,XMAX,XMIN,YMAX,YMIN,INTENS , I SCALE 

COMMON AREA FOR DATA AND VARIABLES TRANSFERRED 
FROM ONE CHAIN LINK TO ANOTHER. 

DISPLAY ARRAYS: 

ADVAL<3003): DATA ARRAY 

GRI D ( 1 000 ) : GRID SETUP FOR DISPLAY 

DEC (100)! USED FOR "CAL" AND "DATA” OPTIONS IN 

LINK "MAIN” AND "2 SEC DIS" IN LINK "CALD. 
XA (5) & YA ( 5 ) ! SEPERATION LINES BETWEEN DISPLAY 
LITERALS. 

STAR (20) 4 SEC ( 1 0 ) MOBILE STAR 

ILINECIC0)! NEEDED TO DRAW A LINE BETWEEN POINTS ON 
CURVE. 

TIME(250>! TEN SECOND TIME SCALE. 


OTHER MOST SI 
A VC AS L: 

SLOPE (150,2 


TI F ( 1 50 , 3 ) ! 

SI GPOS < 1 50) 

SIGNEG( 150) 

TI .MEI F( l 50, 
FI LID ( 1 0 ) : 

FVAL( 1005)1 

C VAL ( 1 600 ) : 

ALT (150,2)! 
IDSKi 


SLOPE (LINKS CALC A 


GNIFICANT VARIABLES! 

AVERAGE CALIBRATION 
CONC). 

) iC a LI BRA T I ON SLOPE AR**Y IN LINK "CALC", 
POSITIVE AND NEGATIVE SLOPE ARRAY IN 
LINK "CALD". 

TIME IN FLIGHT ARRAY USED IN LINKS 
"CALD" AND "MTAH”. 

: POSITIVE CONDUCTIVITY ARRAY USED IN 
LINKS "CONC”, "MTAH", AND "SIGP". 
l NEGATIVE CONDUCTIVITY ARRAY USED IN 
LINKS "CONC”, "MTAH", ANS "SIGP". 

2)! TIME IN FLIGHT ARRAY USED IN LINK "CONC", 
ARRAY WHICH CONTAINS INFORMATION 
PERTAINING TO THE FLIGHT. 

ARRAY CONTAINING DIGITIZED FLIGHT DATA 
READ IN FROM MAG TAPE, USED IN LINK "CALD". 
ARRAY CONTAINING DIGITIZED CALIBRATION 
DATA READ IN FROM MAG TAPE, USED IN 
LINK "CALC”. 

ARRAY FOR STORAGE OF ALTITUDES, USED IN 
LINKS "MTAH" AND "SIGP". 

I DSK - 1 INDICATES DISK OPERATION, OTHER 
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VALUES OF IOSK INDICATE DEC TAPE 
OPERATION. 

COMMON /A/ ADVAL<3000), GRIDU000), DEC (100) /3/ XA(5), YA < A) 

1 /C/ AVCASL /D / STAR (20 ) , SEC(U) /£/ SLOPE<150,2> 

2 /F/ X ( A ) , Y (A) /G/ ISKIP A\ / I1F<150,3) 

3 /I/ TIME (250 ) /J/ ASK ( A ) , RC6), B(2)/K/ SI GPOS ( 150) , 

A SI GNEG (150) /L/ TIME1 F< 1 50,2) /M/ IND, I FLAG, 

5 1E0F , KIN, INDF , NFLG /N / IPSI, INSI 

6 /P/ IL1 N£ ( 100) , RB ( 5) /Q / ALTU50,2) 

7 /R / ASK 1 (2) /3/ ‘FVALUaOS) /T/ CVAL( 1600) 

8 /U/ I, J /DA /IOSK , NUM /DB/ NUM2 , LB /DC/ FILiDC 10) , 

9 TDMPC100) 


SETTJNG OF CONSTANTS. 

DATA XA/0. 0,380.0, 630.0,880.0, 1023.0/, 

1 YA/0.0, 180.0,630.0, 1023.0/, XMAX/1023.0/, XMIN/0.0/, 

2 YMAX/1023.0/, YMIN/0.0/, B/695. 0,958.0/, 

3 R8/380. 0.A33.0, 500. 0,5A8. 0,600.0/, LB/10/, NUM/100/, 

4 IDSK/1/, NUM2/100/ 

END 
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C ******* ********** ********************************************** 

C DATANs asks if dsk or dec tape is being used, 

C AND IF HELP IS NEEDED TO RUN PGM. 

c *************************************************************** 

COMMON /DA/ I DSK 


C 

C QUESTION TYPED TO DETERMINE IF A DESCRIPTION 

C OF THE PROGRAM IS DESIRED. IF NOT CONTROL IS 

C TRANSFERRED TO STATEMENT 60. 

C 

30 TYPE 40 

40 FORMA T ( ' DO YOU DESIRE A PROGRAM DESCRIPTION? 

1 Y OR N *,I) 

ACCEPT 50, ANS 
50 FORMAT (A2) 

I F(ANS .EQ . 'N * ) GO TO 60 
IF (A NS .N£. 'Y *) GO TO 30 
C 

C TYPE OUT EXPLANATION OF ENTIRE PROGRAM 

C 

PRINT 100 

100 FORMAT ( * 1 * , * THIS PRINT OUT IS AN EXPLANATION OF 

1 ALL QUESTIONS WHICH THE PROGRAM WILL ASK. IT IS AIMED AT 

2 GIVING THE USER',/,' ’, *A STEP BY STEP DESCRIPTI 
3 ON OF HOW THE PROGRAM OPERATES .',////, * ********** 
4********* YOU ARE NOW IN IHt STARTING PROGRAM, IIS NAME IS 
5 DATAN +*******************+**’,//, 

6* DATAN IS TO BE USED FOR DATA ANALYSIS ONLY AFETER 
7 DATA HAS BEEN DIGITIZED AND STORED',/,' ON MAG TAPE. 

8F0R A DESCRIPTION OF THE DIGITIZING PROGRAM (DIGPGM) SEE 
9 THE APPROPRIATE',/,' SIDE REPORT.’,////,' \10X,’EXPLA 
1NATI0N OF THE QUESTIONS TO BE ASKED BY THE PROGRAM*,//, 

2’ THE NEXT QUESTION TO BE ASKED PERTAINS TO THE SYSTEM 

3 YOU ARE OPERATING UNDER.',/,' ',10X,’IF THE DISK 

4 SYSTEM IS INOPERABLE YOU ARE OEVIOUSLY RUNNING OFF THE 

5 DEC TAPE SYSTEM, SO TYPE IN’,/,’ ’,10X,'THE NUMBER OF 
6THE DEC TAPE DRIVE THE MAIN TAPE IS ON. THIS NUMBER 
7MUST BE IN THE RANGE OF I THRU 8.’,/,’ ’,10X,*IF YOU 
8ARE PERATING ON THE DISK SYSTEM TYPE THE NUMBER 

9 10. ’,///> 

PRINT 200 

200 FORMAT ( ' THE PROGRAM HAS BEEN CONSTRUCTED TO CHECK FOR 
1ANY WRONG RESPONSES. UPON DETECTION OF AN’,/,’ INCORR 
2 EOT RESPONSE THE PROGRAM WILL ASK THE SAME QUESTION OVER 
3.’,///, ' AFTER ANSWERING THE INITIAL QUESTION THE PROG 
4RAM BEGINS BY ASKING A SERIES OF QUESTIONS PERTAINING 
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5 TO THE STATUS',/,* OF THE DATA REDUCTION. FROM 

6 THIS POINT ALL QUESTIONS WILL BE NUMBERED IN THIS 

7 DESCRIPTION. ’,//) 

PRINT 250 

250 FORMA T ( '- BECAUSE OF THE TREMENDOUS SIZE OF THE PROGRAM 
l IT HAS BEEN BROKEN DOWN INTO SEVEN SEPARATE *,/, * PRO 
2GRAMS CALLED CHAIN LINKS. AT THE CORRECT POINTS IN THE 

3 DATA REDUCTION CONTROL IS AUTOMATICALLY TRANSFERRED 

4 TO THE',/,' APPROPRIATE CnAIN LINK. FOUR OF THE CHAIN 

5 LINKS PRINT OUT THEIR FUNCTION AND YOU HAVE THE OPTION 

6 OF CONTI NUE1 NG THE*,/, 1 PROCESSING OR EXITING. BEFORE LISTING 

7 THE QUESTIONS HERE IS A BRIEF DESCRIPTION OF THE’,/,’ 

8 NAMES AND PURPOSE OF SOME OF THE CHAIN LINKS :',/,*-', l 5X 
9, 'MAIN: QUESTIONING PROGRAM WHICH DETERMINES STATIUS 

1 OF DATA REDUCTION AND',/,’ *,20X. 'LINKS TO THE 

2 APPROPRIATE CHAIN LINK TO BE USED.*,/, * ALL OTHER LINKS 

3 CHAIN BACK TO "MAIN" AFTER PERFORMING THEIR OPERATIONS.', 

4/,' *, 1 5X , ‘CALC : CALIBRATION SLOPE CALCULATIONS 1 ,/, * ' , 1 5X , 
5'CALD: FLIGHT SLOPE CALCULATI ONS 1 , / , ' \15X,'C0NC; COND 
6UCTI VI TY CALCULATIONS 15X, 'MTAH : MATCH TIMES AND 

7 ALTITUDES*,/,’ *,15X,*SIGPt ALTITUDE VS. CONDUCTIVITY 

8 PLOTS *, //) 

EXPLANATION OF QUERTIONS ASKED IN LINK MAIN. 

PRINT 300 

300 FORMAT ( * 1 QUESTION (Is*,/,* *,I5X,*IF YOU HAVE CONDUCTIV 

II TIES CALCULATED AND STORED, DO YOU DESIRE TO’,/,* ' , 1 5X , 
2'MATCH TIME IN FLIGHT WITH ALTI TUDE 7 ’ , /, * *,22X,*P0SS 
3IBLE RESPONSES: "Y"...YES OR ”N" . . .NO. * , /, * * , 1 5X 
4, *A RESPONSE OF "Y " CAUSES CONTROL TC BE TRANSFERRED 
5 TO THE CHAIN LINK "MTAH", IF AND ONLY IF*,/,’ *,15X,*0PER 
6ATI0N IS OFF THE DISK SYSTEM. IF OPERATION IS OFF DEC TAPE 

7 SYSTEM THE PROGRAM REQUESTS’,/,’ ’,I5X,’THE NUMBER OF 

8 THE DEC TAPE DRIVE WHICH CONTAINS PREVIOUSLY REDUCED DATA 
9. A RESPONSE OF "N" CAUSES’,/.’ ' , 1 5X , ’QUESTI ON #2 TO 

1 BE TYPED.’,//, ’ QUESTION #2:’,/,’ \15X,’IF YOU ALREADY 

2 HAVE CONDUCTIVITIES MATCHED WITH CORRESPONDING ALTITUDES 
3,',/,’ * , 1 5X , 'DO YOU DESIRE TO PLOT ALTITUDE VS. CONDUCTI 
4VITY CURVES?’,/, ’ ' , 20X , ’POSSI BLE RESPONSES: "Y"...YES OR 

5 "N"...NO. ’,/, ' ' , t 5X , 'A RESPONSE OF "Y" CAUSES CONTROL 

6 TO BE TRANSFERRED TO THE CHAIN LINK "S1GP", IF AND ONLY', 

7/,' ' , 1 5X , 'IF OPERATION IS OFF THE DISK SYSTEM. IF OPERA 
BTION IS OFF THE DEC TAPE SYSTEM THE PROGRAM REQUESTS') 

PRINT 350 

350 FORMAT!' ',15X,'THE NUMBER OF THE DEC TAPE DRIVE WHICH CONT 
IAINS PREVIOUSLY REDUCED DATATA RESPONSE OF ’’N"',/,* * , 1 5X 
2, 'CAUSES QUESTION #3 TO BE TYPED.',//) 

PRINT 400 

400 FORMAT!* QUESTION *3:',/,' \I5X,'IF YOU HAVE SLOPES 
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1 CALCULATED AND STORED DO YOU DESIRE TO CALCULATE THE 

2 CONDUCTIVITIES * ' , 1 5X , 'OR DO YOU DESIRE TO REDUCE 

3 NEW DATA?',/,* 1 , 20x , ‘BOSS I BLE RESPONSES: "0"...0LD OR "N" 
A. ..NEW.*,/, ’ ',15X,’A RESPONSE OF "0" CAUSEC CONTROL 

5 TO BE TRANSFERRED TO THE CHAIN LINK "CONC".*,/, 

6* \I5X,*TH1S TRANSFERR OF CONTROL TAKES PLACE ON BOTH DEC 
7 TAPE AND DISK SYSTEMS OF OPERATION. A hESPONSE OF "N" 
S',/,’ ', 1 5X, ’RESULTS IN TYPE OUT OF QUESTION »*.',//, 

9 ’ QUESTION it*:',/,' \15X,’D0 YOU KNOW IF THE FILE YOU 

1 DESIRE TO ANALYZE CONTAINS CALIBRATION OR FLIGHT 

2 DATA?',/, ‘ *,20tf, ‘POSSIBLE RESPONSES: "Y"...YES OR "N“ 

3 ... NO , ' , * * , 1 5X , 'A RESPONSE OF "Y" CAUSES QUESTION #6 
A TO BE TYPED. AND A RESPONSE OF "N" CAUSES QUESTION #5 

5 TO BE TYPED. ',//) 

PRINT 500 

500 FORMAT ( ' QUESTION #5:*,/, * ',15X,'TYP£ THE NUMBER OF THE 
I MAG TAPE FILE OF THE DATA WHICH IS TO BE OBSERVED.*,/ 

2,* *,20X, 'POSSIBLE RESPONSES: ANY MAGNETIC TAPE FILE 

3 NUMBER. *,/,* * , 1 5X , 'THE PROGRAM WILL DISPLAY THE DATA 
A CONTAINED IN THE APPROPRIATE MAG TAPE FILE AND AWAIT 

5 A LITE PEN',/, ' ', 15X , 'DETEC TI ON IN THE "CAL” 

6 (CALIBRATION DATA) OR "DATA" (FLIGHT DATA) REGIONS. 

7 CONTROL IS THEN TRANSFERRED',/,' ’,15X.'T0 THE APP 
8R0PRI ATE ChAlN LINK.',//,' QUESTION #6:',/,' ' , 1 5X , 

9'IF THE FILE YOU WISH TO ANALYZE IS A CALIBRATION FILE 
1 TYPE CAL; IF IT IS A FLIGHT DATA FILE TYPE FLT . ' 

2,/,' \20X, 'POSSIBLE RESPONSES; CAL OR FLT.',/,' ' , 1 5X , 

3 'A RESPONSE OF CAL CAUSES CONTROL TO BE TRANSFERRED TO 

4 CHAIN LINK "CALC", AND A RESPONSE OF FLT CAUSES’,/, 

6* ',15X, 'CONTROL TO BE TRANSFERRED TO CHAIN LINK "CALD".' 
7,//) 

PRINT 600 

600 FORMAT! ’- AFTER COMPLETION OF ANALYSIS OF CALIBRATION 

1 OR FLIGHT SLOPES THE FOLLOWING QUESTIONS',/,' ARE 

2 ASKED.',/, '1 ', ' QUESTION #7: ',/, ' ' , I 5X , 'AT THIS 

3 POINT DO YOU WISH TO COMPUTE CONDUCTIVITIES?’,/,' ', 

420X, 'POSSIBLE RESPONSES: ”Y"...YES OR "N" . . . NO. ' , /, 

5* ', 15X, 'IF THE RESPONSE IS "Y” CONTROL IS TRANSFERRED 

6 TO CHAIN LINK "CONC". IF THE RESPONSE IS’,/,' ' , 1 5X , ' 

7"N" QUESTION #10 IS ASKED.',/,’- AFTER CONDUCTIV 

81 TI ES HAVE BEEN CALCULATED IN "CONC" CONTROL IS TRANSF 
9ERRED BACK TO CHAIN LINK "MAIN", AND',/,' QUESTION 

1 #8 IS TYPED.',/,'- QUESTION #8:',/,' \15X,’SAME AS 

2 QUESTION #1 BUT WITH A DIFFERENT RESPONSE TO AN',/,' ', 
315X, 'ANSWER OF "N".’,/,' ’,20X,’A RESPONSE OF "N" CAUSES 
A QUESTION #10 TO BE TYPtD .',/,’- AFTER TIMES HAVE 

5 BEEN MATCHED TO ALTITUDES, IN CHAIN LINK "MTAH" , CONT 
6R0L IS TRANSFERRED BACK TO CHAIN',/, ' LINK "MAIN" AND 

7 QUESTION #9 IS ASKED.’,//) 

PRINT 700 
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700 FORMAT < QUESTION # 9:',/,* \15X,*D0 YOU DESIRE TO 

l PLOT ALTITUDE US. CONDUCTIVITY CURVES?’,/, * \23X, 

2 ’POSSIBLE RESPONSES: "Y"...YES OR " N" . . . NO. ’ , / , ’ ’, 

3 1 5X « ’A RESPONSE OF "Y " CAUSES TRANSFERR OF CONTROL TO 

4 CHAIN LINK "SIGP", AND A RESPONSE OF "N"’,/,’ ’ , 1 5X , * 
5CAUSES QUESTION #10 TO BE ASKED. ’,/,’ QUESTION #10:’, 

6 /,’ ' , 1 5X , ’ARE YOU FINISHED WITH THE PROGRAM?’,/,’ ’, 

720X, 'POSSIBLE RESPONSES: "Y"...YES OR "N” . . .NO. ’,/, ’ ’ 

8 , 1 5X , ’A RESPONSE OF “Y " CAUSES TERMINATION OF THE 

9 PROGRAM BUT A RESPONSE OF "N" TRANSFERRS CONTROL’,/, 

1* ’, 1 5X , ’TO THE START OF THE PROGRAM AND THE QUESTION 
21 NG SEQUENCE AGAIN BEGINS WITH QUESTION #1’,//) 

PRINT 800 

800 FORMAT ( ’ THE ONLY OTHER INFORMATION WHICH THE OPERATOR 

I MUST PROVIDE IS THE NUMBER OF THE MAG TAPE FILE WHICH’,/, 
2* CONTAINS THE DATA TO EE ANALYZED 5 OR, THE NUMBER 

3 OF A DEC TAPE DRIVE ON WHICH DATA IS TO BE STORED OR 

4 READ IF’,/,’ OPERATION IS OFF THE DEC TAPE SYSTEM. 

5',/, *, 5X , ' ********************** 

S************* END DISCREPTION SECTION ************* 
T**********m*********** ’ , /, 

8 1 5X , ’YOU ARE NOW READY TO BEGIN DATA 
9 PROCESSING. :::::::::: GOOD LUCK Ml I 
1 ’I’,/, ’!’,/, ’1 ’,/, ’1 ’) 

C 

C REQUEST FOR INFORMATION AS TO WHETHER THE 
C PROGRAM WILL BE OPERATING OFF THE DISK OR DEC 

C TAPE. IF OPERATION IS OFF DEC TAPE ADD 8 

C TO THE DEC TAPE NUMBER AND THEN CHAIN TO LINK 
C "MAIN”. OTHERWISE "IDSK" WILL EQUAL ONE FOR 
C DISK OPERATION AMD THE PROGRAM WILL AUTOMATICALLY 
C CHAIN TO "MAIN". 

C 

60 IBAD=0 

TYPE 13 

10 FORMAT ( ’ IF YOU ARE RUNNING OFF OF DEC TAPE TYPE 

1 THE # Or THE DTA ' , /, ' OTHERWISE TYPE 10 ’, k) 

ACCEPT 20, IDOT 
20 FORMAT (I) 

IFUDOT.LT.l .0R.1D0T.GT.8) IBAD:1 
IFdBAD.EQ. 1 .AND. IDOT. NE. 10) GO TO 60 
IF (IDOT. EQ. 10) GO TO 70 
IDSK =ID0T+8 

CALL CHAIN(0, IDSK , ’MAIN ') 

END 


70 
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c *************************************************************** 

C MAIN: USE OF DISPLAY SCOPE TO CALCULATE 

C CONDUCTIVITIES 

C ********************************************* ****************** 

COMMON LP, I SHOW, XMAX, XMIN, YMAX , YMIN, INTENS, ISCALE 
COMMON /A/ A D VAL (3 U3U ) , GR I D C 1 000 ) , DEC (102) /B/ XA(5), YA(4) 
l /C/ AVCASL /G/ 1SKIP /N/ IPSI, 1 NSI /DA/ 1 DSK , NUM 
2/DC/ FILIDC12) ,TD'MP (100) /S / FVAL( 1 005 ) 


COMPUTED GO TO REQUIRED TO DETERMINE WHERE CONTROL 
SHOULD BE TRANSFERRED AFTER OTHER LINKS CHAIN 
BACK TO LINK "MAIN". 

GO TO (165, 150, 156, 119, 124, 112, 115), NUM 
I PLT :0 

QUESTIONING SECTION TO DETERMINE ANY PRIOR STATIUS 
OF THE REDUCTION. 

TYPE 7 

7 FORMA T ( ’ IF YOU HAVE CONDUCTIVITIES CALCULATED 

1 & STORED,’/, * DO YOU DESIRE TO MATCH 

2 TIME IN FLIGHT WITH ALTITUDE? Y OR N ’,S) 

ACCEPT 10, MTA 

IF(MTA.EQ.’Y') GO TO 157 
I F (MTA . NE . ‘N ’ ) GO TO 6 
B TYPE 5 

9 FORMAT ( ’ IF YOU ALREADY HAVE CONDUCTIVITIES MATCHED 

1 WITH CORRESPONDING’,/, ’ ALTITUDES, DO YOU DESIRE 

2 TO PLOT ALTITUDE VS. CONDUCTIVITY’,/,’ CURVES 

3 Y OR N? \S) 

ACCEPT 10, PLT 

I F (PLT .EQ . *Y ’) IPLTzl 
IF(IPLT.EQ.l) GO TO 157 
IF(PLT.NE.’N’) GO TO 8 
2 TYPE 5 

5 FORMAT ( ’ IF YOU HAVE SLOPES CALCULATED AND 

1 STORED,’,/, ’ DO YOU DESIRE TO 

2 CALCULATE THE CONDUCTIVITIES',/,’ OR DO YOU 

3 DESIRE TO REDUCE NEW DATA??’,/,’ 

4TYPE "0" FOR OLD..."N" FOR NEW. ’,$) 

ACCEPT 10, ITY 

IFUTY.EQ .’O’) GO TO 155 
IFdTY.NE. ’N ’) GO TO 2 

10 FORMAT (A3) 

30 FORMAT (31) 
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1 LF = 1 

ISCALE=0 
YI =200 • 0 
C 

C INITIALIZATION OF GRID ARRAY; THIS SETS UP 
C THE GRID FOR THE DISPLAY. 

C 

CALL I NTAB (GRI D , 1000 ) 

I NTENS = 5 
ISHOWsl 
C 

C GENERATION OF HORIZONTAL LINES OF GRID 
C 

CALL POINKGRID, 0.0, 160.0) 

CALL VECONT (GRID, 0.0, 160.0, 1023.0, 160.0) 

ADV=0,0 
DO 140 1=1,9 

CALL PCINT(GRID,0.0,ADV+200.0) 

CALL VECONT (GRID,0.0,ADV+200.0, 1023.0, AD 1/4-200.0) 
140 ADV=ADV+100.0 

C 

C GENERATION OF VERTICLE LINES FOR GRID 
C 

ADV=0.0 
DO 180 1=1,1 l 

CALL POINT(GRID,ADV, 160.0) 

CALL VECONT (GRID, ADV.t 60. 0,ADV, 1023.0) 

180 ADV =ADV+ 100.0 

C 

C GENERATION OF EASE LINE IN GRID 
C 

I NTENS =7 
ISHOU= 1 

CALL P01 NT (GRI D , 0.0, 0.0) 

CALL VECONKGRID, 0.0, 0.0, 1023.0, 0.0) 

CALL POI NT (GRID , 0.0, 40.0) 

CALL VECONKGRID, 0.0, 40.0, 1023.0, 40.0) 

CALL POI NT (GRID, 0.0, 99. 9) 

CALL VECONT (GRID, 0.0, 99. 9, 1023.0,99.9) 

C 

C TICK MARK TO SHOW 1 SEC INTERVALS. 

C 

IIMK =100.0 
DO 185 1 =1 ,10 

CALL POINT(GRID,TIMK , 150.0) 

CALL LIN£ (GRID .TIMK , l 50.0, TIMK , 1 65.0) 

185 TIMK =TIMK+ 1 00 .0 
C 

C GENERATION OF SEPERATI ON LINES BETWEEN 
C LITERALS ON DISPLAY SCOPE 
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c 

DO 190 J : 1 , 5 

CALL POINT (GRID, XA(J>, 0.0) 

190 CALL LINE (GRID, XA(J), 0.0, XA(J), 40.0) 

00 195 J:1 ,4 

CALL POINT! OR ID, YACJ), 40.0) 

195 CALL LINE (GR I D , YA (J ) , 40.0, YACJ), 99.9) 

C 

C GENERATION OF LITERALS FOR DISPLAY SCOPE 
C 

CALL I NTA3 (DEC , 100) 

CALL POINT (DEC, 180.0,0.0) 

CALL LINE (DEC, 180.0,0.3, 180.0,40.0) 

CALL CHRGENOEC, 'CAL.', 20.3, 5.0, 4, 0) 

CALL CHRGENCDEC, ’DATA', 230.0, 5.0, 4, 0) 

CALL CHRG£N(GRID, 'ACCEPT*, 430.0, 5.0, 4, 0) 

CALL CHRG£N(GRID, 'REJECT*, 680.0, 5.0, 4, 0) 

CALL CHRGENCGR1D, 'END * , 9 1 0 . 0 , 5 .0, 4 , 0) 

CALL CHRGEN (GRID , 'READY ',2.0,55.0,4,0) 

C 

C GENERATION OF VERTICAL 4 HORIZONTAL 
C ZERO LINES 
C 

CALL POI NT (GRI D , 3.0,160.0) 

CALL VECONT (GRID, 0.0, 160.0,0.0, 1023.0) 

CALL POI NT (GRID ,0.0, 200.0) 

CALL VECONT (GRID, 0.0, 200.0, 1023.0,200.0) 

C 

C DISPLAY OF ARRAY GRID 
C 

CALL DIS°LY (2, GRID , DEC) 

C 

C DO YOU KNOW WHAT TYPE OF DATA IS IN THE FILE 
C YOU WANT TO ANALIZE? IF YES GO TO 121 
C 

85 TYPE 86 

86 FORMA T( ' DO YOU KNOW IF THE FILE YOU DESIRE TO 

1 ANALIZE CONTAINS’,/,' CALIBRATION OR FLIGHT DATA? 

2 Y OR N? ',i) 

ACCEPT 10, ANS 

I F(ANS .EG . *Y ' ) GO TO 121 
IFCANS.NE. 'N ') GO TO 85 
C 

C INPUT OF FILE # FOR THE INPUT OF DATA FROM THE 

C MAG TAPE 4 INITIALIZATION OF MAG TAPE 4 ARRAY ADVAL 

C 

IERCT :0 
TYPE 82 

82 FORMAT! ' TYPE THE » OF THE MAG TAPE FILE OF THE DATA 

1 WHICH IS',/,' TO BE OBSERVED ', S) 
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ACCEPT 30, I FIN 
C 

C INITIALIZATION OF MAG TAPE FOR INPUT ANO INPUT OF 
C ARRAY CONTAINING FILE IDENTIFICATION INFORMATION. 

C IDENTIFICATION INFORMATION IS PLACED INTO ARRAY 
C "FI LID". "FI LID < 1 )" CONTAINS THE NUMBER 
C OF THE MAG TAPE FILE. THIS NUM3ER IS COMPAIRED TO 

C "I FI N" TO SEE IF THE CORRECT FILE HAS BEEN READ. 

C IF AN ERROR OCCURS AN ERROR MESSAGE IS TYPED. 

C THE PROGRAM TRYES SIX TIMES TO READ IN THE FILE. 

C IF THE PROGRAM IS NOT ABLE TO READ IN THE FILE, 

C PROCESSING IS TERMINATED. 

C 

CALL MTINITQF1N, 1 ) 

CALL SETDEC 

81 CALL MTREADCFILIO, 12, 1 ERR) 

I F (I ERR .EQ .0) GO TO 87 

type 200 , ierr 

200 FORMAT ( ’ <--*MTl 0 ERROR ',11,' IGNORED PROCESSING 
ICONTI NUED ---»’,///) 

IERCT =1ERCT+ 1 

IF(I£RCT.LE.5) CALL MTSKIPC-1,0) 

IFUZRCT.LE.5) CALL KTSXIP(1,0) 

IFUERCT.LE.5) GO TO 81 
205 TYPE 210 

210 FORMAT < ’ ERROR WILL NOT CORRECT, PROCESSING STOPPED.’,////) 

GO TO 150 
C 

C TEST TO SEE IF "FlLIO(l)" CORRESPONDS TO "1FIN", 

C IF IT DOES PROCESSING IS CONTINUED THROUGH THE 
C PRGGRAM, Ir THEY ARE LOT EQUAL CORRECTIVE STEFS 
C ARE TAKEN TO FIND THE DESIRED MAG TAPE FILE. 

C 

87 IFF=FILID< 1 ) 

IF(IFF-IFIN) 1000,1010,1020 
1000 IERCT=IERCT+1 

1FCIERCT.LE.5) TYPE 200, IERR 
IFUERCT .GE.5) GO TO 205 
INFzIFIN-IFF 
CALL MTSK I P (I NF , 0 ) 

GO TO 81 

1020 IERCT :IERCT+ 1 

IFCIERCT.LE.5) TYPE 200, IERR 
I F (1 ERC T . GT . 5) GO TO 205 
I NF :I FI N-I FF* 1 
CALL” MTSK IP (I-NF,0> 

CALL MTSK I P ( 1 , 0 ) 

GO TO 81 
C 

C TYPE OUT OF FILE IDENTIFICATION INFORMATION 
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C AFTER THE CORRECT FILE HAS BEEN FOUND. 

C 

1010 TYPE 1100, FILID (2) , F1LIDC3), F1LID(4) , F1LIDC5), 

1 FILIDT6), FI LID ( 7) , FILIDC8), F1LIDO) 

1130 FORMAT ( * DATE OF LAUNCh: \2A5,/,* TYPE OF OA TAt 

1 ',2A 5, /, * TYPE OF DAY: \2A5,/ f * COMMENTS: \2A5, 

2 ////////) 

IERCT = 0 

CALL INTAB<ADVAL,3000) 

READ IN LOOP CONSISTING OF 1000 POINTS FOR DISPLAY 
OF ARRAY ADVAL WHICH CONTAINS DATA FROM A FILE 
0 MAG TAPE, IN ORDER TO DETERMINE IF THE FILE 
CONTAINS CALIBRATION OR FLIGHT DATA. 


GET DATA FROM THE MAG TAPE 

DO 90 K=l,50 
IERCT =0 

490 CALL MTREADITDMP, 100, IERR) 

IF (IERR ,EQ .2) GO TO 100 
IFCIERR .NE.0) GO TO 500 
IF(X.EQ.l) 1 NDX = 1 
DO 88 L = 0, 100,5 
LN=L 

I F(L.E9 .0) L N = 1 
FVAL(INDX)=TDMP(LN) 

88 I NDX =1 NDX+1 

GO TO 90 

500 TYPE 200, IERR 

IERCT =IERCT+ 1 

I F (IERCT .LE .5) CALL MTSKIP(0,-1) 

IFdERCT.LE.5) GO TO 490 
TYPE 210 
GO TO 160 
90 CONTINUE 

C 

C MULTIPLY DATA IN BY 150.0 & ADD 200.0 TO SHIFT ZERO 

C AXIS UP TO LEVEL 200.0 ON THE DISPLAY SCOPE 

C 

XI =0.0 

YI =FVAL( l )* 15O.0+200.0 
IF(YI .LT . 1023.0) GO TO 510 
YI =1020.0 
TYPE 91 
C 

C SETUP OF DISPLAY DATA ARRAY "ADVAL". IF ANY VALUES 
C ARE OVER MAXIMUM VALUE FOR DISPLAY THOSE VALUES 
C ARE TRUNCATED TO FIT THE DISPLAY, AND A MESSAGE 
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C IS TYPED OUT. 

C 

510 CALL POINT (ADVAL.XI ,YI) 

DO 515 K =2, 1000 
XE -K 

YE=FVAL(K)*153. 0+200.0 
IF(YE.LT.t32J.0> GO TO 89 
TYPE 91 

91 FORMAT? * «-— ERROR: VALUE TO BE DISPLAYED IS OUT OF 

1 THE RANGE OF DISPLAY —»*,/,• THE VALUE HAS BEEN 

2 TRUNCATED TO FI? IHE DISPLAY*,//) 

YEr 1020.0 

89 IF?YE-YI .LT. 0.999) CALL POI NT (ADVAL ,XI ,YI ) 

CALL LI NE (ADVAL , XI, YI, XE, YE) 

XI =XE 
515 YI=Y£ 

C 

C DISPLAY OF DATA 1 GRID ARRAYS 
C & ZERO THE ISXIP FLAG 
C 

100 CALL DISPLY (3 , GRID, ADVAL, DEC) 

IFCIERR .EQ.2) TYPE 105 

105 FORMAT? * ***END OF MAG TAPE FILE***’,///) 

ISKIP =0 
C 

C ENABLING OF LITE PEN FOR DECISION MAKING PROCESS 
C 

110 CALL LI TEPN ?IDET , XDET , YDET) 

IF (IDET.EQ.0) GO TO 110 

IF? (XDET. GT.680.0). AND. (XDET. LT. 800.0). AND. ( 
1YOET.LE.A0.0)) GO TO 120 

IF((XDET.GT. 920.0). AND. (YDET.LE. 40. 0)> GO TO 150 
C 

C IF CAL DETECTED FROM DISPLAY SCOPE, CHAIN 
C TO LINK CALC (CALIBRATION COMPUTATION) 

C 

I F( (XDET .LE .0) .OR . (XDET .GE . 1 80 ) .OR • (YDET .GE . 40) ) 

1 GO TO 112 
NUM :6 

CALL CHAIN(0, IDSK, 'CALC*) 

C 

C IF DAT DETECTED FROM DISPLAY SCOPE, CHAIN 
C TO LINK CALD (CALCULATIONS ON FLIGHT DATA) 

C 

1 12 I F( (XDET. LE. 180). OR .(XDET .GE .380) .OR . (YDET .GE . 40) ) 

l GO TO 115 
NUM = 7 

CALL CHAIN(0,IDSX, ’CALD’) 

C 

C FLAG ISKIP TRANSFERS CONTROL TO 125 IF THE FLAG 
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C IS SET By LINK "CALC" OR LINK "CALD". 

C 

115 I F (ISK IP .EG . 1 ) GO TO 125 

GO TO 110 

120 CALL D1SPLY (0) 

GO TO 1010 

C 

C ZERO THE DISPLAY 
C 

160 CALL DISPLY(0> 

C 

C IF THERE HAS BEEN AN ERROR WHICH WILL NOT CORRECT 

C THERE IS A PAUSE AND ONE CAN EITHER CONTINUE 

C OR EXIT THE PROGRAM. 

C 

I F(I ERCT .GT . 5 ) PAUSE 'WHAT DO YOU WANT TO DO?' 

GO TO 125 
C 

C OPERATOR PROVIDES TYPE OF FILE TO BE ANALIZED 
C AND THE CORRECT SUBROUTINE IS CALLED. 

C 

121 ISK IP =0 
TYPE 122 

122 FORMATC IF THE FILE YOU WISH TO ANALIZE IS A 

1 CALIBRATION FILE TYPE CAL',/,' IF IT IS A FLIGHT 

2 DATA FILE TYPE FLT ',$> 

ACCEPT 10, TYPF 

IF (TYPF.NE. 'CAL ') GO TO 123 
NUM=4 
C 

C CHAIN TO LINK "CALC" IF ANSWER TO ABOVE QUESTION 
C IS "CAL". 

C 

CALL CHAIN(0,IDSK, ’CALC") 

119 IF (ISK1P .EG . 1 ) GO TO 125 

123 IFCTYPF.NE. 'FLT ’) GO TO 124 
NUM = 5 

C 

C CHAIN TO LINK "CALD" IF ANSWER TO ABOVE QUESTION 
C IS "FLT". 

C 

CALL CHAIN(0,IDSK , 'CALD ') 

124 IF(ISKIP.NE.l) GO TO 121 
C 

C DECISION TO COMPUTE CONDUCTIVITIES 
C 

125 TYPE 126 

126 FORMATC AT THIS POINT DO YOU WISH TO COMPUTE 

l CONDUCTIVITIES? Y OR N $) 

ACCEPT 10, COND 
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IF(COND.EQ . ‘N ') GO TO 150 
I F CCONO . NE . *Y ' ) GO TO 125 
C 

C LINK TO CONC (CONDUCTIVITY CALCULATION) 

C 

155 NUM:3 

CALL CHA1N(0,IDSK , ‘CONC *> 

C 

C TYPE QUESTION AS TO WHETHER OPERATI OR WANTS 
C TO MATCH TIMES TO CORRESPONDING ALTITUDES. 

C 

156 TYPE 7 
ACCEPT 10, MTA 

I F (MTA .£Q . 'N ’ ) GO TO 150 
IF(MTA.N£. *Y *> GO TO 156 

157 IF(IDSK.EQ.l) GO TO 2000 

C 

C GO TO 2000 IF OPERATING UNDER THE DISK, OTHERWISE 
C TYPE THE FOLLOWING STATEMENT. 

C 

type iss 

158 FORMAT ( * TYPE » OF DTA ON WHICH POSITIVE CONDUCTIVITY 

1 IS STORED 4’,/, * I OF DTA ON WHICH NEGATIVE CONDUC 
2TIVITY IS STORED. *,//, 1 POSITIVE: \J> 

ACCEPT 30, IPSI 
TYPE 159 

159 FORMAT ( * NEGATIVE: *,S> 

ACCEPT 30, INSI 

GO TO 2010 
2000 IPSI =IDSK 

INSI rIDSK 

2010 IF(IPLT.EQ.I) GO TO 175 

NUM = 1 
C 

C CHAIN TO LINK "MTAH" TO MATCH TIMES TO ALTITUDES. 

C 

CALL CHAIN(0,IDSK, 'MTAH ’) 

C 

C QUESTION TO ASCERTAIN WHETHER OPERATOR DESIRES 
C TO PLOT CONDUCTIVITIES ON THE DISPLAY. 

C 

165 TYPE 170 

170 FORMAT ( * DO YOU DESIRE TO PLOT ALTITUDE VS 

1 CONDUCTIVITY CURVES? Y OR N, *,$) 

ACCEPT 10, ANS 
IF<ANS ,EQ . 'N ’) GO TO 150 
IF CANS. ME'. ’Y *) GO TO 165 
175 NUM=2 

C 

C CHAIN TO LINK "SIGP" TO PLOT ALTITUDE VS. 
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C CONDUCTIVITY. 

C 

CALL CHAIN(0,IDSK, 'SIGP *) 

C 

C TYPE OUT THE FINAL QUESTION TO BE ASKED IN THE 
C PROGRAM. 

C 

150 TYPE >30 

130 FORMAT ( * ARE YOU FINISHED WITH THE PROGRAM? Y OR N. 

ACCEPT 10, ANS 
IF (ANS ,EQ . 'N 1 ) GO TO 6 
I F (ANS . NE . 'Y *) GO TO 150 
END 
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C ************************************************************** 

LINK CALCs DETERMINES THE SLOPES OF 
CALIBRATION DATA 

************************************************************** 

COMMON LP, ISHOW, XMAX, XMIN, YMAX, YMIN, INTENS, ISCALE 
COMMON /A/ ADVAL (3.100 } , GR1DCI0I20) /C / AVCASL /D/ STARC20), 

1 SEC, (10) /E/ SLOPE (150,2) /F/X(4),Y(4) /P / ILINE ( 100) , 

2 RB ( 5) /G/ ISKIP /H/ T1F(150,3) /I/ TIME(250) /M/ 

5 I ND, I FLAG , IE0F, KIN, I HDF, NFLG /T/ CVALO600) 

4 /DA/ IDSK /DC/ FI LI D ( 1 0) ,TDMP ( 1 00 ) 

STATEMENT REQUIRING OPERATOR ACTION UPON 
ENTERING SUBROUTINE. 


PAUSE ’CALIBRATION SLOPE CALCULATIONS’ 

AVCASL=0.0 

ADDITION OF THE LITERALS ’SR’,’SL’,4 ’TAKE SLOPE’ WHERE 
SR IS SHI FT@RI GHT & SL IS SHIFT LEFT. 

CALL CHRGEN(GRIl), ’SR ’ ,210.0,62.0,2,0) 

CALL CHRGEN(GR1D, ’SL ’,260.0,62.0,2,0) 

CALL CHRGENCGRID, ’TAKE SLOPE *,680. 0,55. 0,3,0) 

REQUESTS THE # OF THE CALIBRATION FILE BEING 
ANALIZED 4 INITIALIZES NAG TAPE FOR INPUT 

IERCT =0 
NFLGrl 

REQUEST FOR INPUT OF MAG TAPE FILE NUMBER OF THE 
CALIBRATION FILE TO BE ANALIZED. 

1 TYPE 2 

2 FORMA T ( ' TYPE THE MAG TAPE FILE # OF THE CAL FILE YOU 
1 ARE ANALIZING. ’,$) 

ACCEPT 3, I FILE 

3 FORMAT(I) 

TYPE 5000 

5300 FORMA T ( * WHJLE TAK I NG CAL SLOPES ESTIMATE THE TIME 

T REQUIRED',/, ’ FOR SLOPE TO SWEEP FROM MINIMUM 
2T0 MAXIMUM. AT THIS TIME TYPE IN SWEEP VOLTAGE. ’,$) 
ACCEPT 5100, SWVOL 
5100 FORMAT (F 10,5) 

C 

C I NITCALIZATI ON OF MAG TAPE, SETTING DEC MODE 4 
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C READ IN OF FILE IDENTIFICATION INFORMATION. 

C IF THERE IS AN ERROR IN READING IN ThE IDENTIFICATION 

C FILE A MESSAGE IS TYPED & ThE PROGRAM WILL ATTEMPT 

C TO TRY AGAIN. IF THE ERROR DOES NOT CORRECT BY THE 

C SIXTH ATTEMPT THE PROGRAM WILL TRANSFER CONTROL TO 

C LINK ’MAIN*. 

C 

CALL MTINlTdFILE.l) 

CALL SETDEC 

6 CALL MTREAD < FI LID, 13.IERR) 

I F (IERR . EQ .0 ) GO TO 7 
TYPE I I, IERR 

II FORMAT ( * <--- MTI 0 ERROR *,11,* IGNORED PROCESSING 

1 CONTINUED --->',///) 

IERCT :IERCT+ l 

I F (IERC T .LE . 5 ) CALL MTSKIPC-1,0) 

I F (IERCT ,LF. . 5 ) CALL MTSKIP<1,0> 

IFCIERCT .LE.5) GO TO 8 

6 TYPE 9 

9 FORMAT C * ERROR WILL NOT CORRECT, PROCESSING STOPPED. 

GO TO 300 
C 

C SECTION TO DETERMINE IF CORRECT FILE HAS BEEN 
C LOCATED. 'FILID ( 1 ) ' CONTAINS THE » OF 
C THE DESIRED FILE. IF THESE TWO NUMBERS ARE EQUAL 
C THE PROGRAM BEGINS DATA ANALYSIS, OTHERWISE, 

C CORRECTIVE ACTION TO FIND THE CORRECT FILE IS 
C TAKEN. THIS IS DONE BY EITHER ADVANCING 'IFF-IFILE ’ 

C FILES OR BACKSPACING ’IFILE-IFF-1 ' FILES AND READING 

C THE NEW FILE LOCATED. 

C 

7 IFFrFILIDd) 

I FdFF-1 FILE 1500,510,520 
500 IERCT =IERCT+ 1 

I F (IERCT .LE , 5) TYPE 11, IERR 
lFdERCT.GT.5) GO TO S 
1 NF=I FI LE -I FF 
CALL MTSKIP (INF,0) 

GO TO 8 

520 IERCT :I ERCT+ 1 

I F(IERCT .LE . 5) TYPE II, IERR 

IFdERCT .GT.5) GO TO 6 
I nF =1 FILE -I FF - 1 
CALL MTSKIP (I NF, 0) 

CALL MTSKIP Cl ,0) 

GO TO 8 
C 

C TYPE OUT OF FILE IDENTIFICATION INFORMATION 
C 

510 TYPE 515, FILI D (2 ) , FILID<3), FILID<4), FILID<5), 


,////) 


* 
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l FILID(6),FILID(7),FILID<8>,FILID<9> 

515 FORMATC DAIE OF LAUNCH: \2A5,/,* TYPE OF DATA \ 

12A5,/, • TYPE OF DAY: \2A5,/, * COMMENTS: \2A5, 

2 ////////) 

I ERC T =0 
I FLAGsfl 
I ND = 1 
INDFrfl 
IEOF:0 
N = 0 

DO LOOP TO READ IN 1600 DATA POINTS INTO ARRAY *CVAL*. 

DO 20 L=l,80 
IERCT =0 

READ IN OF ONE BLOCK OF DATA FROM MAG TAPE. 

FOLLOWING ’MTREAD* THERE IS AH ERROR DETEDTION SECTION. 


400 CALL MTREAD(TDMP, 100, IERR) 

IFCIERR .EQ.2) TYPE 405 

405 FORMAT < ' *** END OF MAG TAPE FILE ***’,//> 

IF (IERR .EG .2) GO TO 40 
IF (IERR .NE.0) GO TO 410 
IF(L.EQ.I) INDX=1 
DO 5 K=0, 100,5 
KN=K 

IFCK.EQ.0) GO TO 5 
C 

C ADDITION CF SCALINC FACTORS TO EACH ELEMENT OF THE 
C BLOCK OF DATA READ IN FROM MAG TAPE & THEN EACH OF 

C THESE ELEMENTS BECOME ONE ELEMENT OF ARRAY ’CVAL*. 

C FOLLOWING TmIS EACH ELEMENT IS TESTED TO INSURE IT DOES 
C NOT EXCEED DISPLAY SCOPE RANGE. IF A VALUE 
C EXCEEDS THE RANGE IT IS TRUNCATED TO FIT THE DISPLAY 

C 4 A MESSAGE IS TYPED PRIOR TO ACTUAL DISPLAY OF DATA 

C 

CVAL(INDX)=TDMP<KN)* 150. 0+200.0 
IF(CVALCINDX).LE. 1023.0) GO TO 5111 
IF(N.NE.0) GO TO 10 
TYPE 31 
N:1 

10 CVAL(INDX)=1020.0 

5111 IN0X=INDX+1 

i 5 CONTINUE 

GO TO 20 
C 

C SHOULD AN ERROR OCCUR DURING READ IN OF A BLOCK OF 
C DATA. CONTROL IS TRANSFERRED HERE TO ATTEMPT TO CLEAR 
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C THE ERROR 
C 

410 TYPE 11,1 ERR 
IERCT :I ERCT+ 1 

IFUERCT.LE.5) CALL MTSKIPC0.-1) 

I F (IERCT .LE .5) GO TO 400 
TYPE 9 
GO TO 75 

31 FORMAT < ' <-— ERROR: VALUE TO BE DISPLAYED IS OUT OF 

1 THE RANGE OF DISPLAY ---»',/,' THE VALUE HAS BEEN 
2TRUNCATED TO FIT THE DISPLAY',//) 

20 CONTINUE 
C 

C TEST TO SEE IF 'CVAL* CONTAINS 1600 ELEMENTS IN ORDER 
C TO CENTER THE DISPLAY OF DATA CONTAINED IN ARRAY 'CVAL'. 

C 

40 L=1NDX-1 

IF(L.EQ.1600) GO TO 25 
INDF:(L-1000)/2 
IF(INDF) 21,22,23 

21 INDFil 
I =L 

GO TO 26 

22 I NDF : 1 
1:1000 
GO TO 26 

23 I:L-INDF 
GO TO 26 

25 I NDF :300 
I :L-3U0 

C 

C INITIALIZATION AND SET UP OF DATA DISPLAY ARRAY 'ADVAL*. 

C 

26 CALL INTAB(ADVAL,3000) 

XI =0.0 

YI =CVAL(INOF) 

CALL POINT (ADVAL, XI ,YI ) 

DO 30 K :I NDF+ 1 , 1 
XE:XI+1.0 
YE=CVAL(K ) 

IF(YE-YI .LT. 0.999) CALL POI NT (ADVAL ,XI ,YI ) 

CALL LINE(ADVAL,XI ,YI ,XE,YE) 

XI sXE 
30 YI =YE 

C 

C CALL TO SUBROUTINE 'SLOPCA ' TO ANALIZE 
C CALIBRATION DATA. 

C 

CALL SLOPCA 

IF (1EOF.EQ.D GO TO 75 
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GO TO 4 
C 

C OPTION TO CONTINUE PROCESSING IF A MINOR ERROR 
C OCCURS OR EXIT IF A MAJOR ERROR OCCURS. 

C 

75 IF (IERCT .GT .5) PAUSE *WHA T DO YOU WANT TO DO?’ 

C 

C TRANSFER OF CONTROL TO STATEMENT 1203 IF THE PROGRAM 
C IS OPERATING OFF THE DISK SYSTEM; OTHERWISE, THE PROGRAM 
C REQUESTS DATA STORAGE INFORMATION. 

C 

TYPE 5200 

5200 FORMAT ( * TYPE IN ESTIMATED TIME FOR PROBE TO SWEEP 

1 FROM * , /, * FROM MINIMUM TO MAXIMUM. *,$> 

ACCEPT 130, ST 
VDS: SWVOL/ST 
IFdDSK .E8.1) GO TO 1000 
TYPE 76 

76 FORMAT C TYPE THE # OF THE DTA YOU WISH THE DATA TO BE* 
1,/,* STORED. *,S) 

ACCEPT 3, IOTA 
I DTA =1 DTA+8 
GO TO 1010 
C 

C INITIALIZE STORAGE DEVICE TO ACCEPT CALIBRATION SLOPE DATA 
C 

1000 IDTArlDSK 

C 

C OUTPUT OF CALIBRATION SLOPES & AVERAGE CALIBRATION 
C SLOPE TO STORAGE DEVICE. 

C 

1010 CALL OFILEdDTA, *VOST *) 

WRITEdDTA , 130) VDS 
CALL RELEAS (I DTA ) 

CALL OFILEdDTA, *CAL ’) 

DO 90 I =1 ,IND-l 

90 WRITEdDTA, 130) SLOPEd , 1 ) 

CALL RELEAS dDTA) 

DO 85 I = 1 , I ND - 1 

85 AVCASL: A VCASL+SLOPE (I , 1 ) 

AVCASL=AVCASL/dND-l ) 

C 

C INITIALIZE STORAGE DEVICE TO ACCEPT AVERAGE OF CAL SLOPES 
C 

CALL OFILEdDTA, *AVCAS *) 

WRITEdDTA, 130) AVCASL 
CALL RELEAS (I DTA ) 

IFdDTA .EQ.IDSK) GO TO 1020 
IDTA =10TA-8 
C 
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C PRINT OUT OF HEADINGS AND FILE IDENTIFICATION 
C INFORMATION 
C 

PRINT 78 

78 FORMAT < ’ 1 THE CALIBRATION DATA HAS BEEN 

1 STORED ON DTA'.IX,' UNDER THE NAME CAL. DAT.’) 

PRINT 79, IOTA 

79 FORMAT ( '+ ' , 44X , I l ) 

PRINT 1315, FI LID (2) , FILID<3), FILI D ( 6) ,F1 LID ( 7) , 
1FILID(8), FILID (9) 

1015 FORMAT ( *- * ,22X , ’DATE OF LAUNCH t *,2A5,5X, 'TYPE OF DAY: 

1 ', 2A 5, 5X, ’COMMENTS: \2A5) 

GO TO 1040 

1020 PRINT 1230, FI LI D (2 ) , FILID(3), FILID(6>, FILID(7), 
IFILID(S), FILID(9> 

1030 FORMA T C ’ 1 THE CALIBRATION DATA HAS BEEN STORED ON THE 
IDSK UNDER THE NAME CAL. DAT . ’ , / , ’,20X, ’DATE OF THE 

2 LAUNCH *,2A5,5X, 'TYPE OF DAY: * ,2A 5, 5X , 'COMMENTS : 

3 ,2A 5) 

C 

C PRINT THE VALUES OF THE CALIBRATION SLOPES 
C 

1040 PRINT 80 

80 FORiVATC- THE VALUES OF THE SLOPES OF THE CALIBRATE 
1 CURVES FOLLOWS:*,/) 

DO 9 l I=l,IND-l 

91 PRINT 95, I, SLOPE (1 , 1 ) 

95 FORMAT ( '* SLOPE <’,I3, ')= ', F10.5) 

130 FORMAT (F10.5) 

C 

C PRINT THE AVERAGE CALIBRATION SLOPES 

C 

PRINT 120, AVCASL 

120 FORMAT <\ THE AVERAGE CAL SLOPE: ', F10.5) 

300 ISKI P = 1 

IFCIDTA .EQ.IDSK ) GO TO 1050 
I DTA :I DTA + 8 
1050 CALL RELEAS <3 ) 

CALL DISPLY<0> 

C 

C LINK BACK TO 'MAIN* UPON COMPLETION OF THE ANALYSIS 
C OF CALIBRATION DATA . 

C 

CALL ChaI N (0,1D5K , 'MAI N ’ ) 

END 


C **************************** ****** **^* ************** ******** 
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C 

C SUBROUTINE SLOPCA : CALCULATION OF SLOPES OF 

C CALIBRATION FLIGHT DATA. 

C 

C *********************************** ******************** ****** 

SUBROUTINE SLOPCA 

COMMON LP, ISHOW, XMAX, XMIN, YtfAX, YMIN, INTENS, I SCALE 
COMMON /A/ ADVAL(3233) , GRIDO023) /C/AVCASL /D/SIAR(20), 

1 SEC C 1 2 ) /E/ SLOPE( 153,2) /F / X(4), Y ( 4) /P/ ILINE(100), 

2 RB (5) /G/ ISXIP /M/ IND, I FLAG, IEOF, KIN, INDF, NFLG 

3 /T/ C VAL( 1 600) 

C 

C SET UP LITERAL SHOWN FOR DISPLAY ON SCOPE, TO ALLOW 
C USER TO ANAL1ZE CALIBRATE SLOPES. 

C 

ISTAT s0 
JFlP =0 

IF(NFLG.GT.l) GO TO 5 
NFLGsNFLG+1 

CALL CHRGEN(GRID, ’NEXT SET *,300.0,62.0,2,0) 

C 

C INITIALIZATION OF ARRY STAR A SEC. THEN CONSTRUCT 
C MOVEABLE STAR. 

C 

5 CALL I NTAB (STAR, 20) 

CALL I NTAB (SEC , 10) 

ISHOW:0 

CALL POINT(SEC, 155.0,75.0) 

CALL LINE (STAR, 20. 0,20. 0,0. 0,20.0) 

I SHOWs! 

CALL LINE (STAR, 0.0, 20. 0,40.0,20.0) 

ISHOU;0 

CALL LINE (STAR, 40. 0,20. 0,40. 0,0.0) 

ISKOWsl 

CALL LINE (STAR, 40. 0,0. 3, 0.3, 40.0) 

ISHOWsO 

CALL LINE (STAR, 0.0, 40. 0,20. 0,40.0) 

ISHOWsl 

CALL LINE (STAR, 20.0, 40.0,20.0,0.0) 

ISHOW=0 

CALL LINE (STAR, 20. 0,0. 0,0. 0,0.0) 

ISHOWsl 

CALL LINE (STAR, 0.0, 0.0, 40. 0,40.0) 

C - - ------ - - - - 

C MPOINT MOBILIZES THE STAR, THEN ALL ARRAYS ARE 
C DISPLAYED, THEN CALL LITEPN TO ACTIVATE THE LITE PEN. 

C 

CALL MPOI NT (SEC , 155.0,75.0) 

CALL DISPLY (4, SEC, STAR, GRID, AOVAL) 
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10 CALL LI TEPN (NSLOPE , XS , YS ) 

IF ( NSLOPE. EQ .0) GO TO 10 

C 

C DECISION SECTION: 

C REJECT???? 

C 

IFUXS . GT. 683.0). AND. (XS .LT . 800 . 0) .AND . (YS .LE . 40. 0) ) 

1 GO 10 50 
C 

C END???? 

C 

I F ( (XS . GT. 920.0). AND ,<XS .LT . 1 023 .0 > .AND . (,YS .LE . 40.0) > 

1 GO TO 60 
C 

C NEXT SET???? 

C 

I F ( (XS .GT. 300.0) .AND .(XS .LT .380 . 0) .AND . (YS .GT . 62.0) .AND . 
1 <YS .LT .76.0)) GO TO 50 
IF (INDF.EQ . 1 ) GO TO 19 
C 

C SHIFT RIGHT??? 

C 

I F( (JMP .EQ . 1 ) .AND . (ISTAT .EQ . 1 ) ) GO TO 12 
IF((XS.LT.210.0).OR.(XS .GT .230.0) .OR . (YS .LT . 62. 0) .OR. 

1 (YS .GT .76.0)) GO TO 12 
• IFGSTAT.NE.2) GO TO 11 
IL=INDF 
IM=INDF+1000 
ISTATI0 
GO TO 1000 

11 !L=1 

IK : 1 030 
JMP s 1 
ISTAT=l 
GO TO 1000 
C 

C SHIFT LEFT??? 

C 

12 I F ( (XS .LT. 263.0) .OR .(X3 .GT .280.0) .OR . (YS .LT .62.0) .OR. 


l (YS .GT . 76. 

0)) 

GO 

TO 

19 

IFGSTAT.EQ 

.2) 

GO 

TO 

10 

IF (ISTAT. EQ 

.0) 

GO 

TO 

13 


ILslNOF 
IM=INDF+ 1000 
XSTAT:0 
JMP =0 
GO TO 1000 
13 IL=2*INQF 

IK=IL+1000 
ISTAT -2 
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120k) XI =2.2 

CALL 0I5PLY (2) 

INITIALIZATION OF DATA DISPLAY ARRAY 'ADVAL * AND 
BUILD UP OF SAKE TO DISPLAY DATA AFTER A SHIFT LEFT 
OR SHIFT RIGHT. 

CALL IN TAB (A OVAL, 3222) 

Yl rCVALCIL) 

CALL POI NT (ADVAL ,X I ,YI ) 

DO 1213 X=IL+l,I-'!-l 
XE=XI+I .0 
YErCVAL(K) 

1FCYE-Y1 .LT .0.999) CALL POINT CAD VAL ,XI ,YI ) 

CALL LINECADVAL.XI ,YI ,XE,YE) 

XI =XE 
010 YI=YE 

DISPLAY OF ALL ARRAYS TO SHOW SHIFTED SECTION 
OF DATA ARRAY. 

CALL DISPLYC 4, SEC .STAR, GRID, ADVAL) 

GO TO 10 

CALCULATION OF SLOPES 

9 1 F< (XS . LT. 680.2 ). OR. (XS.GT. 830.0). OR.CYS.LT. 55.0). OR. 

1CYS.GT.76.0)) GO TO 10 
ISA V: 1 

ACTIVATE LITEPEN TO ACCEPT FIRST POINI IN CALCULATING 
SLOPES • 

IF POINT RECEIVED INDICATES "ACCEPT” GO TO 30 
OTHERWISE MOVE 1HE STAR. AND INITIALIZE 
ARRAY ILINE. 

I CALL LITEPN(IACC,X2,Y2) 

I F (I ACC .EG.0) GO TO 21 

IF((X2. GT. 380. 0).AND.(X2.LT. 630.0). AND. (Y2.LE. 42.0)) GO TO 30 
CALL MPOI NT (SEC , X2, Y2) 

CALL DISPLYC4, SEC, STAR, GRID, ADVAL) 

CALL I NTAB (ILI N£ , l 00 ) 

X A Y COORDINATES OF DESIRED POINT ON CURVE 

X (ISAV) =X2 
Y (ISAV) =Y2 

CALL POI NT ( ILI NE , X(l), Y(l)> 

IFdSAV.LT. 2) GO TO 21 
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C CHECK TO SEE IF THERE WAS AN ERROR IN SELECTING 
C THE SECOND POINT ON THE CURVE . IF THERE WAS 
C DISPLAY THE ERROR MESSAGE SHOWN, CONNECT THE 
C TWO POINTS SELECTED WITH A LINE, AND DISPLAY 
C ALL ARRAYS. IF NO ERROR JUST CONNECT THE 

C POINTS WITH A LINE A DISPLAY ALL ARRAYS. 

C 

CALL LINE <1 LINE, XU), YU), X<2), Y<2>) 
IFC(YU).LT.Y(2)).AND.(X(1).LI.X(2))> GO TO 25 
C 

C SETUP A DISPLAY OF ERROR MESSAGE IF A NEGATIVE 
C SLOPE IS DETECTED. 

C 

CALL CHRGEN CILI NE , ERROR NEG. SLOPE --*’,210.0,82.0,2,3) 
CALL DISPLY (0) 

DO 26 I E =1 ,4 

CALL DISPLY < 1 ,ILINE> 

CALL WAIT (250 > 

26 CALL DISPLY (2) 

25 CALL DISPLYC5, SEC, STAR, GRID, ADVAL.ILINE) 

GO TO 21 
C 

C WAIT 1 500 MSEC, CLEAR LITE PEN VALUES, A MOVE STAR 
C TO READY POSITION. INCREASE ISAV BY l, A IF ISAV« 

C OR -2 GO TO 50, TO ACCEPT SECOND COORDINATES OF DATA 
C IF ISAV>2 CALCULATE SLOPES. 

C 

30 CALL WAITU500) 

CALL LITEPNUCLE, X4, Y4) 

CALL MPOI NT (SEC , 155.0, 75.0) 

ISAVrISAV+l 

IFUSAV.LE.2) GO TO 21 

SLOPE (I ND, I ): (Y<2)-Y<l))/(XC2)-XCl)> 

IND = IND-t-1 
GO TO 10 
53 IFLAG=0 

GO TO 70 
C 

C UPON COMPLETION OF CALIBRATION DATA ANALYSIS FLAG *IEOF’ IS 
C SET TO ONE WHEN ’END’ IS DETECTED BY DISPLAY. 

C 

60 I EOF : 1 

70 CALL DISPLY (0) 

RETURN 

END 
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C j************************************************************** 

CALD: DETERMINES THE SLOPES OF 
THE FLIGHT DATA 

** *** ******** ******** ************************************* 

COMMON LP, ISH3W.XMAX, XMIN, YMAX, YM1N, INTENS, ISCALE 
COMMON /A/ ADVAL <3223) , GRI D ( I 000) ,DEC ( 1 00 ) 

4 /D / STAR (22 ) , SEC(IO) /DC/ FILID<I0), TDMPU00) 

1 /E/ SLOPE ( 1 53 , 2 ) /F / X(4), Y(4> /G/ 1SKIP. /H/ TIF(150,3> 

2 /I/ TIMEC253) /J / ASK(4), R(6>, B(2) /M/ I ND , I FLAG, 

3 I EOF, KIN.INDF, NFLG /S / FVALU005) /DA/ IDSK 

SET UP THE LITERAL SHOWN TO DETERMINE THE TYPE 
(POS. OR NEG. COND.) OF SLOPE TO BE CALCULATED. 


STATEMENT REQUIREING OPERATOR ACTION UPON 
ENTERING SUBROUTINE. 

PAUSE ’FLIGHT SLOPE CALCULATIONS ’ 

CALL DISPLY(I.GRID) 

NFLG=1 

SET UP THE LITERAL SHOWN TO DETERMINE THE TYPE 
(POS. OR NEG. COND.) OF SLOPE TO BE CALCULATED. 

CALL CHRGEN(GRID, ’SLOPE FOR +OR- C OND ??’, 635 . 0 , 82 . 0 , 2 , 0 ) 
CALL CHE GEN (GRID , ’+', 650.0, 44.0, 5 , 0) 

CALL CHRGEN (GRI D , ’-’,978.0,44.0,5,0) 

ADDITION OF LITERAL TO ARRAY ’GRID’ TO ALLOW MAGNIFICATION 
BY A FACTOR OF 4 & RETURN TO NORMAL DISPLAY. FOLLOWING 
THAT THERE IS INSERTION OF LITERAL SEPARATORS. 

CALL CHRGENCGR1D, *X4 NORMAL ’ , 71 0 .0, 55. 0 ,3 , 0) 

DO 2 J : 1 ,2 

CALL POI NT(GRID,B (J) , 40.0) 

CALL LINE (GRID, 3 (J),40.0,B (J),80.0) 

INITIALIZATION 4 SET UP OF ARRAY ’TIME’ FOR A TEN 
SECOND TIME BASE. 

CALL I NTAB (TI MF , 250) 

XCOOR :65 . 0 
DO 1 1:1,10 
ENC0DE(4, 5, ASK ) , I 
5 FORMAT (14) 
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ASK (2):0.0 

CALL CHRGEN(T1ME ,ASK ,XCOOR ,130.0,2,0) 

1 XCOOR rXCCOR+100.0 

CALL CHRGENIGRID, 'TIME I N * , 900 . 0, I 1 5 . 0 ,2 , 0 ) 

CALL CHRGEN(GRID, * SEC *,900.2,100.0,2,0) 

ACCEPTS THE NAG TAPE PILE # OP THE PLIGHT DATA 
BEING ANALIZEO, S. INITIALIZES MAG TAPE FOR INPUT. 

IERCT=3 
TYPE 20 

20 FORMAT ( * TYPE THE « OF THE MAG TAPE FILE t. WHICH CONTAINS 
1 THE*,/,* FLIGHT DATA YOU ARE ANALYIZING. *,S) 

ACCEPT 25, I FILE 

25 FORMAT (I ) 

CALL KTI N I T (I FILE , 1 ) 

CALL SETDEC 
C 

C READ IN OF FILE IDENTIFICATION INFORMATION. 

C IF THERE IS AN ERROR IN READING IN THE INFORMATION 
C A MESSAGE IS TYPED & THE PROGRAM WILL ATTEMPT 
C TO READ IT AGAIN. SHOULD THE ERROR NOT CORRECT BY 
C THE SIXTH ATTEMPT THE PROGRAM WILL TRANSFER 
C CONTROL TO LINK 'MAIN *. 

C 

19 CALL MTREAD (FI LID , 1 3 ,1 ERR ) 

I F (I ERR .LE . 0 ) GO TO 24 
TYPE 21,1 ERR 

21 FORMAT < ' MTI 0 ERROR Ml,* IGNORED PROCESSING 

1 CONTINUED --->*,///) 

IERCT ;IFRCT+ l 

1FUEPCT.LE.5) CALL MTSKIPC-1,0) 

I F( I ERCT .LE . 5 ) CALL MTSK1P(1,0> 

IF (IERCT .LE. 5) GO TO 19 
4 TYPE 22 

22 FORMAT ( * ERROR WILL NOT CORRECT, PROCESSING STOPPED.*,////) 
GO TO 290 

C 

C SECTION TO DETERMINE IF CORRECT FILE HAS BEEN 
C LOCATED. *FI LI D ( 1 ) * CONTAINS THE # OF THE MAG TAPE 
C FILE THAT HAS BEEN FOUND, WhILE 'IFILE* IS THE » OF 

C THE DESIRED FILE. IF THESE TWO NUMEERS ARE EQUAL 

C THE PROGRAM BEGINS DATA ANALYSIS; OTHERWISE, 

C CORRECTIVE ACTION TO FIND THE CORRECT FILE IS 
C TAKEN. THIS IS DONE BY EITHER ADVANCING *IFF-IFILE* 

C FILES OR BACKSPACING *1 FILE -I FF- 1 ’ FILES AND READING 
C THE NEW FILE LOCATED. 

C 

24 I FF:FI LI D ( 1 ) 

3 F (I FF-I FILE ) 1001,1010,1020 
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1031 IERCT =1 ERCT+ 1 

1 F (IERC T .LE . 5) TYPE 21.IERR 
IFtIERCT.GT.5) GO TO 4 
I NFsIFILE-I FF 
CALL MTSXIP (INF f 0) 

GO TO 19 

1020 IERCT sI£RCT+ I 

IFdERCT .LE.5) TYPE 21.IERR 
IF (IERCT .GT .5) GO TO 4 
INFsIFIN-IFF-l 
CALL MTSKIP (INF ,0)' 

CALL MTSKIP (1,0) 

GO TO 19 

TYPE OUT OF FILE IDENTIFICATION INFORMATION 

010 TYPE I030,FILID(2),FILID(3),F1LID(4>,FILID(5), 

1 FILID(6),FILID(7),FILID(8) ,FILID<9) 

030 FORMAT ( * DATE OF LAUNCH: *,2A5,/, * TYPE OF DATA: \2A5,/, 

l' TYPE OF DAY: ',2A5,/,* COMMENTS: \2A5, ////////) 

IERCT =0 
IFLAGs0 
I ND: 1 
I NDF: 1 
IEOF:0 
KINsl 
MI Ns 1 0 

INITIALIZATION OF ARRAY DEC TO FACILITATE 
THE 2 SEC DISPLAY OPTION USED WHEN ANALYZING 
NEGATIVE CONDUCTIVITY SLOPE. IT OCCURS TWICE 
BECAUSE AFTER THE FIRST DISPLAY OF DATA THE 
PROGRAM NEVER PASSES THE FIRST POINT. 

CALL I NTAB (DEC ,100) 

CALL ChRGEN (DEC , *2 SEC DIS *,20.0, 15.0,2,0) 

GO TO 27 

6 CALL INTAB(DEC, 100) 

CALL CHRGEN (DEC , '2 SEC DIS ' ,20. 0, 1 5.0, 2, 0) 

UP DATING OF TIME BASE DATA ON 
SUCESSIVE DISPLAYS BY INITIALIZING AND 
CHANGING ARRAY 'TIME '. 

CALL I NTAB (TIME , 250) 

XCOOR s65.0 
KINsKI N+ 10 
MI N=MIN+ 10 
DO 16 IsXIN.MIN 
ENCODE (4.15 .ASK ) , 1 
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15 FORMAT (14) 

ASK (2) =0.0 

CALL CHRGENCTIME ,ASK ,XCOOR, 130.14,2,0) 

16 XCOOR=XCOOR+ 100.0 

27 CALL DISPLY<3,GRID, TIME ,DEC ) 

IDISER =0 
INDX=l 
C 

C DO LOOP TO TRANSFER DATA FROM MAG TAPE 
C TO ARRAY 'FVAL ' TO FACILITATE ANALYSIS. 

C 

DO 40 L: I ,50 
IERCT=0 
C 

C READ IN OF ONE BLOCK OF DATA FROM MAG TAPE. 

C FOLLOWING ’MTREAD * THERE IS AN ERROR DETECTION SECTION 
C 

500 CALL MTREAD(TDMP, 100, IERR) 

I F (I ERR .EQ .2) IEOF=l 
IF(IEOF.EQ.l) GO TO 140 
1 F (TERR . E6 .0 ) GO TO 530 
TYPE 21, IERR 
IERCI =IERCT+1 
IFdERCl .LE.5) GO TO 500 
TYPE 22 
GO TO 200 
C 

C ADDITION OF SCALING FACTORS TO EACH ELEMENT OF THE 
C BLOCK OF DATA READ IN FROM MAG TAPE 4 THEN EACH OF 

C THESE ELEMENTS BECOME ONE ELEMENT OF ARRAY ’FVAL*. 

C FOLLOWING THIS EACH ELEMENT IS TESTED TO INSURE IT DOES 
C NOT EXCEED DISPLAY SCOPE RANGE. IF A VALUE 
C EXCEEDS THE RANGE IT IS TRUNCATED TO FIT THE DISPLAY 

C 4 A MESSAGE IS TYPED PUIOR TO ACTUAL DISPLAY OF DATA. 

C 

530 DO 40 1=0,100,5 
IX=I 

I F < I .EQ .0) IX =1 

FVAL(INDX)=TDMP(IX)+150.0+200.0 
I NDX =1 NDX+ 1 

IF(FVAL(INDX-l).LE. 1023.0) GO TO 40 
FVAL (I NDX -1 )= 1020.0 
IDISER =1 DISER+1 
AS CONTINUE 

C 

C I NI TI LIZATI ON AND SET UP OF DATA DISPLAY ARRAY ’ADVAL*. 
C 

140 CALL I NTAB (ADVAL ,3000) 

XI =0.0 
Yl =FVAL( 1) 
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CALL POINT (ADVAL, XI ,YI> 

DO 1003 L - 2 , 1000 
XEsXl + 1 . 0 
YErFVAL(L) 

1F(YE-YI .LT. 0.999) CALL POl NKADVAL ,XI ,Yl > 

CALL LIN£(ADVAL,X1,YI,XE,YE> 

XI =XE 
1000 Y1=YE 
C 

C IF ANY DISPLAY VALUES HAVE BEEN TRUNCATED TO FIT THE DISPLAY 
C TYPE THE FOLLOWING ERROR MESSAGE. 

C 

IFUDISER .NE.0) TYPE 29.1DISER ' 

29 FORMAT C* <--- ERROR: VALUES TO BE DISPLAYED ARE OUT OF 

1 THE RANGE OF DISPLAY —»*,/,* \7X,I4,* VALUES HAVE BEEN 

2 TRUNCATED TO FIT THE DISPLAY*,//) 

C 

C CALL TO SUBROUTINE ’SLOPCA * TO ANALIZE FLIGHT DATA 
C 

CALL SLOPCA 

IF(IEOF.EQ.l) GO TO 200 
GO TO 26 
C 

C OPTION TO CONTINUE PROCESSING IF A MINOR ERROR 

C OCCURS OR EXIT IF A MAJOR ERROR OCCURS. 

C 

200 IF (IERCT .GT .5) PAUSE 'WHAT DO YOU WANT TO DO?* 

C 

C TRANSFER OF CONTROL TO STATEMENT 2000 IF THE 
C PROGRAM IS OPERATING OFF THE DISK} OTHERWISE, 

C THE PRCGP.ftt REQUESTS DATA STORAGE INFORMATION. 

C 

IF(IDSK.EQ.l) GO TO 2000 
TYPE 20 A 

204 FORMAT (' TYPE THE OF THE DTA YOU WANT THE DATA TO BE 
1 WRITTEN. \S) 

ACCEPT 25, IOTA 
I DTA =1 DTA+8 

2000 I F (IDSK .EQ . 1 ) IOTA =IDSK 
C 

c INITIALIZATION OF OU1PUT STORAGE DEVICE TO ACCEPT 
C POSITIVE SLOPES FOR CONDUCTIVITY CALCULATIONS 
C 

CALL OFILE (IDTA , *POS *) 

Jrl 

N=1ND-1 

C ~ ' “ 

C DO LOOP TO WRITE FIRST POSITIVE SLOPE DATA 
C THEN NEGATIVE SLOPE DATA ONTO THE STORAGE CEVICE 
C 
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21 5 DO 220 1=1, N 

IF(J.EG.l) WRITE(IDTA,230) SLOPEU.J), IIF(I,J) 

220 1F(J ,EQ .2) WRITEUDTA ,230) SLOPEU.J), T1F(I,J), TIF(I,3) 
230 FORMAT (FI 0.5) 

CALL RELEAS (IOTA ) 

JsJ + 1 

IFCJ.GT.2) GO TO 270 
C 

C INITIALIZATION OF OUTPUT STORAGE DEVICE TO 
C ACCEPT NEGATIVE SLOPES FOR CONDUCTIVITY 
C CALCULATIONS. 

C 

CALL OFILEUDTA, ’NEG*) 

N=lNOF-l 
GO TO 215 
270 N=I ND- l 

IF(IDSK.NE.t) IDTA =IDTA-8 
J=1 

IF(IDSK.EQ.l) GO TO 2010 
C 

C PRINT OUT OF HEADING A FILE INFORMATION 
C DATA FOR POSITIVE SLOPES 
C 

PRINT 206 

206 FORMAT ( ' 1 THE FLIGHT DATA HAS BEEN WRITTEN ON DTA'.IX, 

I* UNDER THE NAMES POS.DAT (FOR POS . COND.) A 

2 NEG .DAT (FOR NEG. COND.)*) 

PRINT 207, IDTA 

207 FORMAT ( *+ ’ , 40X ,11) 

PRINT 208, FI LID (2) ,FILID(3),FILID(6),FILIQ(7), 
1FIL1D(8),FIL1D(9) 

208 FORMA T( ' - ’ ,20X , 'DATE OF LAUNCH: *,2A5,5X, TYPE OF JAY:' 

1 ,2A 5, 5X, ’COMMENTS: *,2A5) 

C 

C PRINT THE VALUES OF THE FLIGHT SLOPES 
C 

2010 PRINT 2020 

2020 FORMAT ( * THE FLIGHT DATA HAS BEEN WRITTEN ON THE DSK 

1 UNDER THE NAMES POS.DAT (FOR POS. COND.) A NEG. DAT (FOR 

2 NEG. COND.) ') 

PRINT 208,FILID(2),FILID(3),FILID(6),FILID(7), 

1 F1LID(8),FILID(9) 

2030 PRINT 210 

210 FORMAT ( THE VALUES OF THE SLOPES OF THE FLIGHT 

1 CURVES FOLLOWS : ',/) 

280 1 F(J .EQ .2 ) PRINT 245 

I F(J ,GE .2 ) GO TO 249 

C 

C PRINT OUT OF HEADING 4 FILE INFORMATION 
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C DATA FOR NEGATIVE SLOPES. 

C 

PRINT 240 

240 FORMAT ( * FLIGHT SLOPES USED TO CALCULATE POSITIVE 
1 CONDUCTIVITIES FOLLOWS ’,//,’+ \3X , ‘SLOPE ’, 13X, 
2’TIME IN FLIGHT ’,//,*■*■’, 25X, ’IN SEC*,//) 

245 FORMAT! *- FLIGHT SLOPES USED TO CALCULATE NEGATIVE 
l CONDUCTIVITIES FOLLOW: *,//,•+’ ,3X , ’SLOPE ’, 13X, 
2’TIME IN FLIGHT ’+\25X, ’IN SEC’,//) 

C 

C ACTUAL PRINT OUT OF SLOPE EITHER POSITIVE 
C OR NEGATIVE 
C 

249 DO 250 I=1,N 

250 PRINT 260, SLOPE(I,J>, TIFCI.J) 

260 FORMAT!’ \F 10.5, 12X ,F 10.5) 

J=J + l 

IFU.GT.2) GO TO 290 
N:INDF-l 
GO TO 280 
290 ISKIP : 1 

IF(IDSK.NE.l) IDTA=I3TA+8 
CALL DISPLY (0 ) 

CALL R£LEAS(3> 

C 

C CHAIN BACK TO LINK ’MAIN*. 

C 

CALL CHAIN!0,IDSK, ’MAIN’) 

END 


C ************************************************************ 

c 

C SUBROUTINE SLOPCA s CALCULATION OF SLOPES OF 

C FLIGHT DATA. 

C 

C ************************************************************* 
SUBROUTINE SLOPCA 

COMMON LP, ISHOW, XMAX, XMIN, YMAX , YM1N , INTENS, 1SCALE 
COMMON /A/ ADVAL (3000) , GRIDII000) 

4, DEC ! 1 00) /D/ STAR (20) , SEC C 1 3) /J / ASK (4 ) 

1 /£/ SLOPE! 150,2) /F/ X!4), Y!4) /H/ TIFU50,3> 

2 /I/ TIME(250) /M/ IND, IFLAG, IEOF, KIN, I NDF , NFLG 

3 /P/ IL1NEU00), RB ( 5) /S/ FVAL! 1005) /DC/ FILID(i0), 

4 TDMP(ia0) 

C 

C SETUP LITERAL SHOWN FOR DISPLAY ON SCOPE, TO ALLOW 
C USER TO CHOOSE THE » OF CALIBRATE SLOPES TO BE 
C ANALIZED IN ONE GIVEN DISPLAY 
C 
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ITS =0 
IEXPN =0 

IF(NFLG.GT.l) GO TO 5 
NFLG=NFLG+l 

CALL CHRG£N(GRID, '# OF SLOPES 12 3 4', 

1 212.0, 62.0, 2, 0) 

C 

C DO LOOP TO INSERT VERTICAL SEPARATION LINES BETWEEN 
C EACH # ABOVE, TO DEFINE SELECTION AREA FOR EACH #. 

C 

DO 100 1=1,5 

CALL POINTCGRID, RBU), 60.0) 

100 CALL LINECGRID, RBU), 60.0, RBU), 80.0) 

C 

C INITIALIZATION OF ARRY STAR 4 SEC. THEN CONSTRUCT 
C MOVEABLE STAR. 

C 

5 CALL INTAB (STAR ,20) 

CALL INTAB (SEC , 12) 

ISHOW=0 

CALL POINTtSEC, 155.0,75.0) 

CALL LINE (STAR, 20. 0,20. 0,0. 0,20.0) 

ISHOWrl 

CALL LI NECSTAR, 0.0, 20. 0,40. 0,20.0) 

ISKOW=0 

CALL LI NE(STAR, 40. 0,20. 0,40. 0,0.0) 

ISHOWrl 

CALL LINE (STAR, 40. 0,0. 0,0. 0,40.0) 

ISHOW=0 

CALL LINE (STAR, 0.0, 40. 0,20. 0,40.0) 

ISHO'Jrl 

CALL LINE (STAR, 20. 0,40. 0,20. 0,0.0) 

ISHOW=0 

CALL LINE(STAR, 20. 0,0. 0,0. 0,0.0) 

ISHOWzl 

CALL LINE (STAR, 0.0, 0.0, 40. 0,40.0) 

C 

C MPOINT MOBILIZES THE STAR, THEN ALL ARRAYS ARE 
C DISPLAYED, THEN CALL L1TEPN TO ACTIVATE THE LITE PEN. 

C 

CALL MPOI NT (SEC , 155.0,75.0) 

CALL DISPLY ( 5,SEC ,S TAR ,GRI D ,ADVAL .TIME ) 

10 CALL LI TEPN (NSLOPE , XS , YS ) 

IF (NSLOPE.EQ .0) GO TO 10 
C 

C DECISION SECTION! 

C REJECT 7777 

C 

IF((XS .GT • 680.0 ) .AND . (XS .LT . 800.0) .AND . (YS ,LE. 40.0) ) 
1 GO TO 50 
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c 

C END???? 

C 

I F ( CXS .GT. 920.0). A NO. (XS .LT . 1 323 .0) .AND . (YS ,L£ . 40. 0 ) ) 
1 GO TO 60 
C 

C WAT IS THE t OF SLOPES TO BE 
C CALCULATED 7777 
C 

IF< (XS.GT. 380.0). AND. (XS .LI. 433.0). AND. (YS.GT. 60.0). 
IAND.(YS.LT. 80.0)7 NUKSLsl 

I F( (XS.GT. 433.0). AND ,<XS .LT .500 .0) .AND . (YS .GT . 60.0) . 
IAND . (y3 .LT .80.0)) NUP!SL=2 

I F < (XS.GT. 500.0). AND. (XS .LT .548.0) .AND . (YS .GT .60.0) . 
IAND , (YS .LT . 80 .0 ) ) NUMSL=3 

IF( (XS.GT. 548. 3).AND.(XS .LT .600. 0) .AND . (YS .GT . 60.0) . 
IAND. (Y5.LT. 80.0)) NUMSL=4 

IF( (XS . GT. 665.0). OR .(XS .LT .380 .0 ) .OR . (YS .LT . 60.0) . 

I OR . (YS .GT . 80 . 0 ) ) GO TO 10 
C 

C CALCULATION LOOP FOR THE # OF SLOPES TO BE CALCULATED 
C 

DO 40 I=1,NUMSL 
ISA V: 1 
NFST =0 
C 

C ACTIVATE LITEPEN TO ACCEPT FIRST POINT IN CALCULATING 
C SLOPES. 

C 

20 CALL LITEPNdACC ,X2,Y2) 

IF(IACC ,EQ ,0) GO TO 20 
CALL MPOINT(SEC,X2,Y2) 

IF(ISAV.NE.l) GO TO 21 
C 

C SLOPE FOR POSITIVE CONDUCTIVITY??? 

C 

IF((Y2.LT .40,0) ,0R.(Y2.GT .80.0). OR. (X2.LT .650.0) 

1. OR. (X2.GT. 675.0)) GO TO 3000 
ITYSL=l 
GO TO 3010 
C 

C SLOPE FOR NEGATIVE SLOPE??? 

C 

C IF LITE PEN DETECTION DOES NOT FALL 
C WITHIN ABOVE TWO SECTIONS THEN 'GO TO 20' 

C AND WAIT FOR ONE OF THE TWO RESPONSES. 

C 

3000 IF((Y2.LT.40.0).OR .(Y2.GT.80.0).OR .(X2.LT.D78.0) 

1. OR. (X2.GT. 1003.0)) GO TO 20 
ITYSL=2 
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5010 CALL WAITU000) 

CALL LITEPN(ICLE,X4,Y4) 

C 

C MOVE STAR TO READY POSITION 1 4 TEST FLAGS 
C ’1 lYSL * AND ’ITS'. 

C 

15 CALL MPOINKSEC, 155.0,75.0) 

IFUTYSL.EQ .l.AND .ITS .EQ .1) GO TO 2000 
C 

C CALL TO 'Ll TEPN * TO ALLOW MOVING OF STAR 
C ANYWHERE ABOVE 200 RASTER UNITS ON THE 
C DISPLAY SCREEN, CHECK OTHER FLAGS, MOVE 
C STAR TO DETECTED POINT, 4 PLACE X 4 Y 

C COORDINATES DETECTED INTO 'X(ISAV)' 4 'Y(ISAV) 

C AT NO TIME WILL '1SAV ' BE ANY OTHER VALUE THAN 

C ONE. 

C 

21 CALL LITEPN(IACC,X2,Y2) 

IFCIACC .EQ .0) GO TO 21 

IF((Y2.LE.200.0).OR.(ISAV.GT.1).OR.(NFST.EQ.D) GO TO 200 
CALL MPOI NT (SEC ,X2,Y2) 

X (ISA V ) =X2 
Y (ISAV) =Y2 
YLOCrY 2+400.0 
GO TO 21 

200 IF(ITS.EQ.l) GO TO 23 

IF (ISA V.GT . 1 ) GO TO 23 
IF(IEXPN.EQ.l) GO TO 22 
C 

C TESTS TO SEE IF OPERATOR DESIRES TD MAGNIFY 
C DATA BY A FACTOR OF 4. IF NOT ’GO IU 22’. 

C 

IF((Y2.LT.55,0).OR.(Y2.G£.76.0).OR. (X2.LE .710.0) 

1. OR. (X2.GT. 740.0)) GO TO 22 
CALL DISPLY (0) 

I EXP N = 1 
ICHKS : 1 
YL0C2=YL0C 
C 

C PERFORM MULTIPLICATION ON FLIGHT DATA STORED 
C IN ARRAY ’FVAL’, THE INITIALIZE DISPLAY ARRAY 
C 'ADVAL* 4 TRANSFER DATA FROM 'FVAL' TO 
C 'ADVAL' 

C 

DO 1000 L=I , 1001 

iB00 FvaL <L ) = (.!• VALIL) -650. 10 J+4.0+YLOC2 
CALL INTAB (ADVAL, 3000) 

XI =0.0 
Y I =FVAL(l) 

I F ( Y I .LT. 150.0) YI =150.0 
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E 010 
22 


I F< YI .GT. 1020.0) YI =1020.0 
CALL P0INI(ADVAL,XI,Y1) 

00 1010 L=2, 1000 
XE=XI+1 .0 
YE=FVAL(L) 

IF(YE.LT. 150.0) YE:150.0 
IF(YE.GT. 1020.0) YE : 1020 . 0 
IFCYE-YI .LT.3.999) CALL POI NT (ADVAL ,XI ,YI > 
CALL LINE (ADVAL, XI ,Y1 ,XE,YE) 

XI =XE 
Y1=YE 
GO TO 24 

IF(IEXPN.NE.l) G TO 23 


TEST TO SEE IF MAGNIFIED DISPLAY SHOULD RETURN TO NORMAL. 

IF( (Y 2. L T. 55, 0), OR. (Y 2. GE. 76.0). 0R.CX2.LT. 755.0) 

1. OR. (X2.GT. 845.0)) GO TO 23 
CALL DISPLY <0) 

X2:155.0 
Y2=75.0 
ICHKS =0 

REDUCES EXPANDED VALUES TO NORMAL SIZE 
THEN INITIALIZE DATA DISPLAY ’ADVAL' 4 
TRANSFER DATA FROM ARRAY 'FVAL ’ TO ARRAY 
’ADVAL '. 

DO 1020 L:l,1001 

020 FVAL (L)s (FVAL CL )-YL0C2)/4. 0+650.0 

CALL INTAE(ADVAL,30C0) 

XI =0.0 
YI=FVAL(1) 

CALL POI NT (ADVAL ,XI ,YI ) 

DO 1030 L =2, 1 000 
XE=XI+1 .0 
YE=FVAL(L) 

IF(YE-YI .LT. 0.999) CALL POI NT (ADVAL ,XI ,YI ) 

CALL LINE (ADVAL, XI ,YI ,XE,YE) 

XI =XE 
030 YI =YE 

IEXPN =0 
GO TO 24 

TEST FLAGS TO SEE WHICH ARRAYS ARE TO 
BE DISPLAYED. 


23 CALL DISPLY(5,SEC, STAR, GRID, ADVAL, TIME) 

I F( (I TYSL .D8 .2 ) .AND . (I TS .EQ .0 ) .AND . (I EXPN .EQ . 0) ) 
1 CALL DISPLY (6, SEC ,S TAR ,GRID ,ADVAL 
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2,TIME ,DEC ) 

1F((ITYSL.£Q.2.AND.ITS.EQ.1).AND.(IEXPN.EQ.0>) 

I CALL DISPLY (5, 3 EC ,S TAR , GRI D ,ADVAL 
2, DEC) 

C 

C TEST TO SEE IF ’ACCEPT’ GAS BEEN DETECTED| 

C IF SO, 'GO TO 50*. OTHERWISE TEST TO SEE 
C IF ’ITYSL' EQUALS TWO A A TWO SECOND 
C DISPLAY C2 SEC DIS’) HAS SEEN DETECTED. 

C IF NEIlhER OF THE ABOVE TWO OCCUR, TEST 

C TO SEE IF ’NORMAL ’ HAS BEEN DETECTED TO 

C RETURN FROM A TWO SECOND DISPLAY. 

C 

IF((X2. GT. 380.0). AND. (X2.LT. 630.0). AND. (Y2.LE. 40. 0)> GO TO 30 
IF((ITYSL.EG.2).AND.(X2.GT.15.0).AND.(X2.LT.H0.0).AND. 
HY2.GT.15.0) .AND .(Y2.LT .30.0)) GO TO 2000 
1F((1TYSL.EQ.2).AND.(X2.GT .755.0 ). AND. (X2.LT. 845.0). AND. 
I(Y2.GT.55.0).AND.(Y2.LT.76.0).AND.(ITS.EQ.m GO TO 2000 

MOVE STAR TO DETECTED POINT A TEST TO 
SEE WHICH ARRAYS SHOULD BE DISPLAYED. 

CALL MPOI NT (SEC , X2, Y2) 

CALL DISPLY (5, SEC ,S TAR, GRID.ADVAL, TIME) 

IF( (ITYSL. EG .2). AND. (ITS .EG .0) .AND . (IEXPN .EQ .0) ) 

I CALL DISPLY(6, SEC, STAR, GRID.ADVAL 
2, TIME .DEC) 

IF ( (ITYSL. EQ.2).AND. (ITS .EQ . 1 ) .AND . (IEXPN .EQ .0) > 
l CALL DISPLY(5,SEC, STAR, GRID.ADVAL 
2, DEC) 

CALL ! KTAB (ILI NE ,133) 

X 4 Y COORDINATES OF DESIRED POINT ON CURVE 

X (ISA V ) =X2 
Y(ISAV);Y2 

CALL POI NT (ILI NE , X(l>, Yd)) 

IF(ICHKS.EQ.l) NFST = 1 
IFdCHKS .EG .0) NFST =0 
IFUSAV.LT.2) GO TO 21 
C 

C CHECK TO SEE IF THERE WAS AN ERROR IN SELECTING 
C THE SECOND POINT ON THE CURVE. IF THERE WAS 
C DISPLAY THE ERROR MESSAGE SHOWN, CONNECT THE 

C TWO POINTS SELECTED WITH A LINE, AND DISPLAY 

C ALL ARRAYS. IF NO ERROR JUST CONNECT THE 

C POINTS WITH A LINE A DISPLAY ALL ARRAYS. 

C 

CALL LINEdLINE, Xd), Yd), X(2), Y(2)> 
IF((Yd).LT.Y(2)).AND.(Xd).LT.X(2))) GO TO 25 
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CALL CHRGENULINE, ERROR REG. SLOPE — > ',210. 3, 82. id, 2,0) 
CALL DISPLY ( ) 

DO 26 IB=l,4 

CALL DISPLY ( 1 ,IL1NE) 

CALL UAIK253) 

26 CALL DISPLY (0) 

C 

C TEST TO SEE WHICH ARRAYS SHOULD BE 
C DISPLAYED. 

C 

25 CALL DISPLY (6, SEC ,S TAR-,GRI D ,A OVAL ,TIME ,1 LI NE> 

IF( (I TYSL .00 .2) .AND . (ITS .EQ .0 > .AND . (IEXPN .EQ. 0) > 

1 CALL DISPLY(7, SEC, STAR, GRID, 

2A DUAL, TIME, ILINE .DEC) 

IFOITYSL.EQ .2). AND .UTS -EC.)). AND. (IEXPN. £6.0)) 

1 CALL DISPLY (6^SH:,STAU,GR1D, 

2A DUAL ,1 LI RE, DEC) 

GO TO 2i 
C 

C WAIT 15180 MSEC, CLEAR LITE PEN VALUES, & MOVE STAR 

C TO READY POSITION. INCREASE ISAV BY 1, A IF ISAV« 

C OR =2 GO TO 50, TO ACCEPT SECOND COORDINATES OF DATA 

C IF ISAV>2 CALCULATE SLOPES. 

C 

50 CALL WAIT ( 1502) 

CALL LITEPNUCLE, X4, Y4) 

CALL MPOI NT (SEC , 155.0, 75.0) 

ISAV=ISAV+1 

I F ( ISA V .LE .2) GO TO 21 
IFUTYSL.EG.2) GO TO 55 
TIFUND, 1 >=<KIN-1)+X( 1)40.01 
C 

C SLOPE CALCULATIONS TO COMPENSATE FOR 
C FACTOR OF FOUR EXPANSION IN 1 POSITIVE 
C SLOPE CALCULATIONS 
C 

SL0PE(IND,l):(YC2)-Y(l))/(X(2)-X(l)) 

IFUEXPN.EQ.l) SLOPECIND , l )=SLOPE(IND, 1 )/4.0 
GO TO 39 

55 IF(ITS.EG.l) GO TO 56 

TIF(INDF,2):(KIN-t )+X (1 >40.01 

TI F (I NDF, 5 )=<((KIN-1)+X (2)40.0 J >-TIF (I NDF,2>>/2.0 
SL0PE(INDF,2)=(Y(2)-Y (l))/(X(2)-X (1)) 

C 

C TEST TO COMPENSATE FOR FACTOR OF FOUR 
C EXPANSION IN NEGATIVE SLOPE CALCULATIONS 
C 

IFOEXPN.EQ ,1 ) SLOPE(INDF,2)=SLOPE(INDF,2)/4.0 
GO TO 57 
C 
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C STATEMENTS TO COMPENSATE FOR TWO SECONS 
C TIME EXPANSION IN NEGATIVE SLOPE CALCULATIONS. 

C 

36 TIF(INDF,2) =1 TMVA+X ( l )*0 .002 

TlFUNDF ( 3> = (<llMVA+X(2)*0.002)-TIF(INDF,2)>/2.0 
SL0PE(INDF,2):((Yt2)-Y C l > ) /(X (2) -X (I ) ) >*5.0 

37 INDF=1NDF+1 
GO TO 40 

39 IND: IND-H 

40 CONTINUE 

50 IFLAG=0 

GO TO 70 
60 1 EOF : 1 

70 CALL OISPLY (0) 

GO TO 80 
C 

C SECTION TO CHANGE TIME SCALE FACTOR FOR 
C A TWO SECOND DISPLAY. 'ITS* EQUAL ONE IMPLIES 
C THAT THESE ALREADY HAS BEEN A TIME SCALE 
C EXPANSION 4 DISPLAY IS TO RETURN TO NORMAL. 

C 

2000 CALL DISPLY<0> 

IFCITS.EG.l) GO TO 2130 
I TS =1 
ISA V = 1 

2005 ITMVA=XIN-1 
ST=0.0 
EN=1 00.0 
C 

C DO LOOP TO DETERMINE ABOUT WHICH 
C TIME Ea? ANSI ON IS TO OCCUR. 

C 



DO 2018 IN=1,10 
I F (X ( 1) .GE .S I .A ND .X <1 > .lE .EN ) 

GO TO 2030 

2010 

ST=ST+100.0 

EN=EN+100.0 


2020 

TYPE 2020 

FORMAT ( ' ***ERR OR REGION NOT 

FOUND*** \/, * \10X,* 

2030 

ITRYING AGAIN *,//) 
GO TO 2005 
ITMVA =1 TMVA+IN 



c 

C STATEMENTS TO DETERMINE HOW MANY RECORDS 
C THE MAG TAPE SHOULD BACKSPACE • 

C BACKSPACE THE MAG TAPE, AND SET UP 

C NEW EXPANDED TIME SCALE FOR DESPLAY 

C IN ARRAY ’DEC’. 

C 

LRUiST 

IBKSP : ( 1 000-LRU ) /20+ 1 
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IRETsIBKSP-ta 

CALL MTSKIPO.-IBKSP) 

XCOOR =460.0 
DO 2040 I N=0, 1 
ITMVA =1 TMVA+I N 
ENCODE (4 ,2035, ASK ) ,ITP1VA 
2035 FORMAT (I 4) 

ASK (2>=0.0 

CALL CHRG£NOEC,ASK .XCOOR, 130.0,2,0) 

2040 XCOOR =XCOOR+500.0 

ITMVA =1 TMVA-2 
LNTH =10 
I NS T = 1 
C 

C INITIALIZE DATA DISPLAY ARRAY 'ADVAL' AND 
C READ IN THE REQUIRED NUMBER ( ’LNTH ') OF 
C RECORDS TOR DISPLAY. 

C 

2045 CALL I N'TAB (ADVAL , 3000 ) 

I NDX = 1 

DO 2090 NR =1, LNTH 
IDISER =0 
IERCT = 0 

2050 CALL MTREAD (TDMP , 100,1 ERR ) 

C 

C TEST FOR MAG TAPE INPUT ERROR AND IF AN 

C ERROR OCCURS TAKE THE APPROPRIATE 

C BRANCH. 

C 

IF (IERR .£Q .2) IEOF=l 
IF(IEOF.EQ.l) GO TO 2110 
1FUERR .EQ .0) GO TO 2080 
C 

c type out error messages if an error 

C OCCURS . 

c 

TYPE 2060, IERR 

2060 FORMAT ( ' ***MTIO ERROR*, 1 1,' PROCESSING CONTINUED***',//) 
IERCT =IERCT+ 1 
I F (IERCT .LE . 5) GO TO 2050 
TYPE 2070 

2070 FORMAT ( ’ ',8X, 'ERROR WILL NOT CORRECT PROCESSING ENDED',///) 
TECF = 1 
GO TO 70 
C 

C DO LOOP TO PERFORM PROPER^ OPERATIONS ON. 

C DATA-READ IN FROM MAG' TAPE AND PUT THE RESULTS 
C ONTO DATA ARRTY *FVAL ' • THEN TEST 
C DATA TO SEE IF ANY POINTS EXCEED DISPLAY 
C SCOPE RANGE; IF SO, THEN TRUNCATE THE VALUE 
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C TO FIT THE DISPLAY 4 PRINT OUT ERROR MESSAGES. 

C 

2080 DO 2090 LI =0, 100, INST 

IF(LI .EQ.0. AND. INST. ES.l) GO TO 2090 
LB:LI 

IF (LI .EG .0) LB=t 

FVAL<INDX)= TQMPCLB)* 1 50.0-4-200.0 
1NDX=1NDX+1 

I F(FVAL (I NDX-1 > .LE . 1023 .0 ) GO TO 2090 
FVAL(INDX-I): 1020.0 
IDISER sIDISER + 1 
2090 CONTINUE 

1FODISER .EG .0) GO TO 2110 
TYPE 2100, IDISER 

2100, FORMAT ( ' ----ERROR: VALUES TO BE DISPLAYED ARE OUT OF 

1 THE RANGE OF DESPLAY — * ', 7X.I 4, 'VALUES 

2 HAVE BEEN TRUNCATED TO FIT THE DISPLAY.',//////) 

C 

C TRANSFER OF DATA FROM ARRAY *FVAL ' TO 
C DISPLAY ARRAY ’ADVAL'. 

C 

2ii0 xi=a.0 

Yl =FVAL( 1 ) 

CALL POI NT (ADVAL ,XI ,YI> 

DO 2120 NI :2, 1000 
X£=XI + 1 .0 
YE=FVAL(NI) 

IFCYE-YI .LT. 0.999) CALL POI NT CADVAL ,XI , YI ) 

CALL LINE(ADVAL,XI ,YI ,XE,YE) 

XI =XE 
2120 YI=YE 

C 

C DISPLAY OF APPROPRIATE ARRAYS 4 RESET 
C ARRAY DEC WHEN NORMAL IS DETECTED. 

C 

CALL DISPLY C5.SLC ,S TAR ,GRI D .DEC .ADVAL) 

IF(ITS.EG.0) CALL DISPLY C 5 , SEC .STAR, GRID.ADVAL, TIME) 
GO TO 15 
2130 1 TS :0 

CALL INTABIDEC, 100) 

CALL CHRGEN (DEC . '2 SEC DIS ',20.0, 15.0,2,0) 

C 

C DETERMINE HOW MANY RECORDS MUST BE 
C BACKSPACED 70 BRING BACK ORIGINAL TEN 
C SECOND DISPLAY. 

C 

IBKSP :LRU/20+9 
CALL MTSK1P(0, -IBKSP) 

LNTH =50 
INSTrS 
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GO TO 2045 
C 

C IF FINAL SLOPE IN DISPLAY HAS 
C BEEN TAKEN A THERE IS A TWO SECOND 
C DISPLAY IN EFFECT, BACKSPACE THE MAG TAPE 
C 'IRET ’ RECORDS TO BEGIN1NG OF ORIGINAL TEN 
C SECOND DISPLAY. THIS IS NECESSARY TO INSURE 
C NO LOSS IN DATA. THEN RETURN TO CALLING 
C PROGRAM. 

C 

80 IF(ITS.EQ.I) CALL MTSKIP (0.1RET) 

RETURN 

END 
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C *************************************************************** 
C CONCs determines boih the positive 

C AND NEGATIVE CONDUCTIVITIES 

c ************************************************************* 

COMMON /C/ AVCASL /K / SI GPOS < 1 50) , SIGNEG(iS3) 

I /L/ TIMEIF<153,2) /DA/ IDSK 
C 

C STATEMENT REQUIRING OPERATOR ACTION UPON 
C ENTERING SUBROUTINE. 

C 

PAUSE 'CALCULATION OF CONDUCTIVITIES' 

C 

C IF THE PROGRAM IS OPERATED OFF OF THE DISK 

C SYSTEM CONTROL IS TRANSFERRED TO STATEMENT 

C NUMBER 233; OTHERWISE, TYPE CUT THE FOLLOWING 
C QUESTION PERTAINING TO STORAGE OF DATA ON 
C DEC TAPE. 

C 

IFCIDSK.EQ.i) GO TO 203 
6 TYPE 1 

1 FORMAT ( ' TYPE # OF DTA ON WHICH CAL SLOPES ARE 
I STORED '.I) 

ACCEPT 2, IANS 

2 FORMAT <1 2) 

C 

C INITIALIZE DEC TAPE INDICATED FOR READ IN OF 
C THE AVERAGE OF THE CALIBRATION SLOPES. 

C 

IFCIANS.GT.7) GO TO 6 
I DTA =IANS+8 
GO TO 210 
200 IDTA =1 DSK 

C 

C INITIALIZATION OF STORAGE DEVICE FOR INPUT OF 
C AVERAGE CALIBRATION SLOPE AND INPUT CF SAME. 

C 

210 CALL IFILECIDTA, 'AVCAS ') 

READCIDTA, 3) AVCASL 

3 FORMA T (F l 3, 5) 

CALL RELEAS(IJTA) 

C 

C INITIALIZE N, J, & CALRES (CALIBRATION RESISTANCE), 

C DETERMINE SMALL R 4 LARGE R (PROEE DIMESIONS), 

C 4 THE RATIO OF LARGE R TO SMALL R MULTIPLIED BY .5 
C 

N = 1 
J s 1 
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CALRESs 0.2E+12 

211 TYPE 212 

212 FORMAT ( * 00 YOU DESIRE TO CHANGE THE PROBE 

1 RADII? Y OR N \S> 

ACCEPT 213, ANS 

213 F0RMAT(A3> 

IF(ANS .£Q . ’N ') GO TO 214 
1 F (ANS . NE . 'Y ' ) GO TO 211 
TYPE 215 

215 FORMA T ( * TYPE THE COLLECTOR RADIUS IN CM, 

1 EX. 2.5 *,S> ’ 

ACCEPT 216, SMALR 

216 FORMAT (F 13.5) 

TYPE 217 

217 FORMAT ( ' TYPE THE GUARD RING RADIUS IN CM, 

1 EX. 4.6 ’,$) 

ACCEPT 216, LARGR 
GO TO 221 

214 SMALR=(17. 0/32. 01*2. 54 
LARGR=(23. 3/16. 01*2. 54 

221 SMALR =SMALR**2 

RATIO=LARGR /SMALR 
C 

C CALCULATE THE FACTOR TO BE MULTIPLIED BY THE 
C FLIGHT SLOPE. 

C 

CONST = (RA TI 0/(2*CALR£S 1 l/AVCASL 
IF(IDSK.EQ.l) GO TO 220 
5 TYPE 10 

10 FORMAT ( ' TYPE THE » OF THE DTA ON WHICH THE FLIGHT DATA HAS 

l SEEN RECORDED. \ S) 

ACCEPT 20, IOTA 
20 FOR MAT (I) 

IDTA =IDTA+B 
GO TO 230 
220 IDTA =ID5K 

C 

C INITIALIZE STORAGE DEVICE FOR INPUT OF DATA. 

C 

230 CALL IFILEdDTA, ’POS ') 

30 READ (I DTA , 43, END=50) SLOPE, TIFL 

40 FORMAT (F 13. 5) 

C 

C CALCULATE POSITIVE CONDUCTIVITY, INCREMENT N, 

C THEN GO TO 30 TO CONTINUE THE PROCESSING 
C OF DA TA . - 
C 

SIGPOS (N ) =CONST*SLOPE 

T1MEIF(N,1)=TIFL 

N=N+I 
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GO TO 33 

50 CALL RELEAS (IOTA ) 

C 

C STORE POSITIVE CONDUCTIVITIES ON THE STORAGE DEVICE UNDER 
C THE NAME PS1.DAT 4 THEN PRINT OUT THE POSITIVE 
C CONDUCTIVITY. 

C 

CALL OFILEdDTA, 'PSI *> 

DO 63 1=1, N-l 

60 WR ITE ( I DTA ,65) SIGPOS(I), TIflEIF (1,1) 

65 FORMAT(£10.5,F10.3) 

CALL RELEAS (1DTA ) 

C 

c if OPERATION IS OFF DISK "GO TO 240" TO PRINT 
C HEADINGS AND POSITIVE CONDUCTIVITIES} OTHERWISE, 

C PRINT THE FOLLOWING HEADING AND THEN THE POSITITE 
C CONDUCTIVITIES. 

C 

IFdDSK.EQ .1 ) GO TO 240 

ISTDTA=IDTA-8 

PRINT 70, ISTDTA 

70 FORMAT Cl THE POSITIVE CONDUCTIVITIES HAVE BEEN STORED 

1 ON DTA*, IX,* UNDER THE NAME PSI .DAT *,/, ’+’,52X,I 1 , ////, 

2* THE POSITIVE CONDUCTIVITIES FOLLOW: ',///,'+• , 

3 ’CONDUVTI VI TY ’ , 4X , ‘TIME IN FLI GHT *,//,*+• ,2X, * (IN MHOS/CM) \ 
48X, '(IN SEC) ',//) 

GO TO 260 

240 PRINT 250 

250 FORMAT ( *1 THE POSITIVE CONDUCTIVITIES HAVE BEEN STORED 
l ON THE DSK UNDER THE NAME PSI .DAT 10X , 

2' THE POSITIVE rONDUCTI VI TIES FOLLOW: ',///,'+ ', 

3 'CONDUCTIVITY ',4X, 'TIME IN FLI GHT ',//,'+' ,2X, ' (IN MHOS 
4/CM) \8X, '(IN SEC ) ',//) 

2$0 DO 80 1=1, N-l 

83 PRINT 90, SIGPOS(I), TIMEI F (1 , 1 ) 

93 FORMA T ( ' ' , IX ,E 1 0 . 5 , 7X , F 1 0 . 5 ) 

C 

C INITIALIZE STORAGE DEVICE TO READ IN NEGATIVE 
C SLOPES TO CALCULATE NEG. CONDUCTIVITIES, INCREMENT J, THEN 
C GO TO 100 TO CONTINUE THE PROCESSING OF DATA. 

C 

CALL I FILE (IDTA , 'VOST *) 

READdDTA , 40) VDS 
CALL RELEAS (IDTA ) 

CALL IFILEdDTA, ; NEG ; ) 

100 READdDTA, 40, END = 110) SLOPE, TIFL, TE 
SIGNEG(J)=CONST*SLOPE 
TIMEI F (J , 1 ) =TI FL 

TIMEI F(J,2)=(VDS*TE*2.0)/(3. 1414*LARGR) 

J=J + 1 
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QO TO 100 
C 

C STORE NEG. COND. ON THE STORAGE DEVICE UNDER THE NAME 
C NSJ.DAT 4 THEN PRINT OUT IHE NEG. COND. 

C 

110 CALL RELEAS (IDTA ) 

CALL OFILE (IDTA , ’NSI ’) 

DO 120 1=1, J-l 

120 VRITEdDTA, 125) SIGNEGCI), TIKEIFCI,!), TlMEIF(I t 2) 

125 FORMAT (E13.5,F10.5,F10.5) 

CALL RELEAS (IDTA) 

C 

C IF OPERATION IS OFF THE DISK SYSTEM "GO TO 265", 

C PRINT APPROPRIATE HEADING AND THEN NEGATIVE 
C CONDUCTIVITIES. OTHERWISE, PRINT THE FOLLOWING 
C HEADING AND THEN THE NEGATIVE CONDUCTIVITIES. 

C 

IF(IDSK .EQ .1) GO TO 265 
PRINT 130, ISTDTA 

130 FORMAT ( ' 1 THE NEGATIVE CONDUCTIVITIES HAVE BEEN STORED 
1 ON DTA ' , IX , ’ UNDER THE NAME NSI .DAT ’,/,'-*, 1 0X , 

2* THE NEGATIVE CONDUCTIVITIES FOLLOW: *,///,’+ ’, 

3 ’CONDUCTIVITY \4X, ’TIME IN FLI GHT ’ , AX , *E FIELD’,// 

4, *+ ' , 2X , ' (1 N MHOS/CM) \8X, ’(IN SEC ) 6X, ’ UN V/CM)',//) 

GO TO 280 
265 PRINT 270 

270 FORMAT ( * 1 THE NEGATIVE CONDUCTIVITIES HAVE BEEN STORED 

1 ON THE DSK UNDER THE NAME NSI .DAT ’ 1 0X , 

2’ THE NEGATIVE CONDUCTIVITIES FOLLOW: ’,///, •+', 

3 ’CONDUCTIVITY ’ , 4X , ‘TIME IN FLI GHT ’ , 4X , ’E FIELD’,// 

4, ’+ ’ ,2X , ’(IN MHCS/CM) ’»SX, "(IN SEC ) *,6X, ’ (I N V/CM)’,//) 
280 DO 140 I=l,J-l 

140 PRINT 145, SIGNEGd), TIMEIF(1,1), T1MEIF(I,2) 

145 FORMAT ( ’ ’ , IX ,E 1 0 . 5, 7X ,F 10.5 , 6X ,F1 0.5) 

CALL RELEAS (3) 

C 

C CHAIN TO LINK "MAIN". 

C 

CALL CHAI N( 0,1 DSK , ’MAIN ’) 

END 
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c ******************************************************** 
c 

C MTAH t MATCH TIME A ALTITUDE 

C 

C ******************************************************* 

COMMON /H/ TIF(I50,3> /K / SIGPOS(150), SIGNEG(I50) 

1 /N/ IPS!, INSI /Q / ALT (150,2) /DA/ ID§K 

C 

C STATEMENT REQUIRING OPERATION ACTION UPON 

C ENTERING THIS LINK . 

C 

PAUSE ’MATCH TIMES TO ALTITUDES’ 

C 

C IF OPERATION IS OFF THE DISK SYSTEM 
C "GO TO 5" TO INITIALIZE THE STORAGE DEVICE 
C FOR INPUT OF CONDUCTIVITIES t OTHERWISE, 

C PERFORM TWO ADDITIONS BEFORE INITIALIZATION 
C OF STORAGE DETICE FOR INPUT OF DATA. 

C 

lF(IDSK.EQ.l) GO TO 5 
IPSI =IPSI+8 
INSI =1 NSI + 8 

5 CALL IFILEdPSI , *PSI ’) 

1=1 

C 

C INPUT POSITIVE CONDUCTIVITIES AND STORE THEM 
C IN ARRAY "SIGPOS". 

C 

10 READdPSI ,20, END =30) SIGPOSII), TIFd.l) 

I =1 + 1 
GO TO 10 

20 FORMAT (E10.5,FI0.5) 

30 CALL RELEAS (IPSI ) 

C 

C INITIALIZATION OF STORAGE DEVICE FOR INPUT 
C OF NEGATIVE CONDUCTI VI TI ES , READ THE 
C CONDUCTIVITIES AND STORE THEM IN ARRAY 
C "SIGNEG". 

C 

CALL I FILE (INSI , 'NSI ’) 

J=1 

40 READ (I NSI , 2 l .END = 50 ) SIGNEG(J), TIF(J,2), TIF(J,3> 

21 FORMAT(D10.5,F 10.5.FI0.5) 

J=J+l 

GO TO 43 
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50 CALL RELEAS (I NSI ) 

C 

C TYPE OUT INSTRUCTIONS AND HEADING FOR 

C INPUT OF ALTITUDES. 

C 

TYPE 60 

62 FORMAT ( * AFTER COMPUTER PRINTS OUT CONDUCTIVITY A 

l TIME, TYPE IN C ORRESP ONDI NG * , /, * ALTITUDE 

2. EXAMPLE. ..65. 5 MEANS 65 . 5KM .*,///, ’ CONDUCTI VI TY ' 
3 , 4X , 'TIME (IN SEC) *,4X, 'ALTITUDE:????*,//) 

C 

C TYPE POSITIVE CONDUCTIVITIES ALONG WITH 
C CORRESPONDING TIMES AND WAIT FOR INPUT 
C OF ALTITUDES. 

C 

DO 83 L = l ,1 - 1 

TYPE 70, SIGPOS(L), TIF(L,l> 

70 FORMAT ( * ',4X,E10.5,7X,F10.5,6X,S) 

ACCEPT 90, HEI 
80 ALT (L , 1 ) :HEI 

90 FORMAT (F 10. 5 ) 

C 

C TYPE NEGATIVE CONDUCTIVITIES ALONG WITH 

C CORRESPONDING TIMES AND WAIT FOR 

C INPUT OF ALTITUDES. 

C 

DO 100 L = l,J-t 

TYPE 70, SIGNEO(L), TIF(L,2> 

ACCEPT 90, HEI 
100 ALT (L ,2 ) : HEI 

C 

C TYPE QUESTI ON CONCERNING THE NEED FOR 
C ANY CORRECTIONS. 

C 

110 Type 120 

120 FORMA T ( * ARE THERE ANY CORRECTIONS. Y OR N ’,$) 

ACCEPT 130, ANS 
130 FORMA T(A2) 

I FCANS ,EQ . 'N ' ) GO TO 190 

I F (ANS . NE . *Y ' ) GO TO 110 

IF(J.GT.I) K=J 
I F(J .GT.I ) GO TO 135 
K:I 
C 

C REQUESTS INPUT OF INCORRECT ALTITUDE 
C THEN SEARCH DATA ARRAYS FOR THE - 
C INCORRECT ALTITUDE. 

C 

135 TYPE 140 

140 FORMAT ( ' INCORRECT ALTITUDE: ’,3) 


151 


ACCEPT 90, ALTIN 
DO 1*0 L s 1 ,K“ 1 

1F(ALTIN.GT.ALT(L,1>-0.1.AND.ALTIN.LT.ALT(L,1)+0.1) GO TO 170 
150 I F CALTI N .GT .ALT (L, 2)~0. 1. AND. ALTIN. LT. ALT (L , 2)+0 . 1 ) GO TO 180 
TYPE 160 

160 FORMAT ( ’ - ERROR: INCORRECT ALTITUDE NOT FOUND 

!---»*,////) 

GO TO 110 
C 

C SECTION TO TYPE OUT CONDUCTIVITIES (POS. OR 
C NEG.) AND THE TIME CORRESPONDING TO THE 
C INCORRECT ALTITUDE, THEN WAIT FOR THE 
C INPUT OF THE CORRECT VALUE. 

C 

170 TYPE 185 

TYPE 70, SlOPOSa), 73Ffl.,i> 

ACCEPT 90, HE1 
ALT (L, 1 >=HEI 
GO TO 113 
180 TYPE 185 

185 FORMAT ( ' CONDUCTI VI TY * , AX , ‘TIME (IN SEC)’, 

1 AX , 'CORRECT ALTITUDE:???? ' ,//) 

TYPE 70, SIGNEG(L), TIFCL.2) 

ACCEPT 90, HE1 
ALT (L ,2) : HEI 
GO TO 110 
C 

C OUTPUT OF CONDUCTIVITIES, TIMES, AND ALTITUDES 
C TO THE STORAGE DEVICE. FIRST POSITIVE THEN 
C NEGATIVE. 

C 


190 

CALL OFILEUPSI, *PSI ') 
DO 200 L : 1 ,1 ~ 1 


200 

WRITEUPSI ,213) SIGPOS(L), TIF(L,I>, 

ALT (L , 1 ) 

210 

FORMAT(E10.5,F10.5,F10.5> 
CALL RELEAS ( IPSI ) 

CALL OFILEUNSI , 'NSI ') 

DO 220 L = l,J-l 


220 

WRI TECINSI ,21 1 ) SIGNEG(L), TIF(L,2), 

T I F (L , 3 ) , ALT(L ,2) 

2 11 

FORMAT (E10.5,F10.5,F10.5,F10.5) 


C 

CALL RELEAS (1NSI ) 



C PRINT OUT OF THE HEADING PERTAINING TO IKE 
C POSITIVE CONDUCTIVITIES AND THEN PRINT OUT 
C THE CONDUCTIVITIES. 

C 

PRINT 230 

230 FORMAT ( '1 \ 'POSITIVE C ONDUCTI VI TES WITH CORRESPONDING 

1 ALTITUDES*,//,’ *, 5X , ’CONDUCTI VI TY * ,20X , *ALTI TUDE * , //, 

2 ’+ * , 5X » '(IN MHOS/CM) ’,20X, '(IN KM)’,//) 
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00 240 

240 PRINT 250, SIGPOSCL), ALT (L , 1 ) 

250 FORMAT ( * ’ , 6X ,E 1 0. 5 ,20X , F 10 . 5 > 

C 

C PRINT OUT OF THE HEADING PERTAINING TO THE 
C NEGATIVE CONDUCTIVITIES AND THEN PRINT OUT 
C THE CONDUCTIVITIES. 

C 

PRINT 260 

260 FORMATt ’NEGATIVE CONDUCTIVITIES WITH CORRESPONDING 

1 ALTITUDES AND ASSOCIATED E FIELD’,//,* *,5X 

2, ’CONDUCTIVITY ’.20X, ’ALTITUDE ’, //, ’■*■’, 5X, ’(IN MHOS/CM) ’, 
320X, ’(IN KM) ’,2 IX, ’(IN V/CM) ’,//) 

DO 270 L = 1,J-1 

270 PRINT 275, SIGNEG(L), ALT (L ,2) , TIF(L,3) 

275 FORMA T ( * ’, 6X ,E 1 0 . 5 ,20X ,F 10. 5 , 1 8X ,F 10.5 ) 

IF(IDSK .EQ. 1 ) GO TO 2B0 
INSI =INSI-8 
IPSI =IPSI-B 
C 

C CHAIN BACK TO LINK ’’MAIN’’. 

C 

280 CALL CHAIN(0,IOSK, 'MAIN ') 

END 
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C ******************************************************* 

c 

C SIGPl PLOT ALTITUDE VS. 

C CONDUCTIVITY 

C 

C ********************************************************* 

COMMON LP »ISHOW,XMAX ,XMI N ,YMAX , YMI N , I NTENS.ISCALE 
COMMON /A/ ADVAL (3200 ) , GRIDU000) /D/STAR(20), SECO0) 

1 /H/ TIF(153,3> /K/ SIGPOS(I50), SIGN£G(15B> 

2 /Q/ ALT( 150,2) /R/ ASK 1(2) 

3 /DA/ 1DSK /DB/ NUM2, LB /N/ IPSI.INSI /U/ I, J 
C 

C STATEMENT REQUIRING OPERA TIR ACTION UPON 

C ENTERING THIS LINK. 

C 

PAUSE ‘ALUTIQE VS- COtfflTICTI VI TY PLOT* 

1=1 

J=1 

C 

C READ IN CONDUCTIVITIES £ CORRESPONDING 
C ALTITUDES. IF OPERATION IS OFF THE DISK 
C "GO TO 105"; OTHERWISE, PERFORM TWO 
C ADDITIONS SHOWN. 

C 

c 

IF(IDSK.EG.l) GO TO 105 
IPSI =IPSI+8 
INSI =INSI+8 
C 

C INITIALIZE STORAGE DEVICE FOR INPUT OF 
C POSITIVE CONDUCTIVITY FILE AND INPUT 
C THE INFORMATION INTO THE TWO ARRAYS 
C SHOWN. 

C 

105 CALL I FI LE (IPSI , 'PSI ’) 

110 READdPSl , 120,END = 130) SI GPOS (I ) ,TI F (I , l ) ,ALI (1,1) 

120 FORMA T(E10.5,2F10.5) 

I =1 + 1 
GO TO 110 

130 CALL RELEAS (IPSI ) 

C 

C INITIALIZATION OF STORAGE DEVICE FOR INPUT 
C OF NEGATIVE CONDUCIITITY FILE AND INPUT T«E 
C INFORMATION INTO THE TWO ARRAYS SHOWN. 

C 

CALL IFILEUNSI , 'NSI *) 

140 READ (I NS 1 , 1 45 ,END = 1 50 ) SIGNEG(J), TIF(J,2), TIF(J,3), ALT(J,2) 

145 FORMAT(E10.5,3F10.5) 

J=J + 1 
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GO TO 140 

150 CALL RELEAS (INS I ) 

EXPI=10.0**-14 

EXP2=10.0**-9 

EXP3=10.0**-10 

EXP4=10.0**-I2 

C 

C . INITIALIZE ARRAY GRID 

C 

I LP = 1 

1SCALE :0 

CALL I NTAB (GRID , 1000 ) 

INTENSsl 

ISHOWzl 

ADV=300.0 

C 

C GENERATION OF HORIZ . LINES FOR 10 INCREMENTS 

C IN ALTITUDE 

C 

DO 20 N= 1 , 7 

CALL POINT(GRID,215.0,ADV> 

CALL VECONT ( GRID, 215. 0,ADV ( 1 023.0, AD V) 

20 ADV=ADV+ 117.0 

INTENS =5 
C 

C INSERTION OF INTERMEDIATE HORIZ. LINES 
C DENOTING 5KM INCREMENTS IN ALTITUDE 
C 

ADV=358.5 
DO 50 N=l,S 

CALL FOINT(GRID,215.0,ADV) 

CALL VECONT ( GRID, 215. 0,ADV, 1023.0 f ADV) 

50 ADV=ADV+ 117.0 

INTENS = 7 
C 

C GENERATION OF ALTITUDE SCALE FACTORS 
C 

YCOOR =293.0 
N =25 

DO 100 L=1 • 13 
ENCODE (2 , 60 ,ASK 1 ) , N 
ASK 1 (2) =0 ,0 

CALL CHRGEN (GRID, ASK 1 , 165.0, YCOOR ,3,0) 

YCOOR =YCOOR + 58.5 

100 N = Wf5_ , _ _ __ _ - _ . _ 

C 

C GENERATION OF LABLES . 

C 

CALL CHRGENIGRID, 'Z ’, 100 . 0 , 750 . 0 , 5 , 2 ) 

CALL CHRGEN(GR1D, '(IN KM )', 70.0, 71 5. 0,2,0) 
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CALL CHRGENCGRID, 'CONDUCTI VI TY ',340. 0,230.0, 3,0) 
CALL CHRGENCGRID, ‘(IN MHOS/CM) ’,5513.13,230.0,2,0) 
CALL CHRGENCGRID, ’NEXT PHASE \ 10.0,50.0,2,0) 

C 

C DECISION SECTION 
C 

IFCNUM2.NE.i0O) GO TO 970 
I Nl =5 
NI =14 

ADVI=t96.3 
GO TO 990 

970 IFCNUM2 . NE . 1 ) GO TO 9S0 
I N 1 =4 
NI=12 

ADV1 =262. 0 
GO 70 990 
980 INV=6 

NI =14 

ADVI=150.0 

C 

C SECTION TO PLOT POSITIVE OR NEGATIVE CONDUCTIVITIES 
C OR BOTH. 

C 

C SET UP LOG SCALEi MAJOR VERTICLE LINES 
C 

990 ADV=21 5 . 0 

DO 1200 N= 1 ,1 N 1 

CALL POINTCGRID.ADV, 300.0) 

CALL VECONT CGRID, ADV, 300. 0,ADV, 1023.0) 

1000 ADV=ADV+ADVI 

C 

C GENERATION OF MINOR VERTICLE LINES 
C 

IN TENS =6 
AL=2 .0 

ADV=215.0+ALOG10CAL)*ADVI 
1005 DO 1010 N : 1 ,1 N t - 1 

CALL P 01 NT CGRID, ADV, 3 J0.0) 

CALL VECONT CGRID, ADV, 320.0, ADV, 1023.0) 

1010 ADV=ADV+ADV1 

AL=AL+ 1 .0 

ADV =2 1 5 ,0+ALOG 10(AL)*ADVI 
IFCAL.LE .9.0) GO TO 1005 
INTENS =5 

ADV =2 15 . 0+ALOG 10CI . 5 )*ADVI 
IFCAL.LE. 10.0) GO TO 1005 
C 

C GENERATION OF CONDUCTIVITY SCALE FACTORS 
C 


INTENS =7 
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XCOOR =168.0 
ENCODE (2, 60, ASK 1 ) g LB 
60 FORMAT (12) 

ASK 1 (2)=0.0 
00 l 020 N=1,IN1 

CALL CHRG£N(GRID,ASK 1 ,XC00R, 265. 0,3,0) 

1020 XCOOR =XC0CR+ADVl 

XCOOR =203.0 
DO 1030 Nsl, INI 

CALL CHRGENCGRID, XCOOR, 277. 0,2,0) 

1030 XCOOR =XC00R+ADV1 

XCOOR =210.0 
DO 1340 L = 1 « I K 1 
ENC0DE(2,63,ASK1), K1 
ASK 1 (2)=0.0 

CALL CHB GEN (GRID, ASK 1 .XCOOR, 2 77, 0,2,0) 

XCOOh =XCD0JH-A2V1 
1040 NI=NI-1 

I NTENS =7 

INITIALIZATION OF DATA DISPLAY ARRAY "ADVAL". 

CALL INTAB (ADVAL, 3000) 

GO TO (1200, 1300), KUK2 
NUM2: 1 

LOCATING REGION, ON DISPLAY, IN WHICH A 
GIVEN POSITIVE CONDUCTIVITY SHOULD BE 
PLACED. UPON FINDING THE REGION "GO TO 1050". 

IF CONDUCTIVITY DOES NOT FALL IN ANY REGION 
TO BE DISPIAYED "TYPE 60" AND CONTINUE. 

DO 1080 1=1, I-l 

IF(SIGP0S(L).GT.EXP1.AND.SIGP0S(L).LT.EXP3) GO TO 1050 
TYPE 160, SIGPOS(L), ALT(L,l) 

GO TO 1080 
1 XSTART =215.0 

IEXP ~~ 13 

SECTION TO PERFORM NECESSARY OPERATIONS ON 
CONDUCTIVITIES IN ORDER TO PLOT DATA AND 
CONSTRUCT DATA ARRAY "ADVAL". 

DO 1060 K = l . 4 

IF(SIGPOS(L) .LT. !0.0**IEXP ) GO TO 1070 
' XSTART =XSTART+196T3 ' 

1060 IEXP =1 EXP+ 1 

1070 XC =ALOG 1 (SIGPOS(L)*l0.0**-ISXP)*10.0)*l96.3+XSTART-4.0 

YC:294.0+(ALT(L,l)-25.0)*l 1 .7 
CALL CHR GEN (ADVAL , ’+ ' ,XC ,YC ,2,0 ) 
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1080 CONTINUE 
C 

C DISPLAY DATA AND WAIT FOR A "NEXT PHASE" 

C DETECTION , UPON A LITE PEN DETECTION 
C "CO TO 1" A(VD PREPARE FOR PLOT 
C OF NEGATIVE CONDUCTIVITIES. 

C 

CALL DISPLY (2, GRID , ADVAL) 

1090 CALL LITEPN(IACC,X5,Y5> 

IFUACC.EQ.0) GO TO 1093 
CALL WAI T( 500 > 

CALL L1TEPN<ICLE,X4,Y4> 

GO TO l 
1230 NUK2s2 
C 

C LOCATING REGION, ON -DISPLAY, IN WHICH A 
C GIVEN NEGATIVE CONDUCT! VI TY SHOULD BE 
C PLACED. UPON FINDING IHE REGION "GO TO 1210". 

C IF CONDUCTIVITY DOES NOT FALL IN ANY REGION 
C TO BE DISPLAYED' “"TYPE 63" AND CONTINUE. 

C 

DO 1240 L = 1 ,J- 1 

1F(S1GNEG(L) .GT ,EXPA .AND .SIGNEG (L> .LT.EXPZ) GO TO 1210 
TYPE 160, SIGNEG(L), ALT CL ,2) 

GO TO 1240 
1210 XSTART =21 5.0 
IEXP=-11 
C 

C SECTION TO PERFORM NECESSARY OPERATIONS ON 
C CONDUCTIVITIES IN ORDER TO PLOT DATA AND 
C CONSTRUCT DATA ARRAY "ADVAL" . 

C 

DO 1220 K = 1 ,3 

I F(SI GNEG (L) ,LT .1 3.0**I£XP ) GO TO 1230 
XSTART :XSTART+262.0 
1220 IEXP=IEXP+l 

1230 XC =A LOG 1 0 ( (SI GNEG (L )* 1 0 . 0***1 EXP )* 10.0) *262. 0+XSTART-4.0 

YC=294.0+(ALT(L,2)-25.0>*11.7 
CALL CHRGEN (ADVAL , ',XC ,YC ,2,0) 

1240 CONTINUE 
C 

C DISPLAY DATA AND WAIT FOR A "NEXT PHASE" 

C DETECTION. UPON A LITE PEN DETECTION 
C "GO TO 1” AND PREPARE- FOR PLOT 
C OF BOTH POSITIVE AND NEGAIIVE CONDUCTIVITIES. 

C 

CALL 0ISPLYC2, GRID, ADVAL) 

CALL LITEPN(IACC,X5,Y5) 

IF (IACC . EQ .0 ) GO TO 1250 
CALL WAITC5Z3) 


1250 
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CALL LITEPNCICLE ,X4,Y4) 

GO TO 1 
C 

C SECTION TO PLOT BOTH POSITIVE AND NEGATIVE 
C CONDUCTIVITIES TOGETHER. 

C 

c 

C PLOT ROUTINE FOR POSITIVE CONDUCTIVITY. 

C 

C 

C LOCATING REGION, ON DISPLAY, IN WHICH A 
C GIVEN POSITIVE CONDUCTIVITY SHOULD BE 
C PLACED. UPON FINDING THE REGION "GO TO 170". 

C IF CONDUCTIVITY DOES NOT FALL IN ANY REGION 
C TO BE DISPLAYED "TYPE 60" AND CONTINUE. 

C 

1300 DO 200 L=l,l-1 

IFCSIGPOS CL) .GT .EXP I .AND .SI GP OS (L ) .LT . EXP2) GO TO 170 

TYPE 160, SI GPQSXL ), ALT CL , 1 ) 

160 FORMAT ( * *,*<--- ERROR: VALUE IS OUT OF DISPLAY 

l SCOPE RANGE *,5X, , C0NDr *,E10.5,5X, 

2'ALTITUDE: ',F10.5,//, * VALUE WAS IGNORED, 

3 PROCESSING BEING CONTINUED*,///) 

GO TO 200 

170 XSTART:215.0 

IEXP s- 13 
C 

C SECTION TO PERFORM NECESSARY OPERATIONS ON 
C CONDUCTIVITIES IN ORDER TO PLOT DATA AND 
C CONSTRUCT DATA ARRAY "ADVAL" . 

C 

DO 180 N: 1 ,5 

IF(SIGPOS(L).LT.10.0**IEXP> GO TO 190 
XSTART =XSTART+ 150.0 
180 IEXP=IEXP+1 

190 XC=ALOG10(<SIGPOS(L)*l0.0**-IEXP>*l0.0>*l50.0+XSTART-4.0 

YC =2 9 4.0+ (ALT (L,l)-25.0)*ll.7 
CALL CHRGEN (ADVAL , *+ *,XC ,YC ,2,0) 

200 CONTINUE 

C 

C PLOT ROUTINE FOR NEGATIVE CONDUCTIVITIES. 

C 

C 

C LOCATING REGION, ON DISPLAY, IN WHICH A 
C GIVEN POSITIVE CONDUCTIVITY SHOULD BE 
C PLACED. UPON FINDING THE REGION J’GO^TO 230". 

C IF CONDUCTIVITY DOES NOT FALL IN ANY REGION 
C TO BE DISPLAYED "TYPE 60" AND CONTINUE. 

C 


DO 260 L=I,J-1 
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IF<SIGN£G(L).GT.£XPI.AND.SIGNEG(L>.LT.EXP2) GO TO 230 
TYPE 160, SIGNEG(L), ALTCL.2) 

GO TO 260 

230 XSTART=215.0 
IEXPs-13 
C 

C SECTION TO PERFORM NECESSARY OPERATI ONS ON 
C CONDUCTIVITIES IN ORDER TO PLOT DATA AND 
C CONSTRUCT DATA ARRAY "ADVAL". 

C 

DO 240 N= i , 5 

IF(SIGNEG(L).LT.10.0**IEXP) GO TO 250 
XSTART=XSTART+150.0 
240 1EXP=IEXP + 1 

250 XC=ALOG1MYSm’EGYL)*l0.0**-IExP)*i0.0)*150.0+XSTART-4.0 

YC =294. 0+ tALT a^.>-25^ 0 >♦ 1 
CALL CWGEtKADVAL,*-*,XC,YC,2,0) 

260 CONTINUE 

C 

C DISPLAY OF DATA ARRAY '"ADVAL" AND WAIT 
C FOR "NEXT PHASE" DETECTION. UPON A LITE 
C PEN DETECTION TERMINATE THE PROGRAM. 

C 

CALL DISPLY<2, GRID, ADVAL) 

267 CALL LITEPNdACC ,X5,Y5) 

IF (IACC .EQ .0) GO TO 267 

IF(X5.GT. 10.0. AND .X5.LT . 110.0.AND .Y5.GT .50.0. AND . 

IY5.LT. 65.0) GO TO 264 
GO TO 267 

264 CALL VAI T (500) 

CALL LITEPN(ICL£,X4,Y4) 

C 

C TYPE OUT TERMINATION STATEMENT AND 
C END THE PROGRAM. 

C 

TYPE 280 

280 FORMAT C * PROGRAM COMPLETED A TERMINATED . \ /////) 

END 
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c ***************************************************** 
c 

C PROGRAM TO CALCULATE DENSITIES 

C 

C *************** ************************************** 

DIMENSION ALTMP (651 ,3 ) , COND(150,8), DEN(150,8), EOPAB( 2000,3) 

J=l 

1=1 

C 

C TYPE OUT OF QUESTION REQUESTING THE 

C NAME OF THE LAUNCH FACILITY. FOLLOWING 

C THAT ARE THE STATEMENTS WHICH DETERMINE 

C WHAT DATA FILE SHOULD EE READ, THIS 

C DEPENDS ON THE RESPONCE TO THE QUESTION 

C ASKED. 

C 

1 TYPE 2 

2 FORMAT C LAUNCH SITE, WSKR OR W1 .-..*,$) 

ACCEPT 3, WHR 

3 FORMAT (A A) 

I F (WHR .EQ. *WI ’) GO TO 4 
IF(WHR .NE . 'WSMR * ) GO TO 1 
CALL I FILE ( 1 , *ZTP ' ) 

GO TO 5 

4 CALL I FILE ( 1 , 'WIZ TP *) 

5 READd, 10,END=20) ALTMP(I,1>, ALTMP(I,2>, ALTMP(I,3> 

10 FORMAT(2F10.5,E10.5) 

1=1 + 1 
GO TO 5 

20 CALL RELEAS(l) 

C 

C INPUT OF FILE CONTAINING 

C POSITIVE CONDUCTIVITY VALUES. 

C 

CALL I FILE ( I f PSI # ) 

25 READ( l ,30,END = 4U) COND(J,i>, TIF, CQND(J,2) 

30 FORMAT(E10.5,2F10.5) 

J=J+1 
GO TO 25 

40 CALL RELEAS ( 1 ) 

I =1-1 
J=J-1 

CONST =293. 0*10.0** 13/(1 .8*1.01*1.6) 

K = 1 

C . , * - 

C DO LOOP WHICH (!) MATCHES ALTITUDES 
C OF CALCULATED CONDUCTIVITIES TO ALTITUDES IN 
C TABLE WHICH CONTAINS TEMPERATURE 
C S, PRESSURE DATA, AND ( 2 ) CALCULATES THE POSITIVE 
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C ION DENSITY. 

C 

45 DO SO N = I,I 

IF<<COND<K ,2).LT.ALTMP<N,l)-.3t).0R.<C0ND<K,2).GT. 
IALTMP<N,1)+.0I)) GO TO 50 

DEN<K,l)=C0ND<K,l)*C0NST*<ALTMP<N,3)/ALTMP<N,2)) 
DEN<K,2)= C0ND<K,2) 

GO TO 55 
50 CONTINUE 

55 K=K+1 

I F <K .LE mJ ) GO TO 45 
C 

C OUTPUT POSITIVE ION DENSITY TO DISK UNDER THE 
C FILE NAME PDEN.DAT AND THEN PRINT OUT DENSITIES 
C 

CALL OFILE < 1 , 'PDEN ’ ) 

DO 60 K = 1 ,J 

60 WRI T£ < 1 , 70 ) DEN<K,1), DENCK ,2) 

70 FORMAT<E10.5,F 10.5) 

CALL RELEAS < 1 ) 

PRINT 80 

80 FORMA T< *1 ',5X, ‘THE POSITIVE DENSITIES FOLLOW* ',///, 

I* ', 5X « "DENSI TY * ,20X , ’ALTITUDE*,//, *+ *,2X, * <+ ION/CM 
2**3) ’ , 18X, *<IN KM) *,//) 

DO 90 K = 1,J 

90 PRINT 100, DEN<K, 1 ) , DEN <K ,2) 

100 FORMAT < * \2X,E10.5,20X,F10.5) 

Ksl 

C 

C INPUT OF FILE CONTAINING NEGATIVE 
C CONDUCTIVITIES. 

C 

CALL I FILE < 1 , ’NSI ') 

110 READ< 1 , 120, END =150) C0ND<K,3), TIF, C0ND<K,5), C0ND<K,4) 
120 FORMAT <£10.5,3F10.5) 

K :K+1 
GO TO 110 
130 K =K - 1 

L = 1 
C 

C LOOP TO CALCULATE E/P 
C 

136 DO 140 N:l,I 

I F< <COIJD <L , 4) .LT .ALTMP <N , 1 )- .0 1 ) .OR . <COND <L , 4) ,GT . 
1ALTMP<N,1)+0.0D) GO TO 143 

CONO <L, 6) =<CONO<L, 5)* l .01*10.0**4) /<ALTMP <N,3)*7.6) 

L=L+l 

1 F <L.LE.K ) GO TO 136 
IF<L.GT.K) GO TO 150 
140 CONTINUE 
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c 

C INPUT OF FILE CONTAINING VALUES OF 
C ALPHA 4 BETA FOR GIVEN VALUES OF EIP. 

C 

150 CALL IFILEU, ’ALBA ’) 

1 EOP = I 

160 READ ( 1 , 1 70, END = 1 80) EOPAB (IEOP , 1 ) , EOPAB (IEOP ,2) , E0PAB(IE0P,3> 

170 FORKAT(F10.5,2E12.5> 

lEOPsIEOP + 1 
GO TO 160 
C 

C SECTION TO DETERMINE UE 
C 

180 CALL RELEAS(l) 

l£OP=I£OP-l 
LA =1 

185 DO 190 I A = 1 ( 1E0P 

I F < (CCND (LA , 6) .LT .EOPAB ( I A , 1 >-0 . 005 > .OR . (COND (LA , 6) .GT . 

I EOPAB (IA f 1) +0.005)) GO TO 187 

D£N(LA , 3 ) = ( 1 /C OND (LA ,5))*((EOPAB(IA , 2)*E0PAB (IA , 1 ) J+EOPAB (I A ,3 ) ) 
LA =LA+l 

I F (LA .GT.K ) GO TO 200 
GO TO 185 

187 IF (IA .LT .IEOP ) GO TO 190 

TYPE 186, COND (LA, 6) 

186 FORMA T ( ' <---ERRORj THE E/P 

1 RATIO OF ’.Fta.S,* WAS OUT*,/, 

2’ OF THE RANGE OF THE TABLES 
3 — -»*,///////> 

LA =LA+ 1 

IF (LA .GI.K J GO TO 200 
GO TO 185 
190 CONTINUE 

C 

C CALCULATIONS OF U- 
C 

200 LA:1 

210 DO 220 I A =1 , 1 

IF((C0ND(LA,4).LT.ALTMP(IA, 1 ) -0.01). OR. (COND (LA, 4). GT. 

1ALTMP (IA, 1)4-0.01)) GO TO 220 

DEN(LA,4):(2.3*I ,0 1* 1 0.0**6*ALTMP (I A ,2))/(293.0*ALTMP (IA ,3} ) 

LA sLA + 1 

IF (LA .GT .K ) GO TO 230 
GO TO 210 
220 CONTINUE 

230 LA~=P " 

I NC =J - 1 

231 ALT1NC=0.0 

DO 280 I A = I , I NC 

IF(COND(LA, A).LT.DEN'(1A+1,2>) GO TO 250 
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1NCALT=(DEN(IA,2)*1 0 .3-DEN (IA+1 ,2)*l0.0)+2 

00 240 IB = 1 ,1 NCALI 
CHX =DEN(1A ,2)-ALTINC 

1 F(C0N'D (LA , 4) .LT .CHK -3. 04.OR .COKD (LA , 4) ,GT« 
iCHK+0.04) GO TO 235 

FACTOR :(D£N(IA+l , I )-DEN (IA , 1 ) ) /FLOAT (INC ALT) 

FACTOR =DEN(IA,1 )+FAC TOR* FLOAT (IB) 

C 

C ELECTRON DENSITY CALCULATION 
C 

DEN(LA,5)=(I.0/DEN(LA,3))*((COND(LA,3>*10.0**19/1.6) 

1 -FACTOR*DEN (LA , 4 ) ) 

C 

C CALCULATION OF NEG ION DENSITY 
C 

DEN(1A« 6)=FACT0R— DEN (LA.5 > 

LA rLA+1 

IF(LA .GT.K ) GO TO 2S0 
GO TO 231 

235 ALTINC=ALTINC+0.1 

240 CONTINUE 

C 

C Type OUT OF ERROR MESSAGE WHEN A VALUE 
C WAS NOT ACCURATELY CONVERTABLE TO A 
C DENSITY. 

C 

TYPE 245 

245 FORMA T( *- «— ERROR IN PROCESSING, PROGRAM CONTINUING.— » ’,//) 

ALTINC=0.0 
GO TO 255 

250 IF(IA.LT.INC) GO TO 280 

255 TYPE 260, C0ND(LA,3>, C0ND(LA,4> 

260 FORMAK’ THE NEGATIVE CONDUCTIVITY OF ',E10.5,' AT AN 

1 ALTITUDE OF’,/, ’ ’.FIO.S,’ WAS NOT ACCURATELY CONVERTABLE TO A 

2 DENSITY; ’ THEREFORE, IT WAS IGNORED. ’//////////> 
DEN(LA,5)=0.0 

DEN (LA , 6) = 0.0 
LA ;LA + 1 

IF(LA.GT.K) GO TO 290 
GO TO 231 
280 CONTINUE 

C 

C OUTPUT OF DATA TO DISK UNDER THE 
C FILE NAMED NENI .DAT. 

C 

290 CALL OFILEd, ’NENI ’) 

DO 300 LA =1 ,K 

300 WRITE (1,310) DEN(LA ,3 ) , DEN (LA , 4 ) , DEN (LA , 5) , DEN(LA,6) 

1, COND(LA , 4) 

FORMAT (4E 10. 4, F 10.5) 


310 



oooo 
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CALL RELEAS(l) 

PRINT OUT HEADING SHOWN BELOW 4 CORRESPONDING 
DATA. 


PRINT 320 

320 FORMA T ( ' 1 MOBILITIES AND NEGATIVE ION DENSITIES FOLLOW 

1»\//, * \2x, 'ELECTRON MOBILITIES ",5X, ’NEG ION MOBILITIES’, 
25X , 'ELEC DEN’,I0X,’- ION DEN 1 0X , ’ALII TUDE 
3 AX , '(IN CM**2/VS) ',10X, '(IN CM**2/VS > ' ,6X , ' (ELEC/CM**3) ' , 
A7X, '(I0N/CM**3) *;10X, '(IN KM) ',//) 

DO 330 LA=1,K 

330 PRINT 3 A0, DEN (LA , 3 ) , DEN (LA , 4) , DEN (LA , 5) ,DEN (LA , 6) , 

1 COND(LA.A) 

340 FORMAT ( ' \6X,EI0.5,l4X,E10.5,7X,E10.4,llX,El0.4 f 8X,FI0.5> 
STOP 
END 
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c *************************************************** 
c 

C PROGRAM TO LINEARLY INTERPOLATE TEMPERATURE 

C AND PRESSURE DATA FROM CIRA 1965. 

C 

C *************************************************** 

DIMENSION ALTMP (651,3), FIND(2,3) 

Mil 
L:50 
K=0 
TYPE 2 

2 FORMA T ( ' FOR WSKR USE LATITUDE OF 30N 4 FOR 
l WI USE 40N*,/,* TYPE WSMR OR WI *,$) 

ACCEPT 3, WHR 

3 FORMAKA4) 

I =1 

5 TYPE 10 

10 FORMAT ( * ALTITUDES *,$> 

ACCEPT 30, ALT 

30 FORMAT (F 10.5) 

TYPE 20 

20 FORMAT ( * TEMPERATURE (INDEGK)s ’,%) 

ACCEPT 30, TEMP 
TYPE 35 

35 FORMAT! * PRESSURE (AS 0.000E+00, IN DYNES/CM**2)s ' ,i) 

ACCEPT 40, PRES 

40 FORMAT (E10.5) 

FIND (I , 1 )s ALT 
F1NDO ,2): TEMP 
FIND<l,3)s PRES 
I F C 1 . G£ . 2 ) GO TO 50 
1=1 + 1 
GO TO 5 

50 ALTMP (M, l ) : F1ND(I-1,1> 

ALTMP (M ,2 ) : FINDCI-1,2) 

ALTMP (M ,3 ) : FIND(I-t,3> 

ALT1NC: (FIND (I , 1 >-FI ND (I -1 , 1 > ) /50.0 
TMPINC: (FIND (I .2J-FIND (I -1 ,2) 1/50.0 
PREINC: (FIND (1 , 31-FIND Cl -1,3)) /5B.0 
DO 60 N:M,L 

ALTMP (N+ 1,1): ALTMP (N , 1 )+ALTI NC 
ALTMPCN+l ,2): ALTMP (N ,2 )+TMPI NC 
A L IMP ( N+ i , 3 > : AL W C .« ; -3 ) +P RE I NC 
I F (ALTMP (N, D.GT.B0.0) ALTMP (N+ 1 ,2) : 1 99.0 

IF (ALTMP (N , 1 ) .GT .80.0 .AND .WHR .EQ . *WI ') 

1 ALTMPCN+l, 2):194.0 

60 IF(ALTMP(N, D.GT.80.0) ALTMP (N+ l ,3 ) : . I 1400E+0 

12*EXP(-( (ALTMP (N, 1 >-80.0>/5.59>> 

I F (ALTMP (N, D.GT. 80.0. AND. WHR. EQ . 'WI ') 
i ALTMP (N+ 1,3):. I 1300E+ 02+EXP (- ( (ALTMP (N , l )-80.0>/5. 59)) 
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FIND<t,l)=Fltn<2,t) 

FINDCl,2)=FIN0<2,2) 

FI ND(1 ,3 )=FIND(2,3 > 

n =M+50 

L=L+50 

IFtM.LE.650) 00 TO 5 
IFCVIHR.EQ.'Wl ') GO TO 85 
CALL OFILEU, # ZTP ’) 

65 DO 70 1=1,651 

70 WRI TEC 1 ,80) ALTMP(I.I), ALTMP<I,2), ALIMP<I,3) 

80 FORMAT(2F10.5,E 10.5) 

CALL RELEAS(l) 

GO TO 89 

85 CALL OFILEU, ’WIZTP ’) 

GO TO 65 

89 PRINT 90 

90 FORM T ( T ', ’ ALTITUDE \I0X, TEMPERATURE \10X, 'PRESSURE', 
l// t •+*,1X, *<IN KM) *,11X, ’(IN DEG 

2 (I N DYNES /CM*#2) ',///) 

DO 130 1 =1,651 

100 PRINT 110, ALTMP (1,1), ALTMP (1,2), ALTMP(I,3> 

110 FORKATC • t F10.5,8X,F10.5,8X,El0.5) 

STOP 

END 
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C *******************************.<(******************** 

c 

C PROGRAM TO DETERMINE ALPHA AND BETA FOR KNOWN 

C VALUES OF E/P. 

C 

C ******** ******♦**>********************* ******** ********* 

DIMENSION EOPAB (2300,3) 

M=2 
N = 5 
1=1 
C 

C STATEMENTS REQUESTING INPUT OF DATA 

C FROM TABLES. 

C 

5 TYPE 13 

10 FORMAT! ' E/P = \$> 

ACCEPT 15, EOPAB (1,1) 

15 FORMAT (F 1 0 . 5 ) 

TYPE 20 

20 FORMAT! * ALPHA !IN THE FORM ,00000E+00>= *,$> 

ACCEPT 30, EOPAB (I, 2) 

30 FORMAT !E 1 0. 5) 

TYPE 40 

40 FORMAT! ' BETA !IN THE FORM .00000E+00) = ’,*> 

ACCEPT 30, EOPAB (I ,3 ) 

IFU .GE.2) GO TO 45 
I =1 + 1 
GO TO 5 
C 

C DETERMINATION OF INCREMENTAL VALUES. 

C 

45 CEOP = ! EOPAB ! 1 , 1 ) -EOPAB !I-1,I))/0.01 

AI NC =! EOPAB !1 ,2)-E0PABU - 1 ,2) ) /CEOP 
BI NC =! EOPAB !1, 3) -EOP ABU-1,3)) /CEOP 
C 

C DO LOOP TO INTERPOLATE VALUES AND INSERT 
C THEM IN THE ARRAY EOPAB. 

C 

DO 50 K=M,N 

EOPAB !K,l)=EOPAB(K-l,l)+0. 01 
EOPAB !K ,2)=E0PAB(K- l ,2>+AINC 
50 EOPAB!K,3)=EOPAB(K-1,3)+BINC 

I =N+ 1 
M =N+ 1 

IF!N.EQ.5.0R .N.EQ . 10. OR .N.EQ .15) N=N+5 

IFIM. EQ.I6) GO TO 5 

IFIN. GT .15. AND .N.LE.200) N=N+20 
IFIN.GT.200) N=N+200 

I F IM. EQ .20 1 ) N=N-20 

IFIN. GT. 2000) GO TO 60 
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GO TO 5 
C 

C OUTPUT OF DATA INTO FILE NAMED ALBA .DAT . 

C 

62 CALL OFILECI, ’ALBA *) 

DO 65 L=l ,2030 

65 WRI TE ( 1 , 67) E0PA8CL,1>, E0PAB(L,2>, EOPAB(L,3) 

67 FORMAT(F10.5,2X,E 10.5, 2x,E 10.5) 

CALL RELEAS(l) 

C 

C PRINT OUT OF CALCULATED VALUES. 

C 

PRINT 70 

70 FORMAT! M TABLE TO FIND ALPHA £ BETA WHEN E/P IS 

1 KNOWN: *,/, '-',6X, 'E/P ',2 0X, 'ALPHA *,20X, 'BETA ',//, 
2'+V<IN V/CM-MMHG) *,SX, '(IN CM**2-MMHG/VS ) ' , 10X, 

3 '(IN CM/S) ',//) 

DO 80 L: 1 ,2000 

80 PRINT 90, EOPAB (L , 1 ) , EOPAB(L,2), E0PAB(L,3) 

90 FORMAT ( ' *,3X,F10.5,14X,E10.5,15X,E10.5) 

STOP 

END 



